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Abstract 


Every  mobile  robot  needs  to  know  where  it  is  so  everyone  who  does  mobile 
robotics  wants  to  know  how  inertial  and  satellite  navigation  works.  These  two 
technologies  are  the  most  sophisticated  forms  of  navigation  systems  available 
and  they  complement  each  other  very  well.  Applications  requiring  indications  of 
highly  dynamic  3D  motion,  excellent  relative  accuracy  and  high  update  rates 
benefit  from  inertial  systems.  Applications  requiring  bounded  absolute  accuracy 
for  extended  excursions  or  position  estimates  relative  to  the  earth  itself  benefit 
from  satellite  navigation. 

This  decade  has  seen  the  development  of  inexpensive  handheld  global  position¬ 
ing  systems  based  on  reception  of  the  GPS  satellite  signals,  and  advances  in  opti¬ 
cal  gyroscopes  and  integrated  circuit  accelerometers  which  promise  to  put  an 
inertial  navigation  system  in  every  automobile  before  the  end  of  the  next.  This 
report  is  a  detailed  tutorial  which  explains  the  principles,  practice,  and  issues  of 
using  these  new  technologies. 
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1.  Introduction 

This  century  has  seen  two  breakthroughs  in  navigation^  technology,  equal  perhaps  in  importance 
to  the  discovery  of  the  compass.  The  latest  is  the  ^obal  positioning  system,  or  GPS,  which  has 
recently  come  on  line.  Based  on  modified  triangulation  of  a  network  of  satellites  in  earth  orbit,  it 
allows  anyone,  anywhere,  to  fix  position  through  the  use  of  a  small  hand  held  receiver.  The  other 
breakthrough,  arriving  some  years  earlier,  is  inertial  navigation  systems  (INS).  These  systems, 
based  on  dead  reckoning^,  have  been  of  considerable  interest  to  the  military  and  aerospace 
industries  for  many  reasons.  This  report  will  discuss  the  basic  principles  underlying  the  operation 
of  inertial  and  satellite  navigation  systems  from  both  a  theoretical  and  a  practical  perspective. 

li.  Cttfflincntaiy 

One  of  the  distinguishing  features  of  the  study  of  navigation  is  that,  as  old  as  it  is,  it  seems  that  it 
has  always  been  the  province  of  a  small  group  of  specialists.  Most  other  ancient  technologies  like 
flint  tools,  fire,  the  wheel,  animal  husbandry  and  agriculture,  had  an  obvious  utility  to  the  ordinary 
person,  and  were  intellectually  accessible  to  the  masses.  Literacy,  while  once  enjoyed  only  by  the 
ruling  classes,  has  now  become  available  to  at  least  most  people  in  the  industrialized  world.  Of  all 
of  the  achievements  of  the  ancients,  navigation  seems  to  stand  alone  in  this  regard. 

This  project  started  as  an  attempt  to  understand  some  of  the  peculiarities  of  the  INS  as  it  is 
employed  on  mobile  robots  at  CMU.  In  particular,  the  issue  was  why  it  is  that  the  Litton  PADS 
inertial  navigation  system  then  in  use  on  the  NAVLAB  n  apparently  could  not  generate  a  reliable 
estimate  of  the  vehicle  z  coordinate.  This  has  become  a  very  important  question  because 
sophisticated  map  matching  algorithms  still  consume  too  much  time  for  high  speed  navigation.  The 
answer  to  this  question  turned  out  to  require  nothing  less  than  the  fundamental  postulate  of  general 
relativity. 

I  ended  up  doing  mote  work  than  I  intended,  and  came  to  answer  a  number  of  other  questions  which 
are  more  fundamental  prerequisites.  I  wrote  down  what  I  learned  in  the  belief  that  others  would 
find  the  information  useful.  I  had  a  hard  time  finding  any  teal  information  myself.  There  seems  to 
be  three  kinds  of  information  out  there,  industrial  product  brochures  which  teach  you  nothing, 
popular  pseudoeducational  settunars  which  sanitize  the  math  to  the  point  of  uselessness,  and  the 
academic  literature  which,  after  sixty  years  of  development  of  the  technology  behind  closed  doors, 
is  now  impenetrable.  What  was  intended  as  a  five  page  analysis,  has  become  twenty  times  that  size 
because  tte  more  1  looked,  the  more  I  did  not  understand. 

Ask  the  community  here  how  these  devices  function  and  you  are  likely  to  get  the  answer  “black 
magic”  as  often  as  something  techiucal.  Like  many  other  devices  whose  major  current  use  is  the 
military,  there  is  an  aura  of  secrecy  about  them.  The  lore  is  that  they  do  really  neat  stuff  like 
measuring  the  earth’s  rotation  and  that  tl^y  incorporate  those  seemingly  ultra  high-tech 
components  called  gyroscopes  which  are  synonymous  with  space  travel  and  decades  of  painstaking 
precision  engineering.  Because  they  use  the  gyroscope,  I  felt  I  could  do  this  project  for  academic 
credit  in  Matt  Mason’s  martipulation  course  and  get  a  double  benefit 


1.  Ihe  term  navigation,  wbidh  is  two  dimeasional,  is  distinguisbed  from  guidance,  vdiich  is  three  dimen- 
tional.  The  word  is  derived  from  the  latin  ”navis”  meaning  ship  and  “agere”  meaning  to  move  or  direct 

2.  Other  navigation  techniques  include  landmark  recognition  and  triangulation.  The  origin  of  the  term  is  “ded¬ 
uced  reckoning’’. 
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It  turns  out  that  the  devices  are  not  that  difficult  to  understand  at  aU  -  there  just  is  no  information 
available.  I  knew  there  was  some  deep  stuff  involved  because  I  had  heard  that  the  motion  of  the 
earth  and  the  gravitation  of  the  planets  had  to  be  modelled,  and  that  you  needed  general  relativity 
to  understand  the  new  optical  gyros,  but  this  is  overstated. 

In  the  course  of  the  work,  it  became  clear  that  integrated  INS-GPS  systems  are  the  wave  of  the 
future,  and  that  the  reason  for  this  is  the  long  recognized  symbiotic  relationship  between  dead 
reckoning  and  triangulation.  Therefore,  a  complete  effort  required  a  study  of  the  new  GPS  system 
as  well. 
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PART  I: Relevant  Concepts 
from  Theoretical  Mechanics 


In  order  to  present  the  implementation  of  a  working  INS,  some  aspects  of  physical  mechanics  that 
are  central  to  the  discussion  must  be  covered.  This  section  will  overview  ideas  from  the  branch  of 
physics  called  theoretical  mechanics,  rather  than  its  pragmatic  engineering  counterpart.  The 
reader  will  find  that  later  sections  will  refer  heavily  to  the  material  presented  here. 

1.  Context  of  Measurements  of  Physical  Quantities 

Physical  quantities  can  be  classified  into  tensors,  vectors  and  scalars,  which  are  all  familiar^. 
Consider  a  measurement  made  by  an  observer  of  a  velocity  of  a  train  relative  to  a  car  of  6.32  meters 
per  second  in  a  northeasterly  direction.  When  the  observer  announces  the  measurement,  there  is  an 
entire  context  in  which  it  is  made,  part  of  it  explicitly  stated,  part  understood  by  convention.  That 
context  involves  five  components: 

•  a  unit  system  (meters,  seconds) 

•  a  number  system  (6.32) 

•  a  coordinate  system  (north,  east  etc.) 

•  a  reference  frame  to  which  the  measurement  is  ascribed  (the  train) 

•  a  reference  frame  with  respect  to  which  the  measurement  is  made  (the  car) 


3.  Scalars  possess  nuignitmte  only,  vectors  m^nitode  in  a  partkiilar  diiection.  Ibosors  socb  as  solid  body 
stress  (not  fluid  pressure  -  wbicb  is  a  scalar)  aiHl  inertia  possess  magnitude  in  all  diiectioiis.  1bns(»s  possess 
tliree  principal  axes  of  magnitude,  from  wbich  all  others  can  be  determined  •  hence  tbe  use  of  the  inertia  ellip- 
soid  in  rigid  body  medianics  and  of  Mohr’s  circle  of  stress,  strain,  and  inertia  in  solid  mechanics. 
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UL  Unit  Systems 

The  unit  system  is  the  conventional  meaning  of  a  meter  and  a  second.  If  the  velocity  were  given  in 
fathoms  per  decade,  it  would  change  the  number  6.32  to  something  else,  but  it  would  still  reflect 
the  same  measurement 

L2.  Nwinfr^r  Sygt^ips 

The  number  system  is  the  conventional  meaning  6.32.  The  Arabic  number  system  is  a  base  ten 
wd^ted  positional  number  system.  Other  systems  include  the  Roman  system  which  is  not 
entirely  weighted  positional  and  has  no  regular  base.  If  the  number  6.32  were  given  in  Roman 
notation  (assuming  there  is  a  way  to  represent  decimals  at  all),  it  would  look  quite  different,  but  it 
would  still  reflect  the  same  measurement.  Hence  unit  systems  and  number  systems  are  coupled 
together  and  are  used  to  express  quantity. 

L2.  Coordinate  Systems 

The  coordinate  system  in  the  above  example  is  a  conventional  two  dimensional  cartesian  system 
based  on  a  locally  flat  plane  tangent  to  the  earth.  It  is  used  to  represent  the  direction  of  the  velocity. 
Other  systems  include  the  spherical  polar  coordinate  system,  cylindrical  polar  coordinates,  and 
curvilinear  coordinates,  and  there  are  many  others.  If  another  coordinate  system  were  used,  it 
would  have  changed  the  word  ‘‘northeasterly”  to  something  else,  but  it  would  not  change  the 
measurement’s  direction  -  it  changes  only  how  the  direction  is  expressed.  A  coordinate  system  is 
therefore  an  artifice  used  to  express  direction  of  a  measurement  For  a  scalar,  the  coordinate  system 
is  just  the  convention  for  assigning  a  plus  or  minus  arithmetic  sign. 

L£  Reference  Frames 

The  idea  of  a  reference  firame  is  deeply  imbedded  in  mechanics.  Reference  frames  are  real  objects 
as  expressed  by  Newton  himself.  Reference  frames  are  distinguished  from  each  other  by  their  state 
of  motion.  Measurements  always  express  a  quantity  with  respect  to  two  reference  frames,  since 
even  such  scalars  as  mass  vary  with  reference  frame  in  general.  In  conversation,  the  phrase 
<quantity>  of  <teference  frame  1>  relative  to  <reference  frame  2>  is  often  used  to  express  these 
two  fiames  of  reference.  Using  velocity  as  an  example,  the  frame  in  which  a  velocity  is  measured 
expresses  the  agreed  upon  notion  of  zero  velocity,  for  that  particular  measurement  The  frame 
possessing  the  velocity  is  the  thing  which  is  considered  to  be  moving,  for  that  particular 
measurement  Only  dungs  have  a  velocity.  The  velocity  of  emp^  space  is  meaningless. 

Reference  Frames  vs.  Coordinate  Systems 

One  often  overlooked  fact  which  one  must  face  squarely  when  doing  advanced  dynamics  is  that  a 
reference  frame  and  a  coordinate  system  ate  different  The  former  is  a  concept  of  physics,  the 
second  a  concept  of  mathematics,  and  they  are  absolutely  not  the  same  thing.  There  are  at  least  four 
common  practices  today  which  obscure  this  distinction.  The  first  is  the  use  of  the  same  term 
“vector^’  in  both  linear  algebra,  and  in  physics,  when  two  different  concepts  are  involved.  Ihe 
second  is  the  lack  of  emphasis  of  the  distinction  between  the  tensor  of  second  order  and  the  matrix. 
The  third  is  the  hybrid  term  “coordinate  frame”  as  it  is  used  in  computer  graphics  and  robotics.  The 
fourth  is  the  wid^pread  and  sometimes  confusing  practice  of  the  use  of  random  expressions  lil^ 
“coordinate  reference  system”,  “reference  coordinate  system”,  “coordinate  reference  frame”,  etc., 
in  many  texts.  This  may  be  justified,  however,  because  the  issue  can  do  mote  to  confuse  the  student 
than  to  illuminate  the  basics.  In  the  work  here,  the  issue  is  of  central  concern. 
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lA  Pcinciple  of  Coordinate  System  Invariance 

In  linear  algebra,  when  applied  to  physics,  its  vectors  are  often  the  mathematical  representation  of 
physical  vectors  in  a  particular  cartesian  coordinate  system.  In  physics,  the  vector  notation  is 
designed  to  emphasize  that  the  laws  of  physics  hold  irrespective  of  the  coordinate  system  in  which 
they  are  represented.  This  is  called  the  principle  of  coordinate  system  invariance‘s.  Newton’s 
second  law  for  instance  says  that  force  and  acceleration  are  parallel  vectors  in  the  physics  sense, 
whose  amplitude  ratio  is  the  mass,  irrespective  of  coordinate  system,  provided  they  are  both 
measured  with  respect  to  an  inertial  reference  frame  and  provided  they  are  the  net  force  on  and 
actual  acceleration  of  the  same  body.  Hence  the  choice  of  different  coordinate  systems  gives 
different  vectors  in  the  linear  algebra  sense,  but  the  same  vector  in  the  physical  sense. 

LL  Transformation  Laws  and  Conversion  Formulae 

In  parallel  to  the  coordinate  system,  reference  frame  distinction,  there  is  a  related  distinction 
between  transforming  vectors,  and  transforming  their  representations  in  a  particular  coordinate 
system.  The  first  involves  physics,  the  second  mathematics. 

Consider  an  experiment  where  a  person  stands  on  a  moving  train  and  throws  a  baseball  toward  the 
front  of  the  train  as  shown  in  Figure  1: 


Let  the  person  in  the  train  be  called  the  moving  observer  and  the  one  on  the  side  of  the  track  be 
called  the  stationary  observer.  For  the  stationary  observer  measuring  the  ball’s  velocity,  it  is  the 
laws  of  physics  which  permit  a  prediction  of  the  measurement  made  by  that  observer  based  on  the 
speed  of  the  train  and  die  speed  of  the  ball.  In  Newtonian  mechanics  it  is  known  that: _ 


'bs 


^bm  ^ms 


Now  this  is  a  law  of  physics  which  relates  measurements  made  by  observers  in  reference  frames. 
It  turns  out  that  it  involves  the  mathematical  operation  of  addition,  as  defined  on  physical  vectors. 
However,  as  velocities  approach  the  speed  of  light,  a  different  law  of  physics  applies,  namely  the 
relativistic  velocity  addition  equation: 


4.  This  is  why  you  can  aibitrarily  “assume”  a  direction  of  a  scalar  quantity,  and  it  will  come  out  negative  if 
ymi  “assumed”  wrong .  You  actually  are  assuming  nodting,  but  rather  are  establishing  a  scalar  coordinate  sys¬ 
tem. 
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allows  prediction  of  a  measurement  of  force  taken  in  an  inertial  frame,  from  a  measurement  of 
mass  and  of  acceleration  taken  in  the  same  inertial  frame.  For  physical  vectors,  change  of 
coordinate  system  does  not  change  the  reference  frames  involved  or  the  quantity  measured  in  any 
way. 

Now  the  coordinate  system  conversion  formulae  of  algebra  solve  an  entirely  different  problem. 
They  permit  conversion  between  two  coordinate  systems  for  a  single  measurement  taken  in  a  single 
frame.  In  summary: 

•  The  laws  of  physics  permit  prediction  of  measurements  based  on  other  measurements 
provided  all  measurements  are  taken  in  the  appropriate  reference  frames  for  the  law  to  tq)ply. 
Different  laws  impose  different  limitations  on  reference  frames. 

•  The  mathematical  manipulation  of  equations  involving  physical  quantities  require  that  all 
quantities  first  be  expressed  in  a  common  coordinate  system,  but  that  choice  of  coordinate 
system  is  arbitrary. 

This  material  is  presented  here  because  a  naive  implementation  of  inertial  navigation  leads  to  two 
basic  problems.  These  problems  come  down  to  handling  coordinate  systems  incorrectly,  and 
handling  reference  frames  incorrectly.  The  solution  to  the  problem  is  twofold.  First  determine  die 
laws  of  physics  which  allow  a  conversion  from  whatever  an  accelerometer  measures  to  the 
measurement  required.  Second,  express  all  such  measurements  in  a  common  coordinate  system. 
Then,  a  straightforward  integration  provides  the  vehicle  position. 

In  order  to  discover  what  an  accelerometer  measures,  it  is  necessary  to  find  a  particular  fimne  of 
reference  called  an  inertial  frame  of  reference,  and  understand  how  it  is  affected  by  gravitational 
fields.  In  order  to  do  that,  it  is  necessary  to  confront  the  basic  postulate  of  the  theory  of  general 
relativity.  \^th  these  tools,  it  will  be  possible  to  address  the  first  half  of  the  problem.  The  first  step 
in  that  process  is  a  brief  discussion  of  the  effects  of  motion  of  observers  on  the  measurements  of 
physical  quantities. 


Notice  that  mis  law  of  pnysics  permits  combmation  ot  measurements  taken  m  two  dittet^t 
reference  frames  in  order  to  predict  a  third  measurement  that  was  never  made.  Newton’s  second 
law: 


F  =  ma 
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2  Frames  in  Relative  Motion 

Accelerometers  intrinsically  measure  with  respect  to  an  inertial  frame  of  reference  because  they 
are  based  on  Newton’s  Laws.  This  section  will  investigate  the  notion  of  an  inertial  frame  and  its 
implications  to  the  navigation  problem. 

2J.  Inertial  Eramss 

Experience  with  physical  measurements  points  to  the  existence  of  certain  distinguished  frames  of 
reference  -  the  inertial  flrames  -  which  are  required  for  the  expression  of  certain  laws  of  physic 
Newton  believed  in  absolute  space  and  time: 

“Absolute  space,  in  its  own  nature,  without  relation  to  anything  external,  always  remains 
similar  and  immovable” 

but  general  relativity  has  changed  ideas  of  space  and  time  since  then.  It  is  customary  now  to  speak 
of  inertial  frames  rather  than  absolute  frames.  There  are  an  infinity  of  inertial  frames  and  they 
pervade  all  of  space  just  as  did  Newton’s  absolute  ones. 

If  the  existence  of  inertial  frames  and  Newton’s  laws  are  taken  as  axiomatic,  then  inertial  frames 
can  be  defined  as  those  frames  in  which  Newton’s  laws  apply.  Under  this  operational  definition,  the 
earth  is  not  an  inertial  frame  because  an  object  dropped  by  an  observer  falls  to  the  ground  along  a 
nonstraight  trajectory.  However,  a  spacecraft  in  orbit  is  such  a  frame,  as  is  a  free-faUing  aircraft,  as 
is  a  point  in  space  far  removed  from  any  massive  object 

The  simple  existence  of  a  gravitational  field  apparently  does  not  disqualify  a  frame,  for  while  the 
observer  of  a  dropped  object  is  not  in  an  inertial  frame  on  earth,  the  free-falling  ball  itself  is  an 
inertial  frame^.  If  one  inertial  frame  exists,  then  all  other  frames  translating  uniformly  with  respect 
to  it  are  also  inertial^.  Any  frame  rotating  with  respect  to  it,  carmot  be  inertial  since  rotation  about 
any  point  but  the  center  of  mass  requires  acceleration. 

Another  definition  of  an  inertial  frame  is  one  where  the  sum  of  all  accelerations  and  gravitational 
attractions  are  zero.  In  the  language  of  general  relativity,  the  acceleration  is  “transformed  away”. 
By  and  large,  in  an  inertial  frame,  the  gravitational  attraction  of  the  sun  and  other  heavenly  bodies 
can  be  ignored.  These  attractions  can  be  ignored  because  the  motions  of  any  bodies  close  to  the 
earth  are  affected  by  the  sun  just  as  the  earth  is,  so  the  effect  cancels  out.  However,  as  bodies  move 
further  from  the  earth,  this  is  no  longer  so.  For  this  reason,  general  relativity  admits  “local  inertial 
frames”,  where  die  laws  of  physics  apply  in  a  local  region.  Such  frames  are  central  to  the 
description  of  both  mechanical  and  electromagnetic  phoiomena. 


5.  Air  fiictioo  is  ^noied  tbroashontthe  docamenL 

6.  Bat  there  ate  still  others,  see  die  next  section 
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A  fundamental  axiom  of  physics  is  the  principle  of  relativity.  This  principle  originally  held  that 
the  laws  of  mechanics  are  the  same  in  all  of  Newtons  “absolute”  frames.  Newton  meant  it  to  apply 
only  to  mechanics: 


*The  motions  of  bodies  included  in  a  given  space  are  the  same  among  themselves  whether 
that  space  is  at  rest  or  moves  uniformly  forward  in  a  straight  line.” 

Another  often  used  way  of  saying  this  is  that  all  absolute  frames  ate  equivalent  or  that  the  laws 
themselves  are  invariant  Einstein  called  this  the  restricted  sense  of  the  principle.  In  his  words: 

“Even  though  classical  mechanics  does  not  supply  us  with  a  sufficiently  broad  basis  for  the 
theoretical  presentation  of  all  physical  phenomena,  still  we  must  grant  it  a  considerable 
measure  of  truth  since  it  supplies  us  with  the  actual  motions  of  the  heavenly  bodies  with  a 
delicacy  of  detail  little  short  of  wonderful  The  principle  of  relativity  (in  the  restricted 
sense)  must  therefore  apply  with  great  accuracy  in  the  domain  of  mechanics.” 

He  was  getting  at  the  fact  that  it  ^patently  does  not  apply  to  electromagnetic  phenomena. 
Carefully  examination  of  the  word  apparently  ted  to  Special  Relativity.  Careful  examination  of  the 
requirement  for  absolute  frames  led  to  General  Relativity.  In  his  General  Relativity  theory,  Einstein 
asserts  that: 


“we  wish  to  understand  by  the  general  principte  of  relativity  the  following  statement:  all 
bodies  of  reference  are  equivalent  for  the  description  of  natural  phenomena,  whatever  may 
be  their  state  of  motion” 


Here,  he  is  getting  at  the  fact  that  the  principle  holds  in  inertial,  not  absolute  frames.  For  the  present 
purpose,  the  point  is  that  inertial  and  gravitational  mass  are  equivalent,  and  therefore  frames  in 
gravitational  free-fall  are  inertial. 

Einstein  argued  this  with  a  thought  experiment  consisting  of  a  man  in  an  elevator  out  in  space  and 
some  other  agent  pulling  with  constant  acceleration  “up”  on  the  entire  elevator  with  the  man  in  it 
as  shown  in  Figure  2. 


Figure  2  Einstein’s  Elevator  Thought  Ejqpeiiment 


I^m  the  point  of  view  of  the  man  in  the  elevator,  he  seems  to  be  experiencing  a  gravitational 
attraction  pulling  him  toward  the  bottom  of  the  elevator.  He  feels  that  he  has  some  wdght  The 
principle  says  that  no  acperiment,  conducted  in  Uu  frame  of  r^erenee  of  the  man  in  Ae 
etemOof,  canted  the  Terence. 
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The  fact  that  astronauts  training  in  a  free-falling  aircraft  experience  weightlessness  is  not  an  aspect 
of  human  subjectivity.  Even  with  the  best  instruments  they  carmot  tell  the  difference,  and  since  in 
physics,  things  that  carmot  be  measured  do  not  exist,  they  are  indeed,  in  a  very  deep  sense,  actually 
weightless.  With  this  final  idea,  it  is  possible  to  find  the  inertial  frame  necessary  for  expressing  the 
equations  of  inertial  navigation. 

2.1  Hitgrem  ftf  Cffri<>lls 

It  is  clear  now  that  accelerometers  on  an  INS  that  is  stationary  with  respect  to  the  surface  of  the 
earth  will  register  an  acceleration  due  to  the  rotation  of  the  system  with  the  earth.  In  order  to 
investigate  this  mote  fully,  the  amount  of  this  acceleration  will  be  quantified. 

Consitter  two  frames  of  reference  that  are  rotating  with  respect  to  each  other  at  some  instantaneous 
angular  velocity  OH,  as  measured  in  the  first  frame,  and  that  have  an  instantaneous  relative  position 
of  p  as  measured  in  the  first  frame.  The  first  one  will  be  caSltdfixed  and  the  second  moving  -  though 
this  choice  is  completely  arbitrary.  Quantities  measured  in  the  first  frame  will  be  indicated  by  an  f 
subscript  and  those  measured  in  the  second  will  be  indicated  by  an  m  subscript  Suppose  an 
observer  measures  the  position  (r),  velocity  (v),  and  acceleration  (a)  of  an  object  in  the  moving 
frame,  and  it  is  necessary  to  know  what  a  second  observer  would  measure  in  the  fixed  frame  for 
the  motion  of  the  same  object  The  situation  is  depicted  in  Figure  3. 


Cteariy,  the  positions  add  by  simple  vector  addition: 


In  order  to  compute  the  time  derivative  a  very  important  general  formula  can  be  used  to  relate 
moving  reference  frames,  as  t^Ued  to  any  vector  v.  It  is  a  tedious  but  straightforward  exercise  to 
show  that: 


The  derivation  proceeds  by: 

•  assigning  coordinate  systems  to  the  fixed  and  moving  frames 


7.  That  is,  be  cannot  Icxdc  om  a  window. 
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•  writing  the  vector  v  in  terms  of  the  unit  vectors  of  both  frames 

•  using  the  chain  rule  of  differentiation  noting  that  the  unit  vectors  of  the  moving  frame  are 
time  varying 

This  formula  will  be  called  the  theorem  of  Coriolis^.  It  is  imperative  to  keep  the  reference  frames 
straight  when  using  this  equation.  The  equation  is  used  to  relate  time  derivatives  of  the  same  vector 
computed  by  two  observers  in  relative  rotational  motion.  The  vector  is  the  seme  in  all  its  instances 
in  tlie  equation.  The  derivatives,  however,  are  different  The  vector  may  be  the  position,  velocity, 
acceleration,  force,  field  strength,  or  anything  else  of  anything  with  respect  to  anything  else.  If 
there  are  physical  transformation  laws  associated  with  the  translational  motion  of  the  two  frames, 
these  are  considered  outside  this  equation. 

The  intuitive  meaning  is  that  die  derivative  computed  by  the  moving  observer  is  the  derivative  seen 
by  the  stationary  observer  plus  the  component  caused  by  his  own  rotational  motion.  The  angular 
velocity  is  the  angular  velocity  of  the  moving  pome  respect  to  the  fixed  one,  and  the  last 

instance  of  the  vector  in  the  equation  is  the  original  vector  again. 

The  first  step  is  to  differentiate  the  position  relationship  in  the  fixed  frame.  This  step  is  a 
straightforward  vector  time  derivative  which  has  nothing  to  do  with  the  above  derivative 


The  fixed  observer  sees  an  extra  two  components  in  addition  to  the  velociQr  seen  by  the  moving 
observe.  The  second  term  is  due  to  the  relive  linear  velocity  of  the  moving  frame.  The  third  term 
is  the  motion  of  the  vector  due  to  the  angular  velocity  of  the  moving  frame.  If  the  fixed  observer 
observes  his  position  vector  to  be  fixed,  the  moving  observer  wiU  observe  circular  motion  due  to 
the  motion  of  the  moving  frame.  Application  of  the  theorem  to  die  velocity  relation  generates: 

There  are  four  extra  components  of  acceleration  in  the  equation.  The  Einstdn  acceleration  is  due 
to  d»  linear  acceleration  of  the  moving  frame.  It  is  the  same  as  the  sqiparent  gravity  sensed  by  the 
observo^  in  die  elevator.  The  Eulo*  acceleration  is  due  to  the  angular  acceleration  of  the  moving 
frame. 

The  Coriolis  acceleration  is  due  to  the  angular  velocity  of  the  moving  frame  and  the  velocity  of  the 
objea  as  seen  in  the  moving  frame.  Objects  do  not  fall  in  straight  lines  on  the  earth  because  of  this 


8.  AUiongb  tbe  tenn  apfriies  strictly  only  to  the  accelefation  foonula,  it  will  be  used  to  desoibe  any  tiansfor- 
mation  of  a  vector’s  derivative. 
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acceleration  seen  in 
fixed  frame 
acceleration  seen  in 
moving  firame 

Einstein 
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tuxeleraiion.  me  uentnpetai  acceleration  is  also  called  Centnlugai.  it  is  due  to  me  angular  velocify 
of  the  moving  frame  and  the  offset  position  of  the  object  in  the  moving  frame. 


In  order  to  see  which  observer  sees  the  extra  accelerations,  identify  the  fixed  frame  with  an  inertial 
frame,  and  the  moving  firame  with  a  spinning  firame.  The  inertial  observer  is  the  only  one  who 
observes  motion  according  to  Newton’s  second  law,  so  his  measure  of  the  acceleration  is  the  one 
that  must  be  used  in  that  law.  Assume  that  there  are  no  forces  on  the  object,  so  that  its  inertial 
acceleration  is  zero,  and  that  its  velocity  is  zero.  If  this  is  so,  tl^n  the  fixed  observer  sees  a 
stationary  object  and  the  moving  observer  observes  all  of  these  extra  terms  of  acceleration,  even 
though  there  are  no  forces  on  the  object  He  would  naturally  assume  diat  there  are  forces  causing 
these  accelerations.  These  forces  are  called  i^tparent  frarces.  Sometimes  they  are  referred  to  as 
inortial  f<Nroes  ■  which  can  lead  to  confusion,  fainause  they  are  not  measured  in  an  inertial  frame. 
They  are  caused  by  the  inertial  properties  of  matter  as  seen  by  an  observer  in  a  non-inertial  frame, 
and  not  by  the  interaction  of  otter  bodies. 

Li  an  inertial  firame,  the  only  forces  on  a  body  required  in  Newton’s  laws  are  titose  arising  from 
the  interaction  of  other  material  bodies,  namely  those  caused  by  contact  (which  is  essentially 
electrostatic),  electricity,  magnetism,  or  gravit^on.  If  this  is  so,  it  is  natural  to  a^  why  it  is  that 
gravitation  seems  to  have  a  distinguished  rote  among  the  fundamental  forces  of  nature.  For 
Einstein  did  not  say  that  firames  in  free-fall  under  the  influence  of  an  electric  or  magnetic  field  are 
initial. 


The  simple  answer  is  diat  Newton’s  laws  say  so,  but  this  does  not  address  the  why  of  the  questioiL 
According  to  Newton’s  laws,  even  though  the  forces  on  a  body  may  result  firom  the  charged 
particles  that  it  contains,  another  body  in  contact  with  it,  the  current  flowing  through  it,  or  the 
gravitational  field  that  attracts  it,  it  is  the  mass  of  the  body  that  determines  its  motion.  There  are 
laws  of  magnetic,  electric,  and  gravitational  intoaction  for  determining  forces,  but  there  is  only 
one  law  for  determining  accelerations,  and  that  law  uses  tte  mass,  not  the  electric  charge  or  current 

Consider  an  electron  in  a  television  set  that  is  bdng  bent  firom  its  natural  trajectory  in  order  to  have 
it  hit  the  cmrrect  place  on  tte  screot  It  is  die  charge  of  tte  etectron  that  determines  how  roach  farce 
can  be  exerted  on  it  by  an  applied  voltage.  It  is  the  mass  of  tte  electron  which  determines  how 
much  acceleration  it  experience  to  bend  it  before  it  hits  tte  screen. 
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3  Rotational  Dynamics 


This  section  will  investigate  the  properties  of  the  spinning  rotor,  a  body  which  exhibits  axial 
symmetry  and  spins  about  that  axis,  liiese  ideas  form  the  basis  of  the  operation  of  the  gyroscope. 

2LL  Fwie  Rotatigg  rf  An  Arbitrary  Body 


Consider  Newton’s  second  law  as  ^plied  to  a  purely  rotating  mass,  and  expressed  in  an  inertial 
frame  (i  subscript)  in  the  angular  momentum  form.  Often  called  Euler’s  equation,  it  applies  only 


ddiS) 

“  dt 

Net  Torque  N 

Angular  Momentum 

L 

Inertia  Tensor  I 

Angular  Velocity 

Undo:  the  assumption  that  a  body  experiences  no  net  torques,  one  would  expect,  in  analogy  to 
translational  dynamics,  that  the  motion  is  one  of  constant  angular  velocity.  This  is  not  the  case.  The 
equation  becomes  on  integration: 

=  IjC5  =  const 


Again  this  only  holds  in  an  inertial  frame,  which  is  a  problem,  because,  for  a  body  of  arbitrary 
sh^,  the  inertia  tensor  of  the  rotating  body  will  not  be  constant,  even  under  a  constant  angular 
velocity.  If  the  inertia  tensor  is  not  constant,  the  angular  velocity  vector  must  also  change  in  order 
for  the  product,  the  angular  momentum,  to  remain  constant  The  conclusion  is  that  the  general 
torque  ftee  motion  of  a  rigid  body  is  a  tumbling  motion.  This  is  why  arbitrarily  shaped  objects 
tumble  -  even  in  space  -  instead  of  ginning  neatly. 


In  the  very  special  case  of  a  line  qmimetric  body,  rotating  about  its  axis  of  symmetry,  the  inertia 
tensor  is  constant  in  inertial  space  .  Under  these  conditions,  it  can  be  conclud^  that 


(0  =  const 


It  has  been  shown  that  a  body  under  no  net  jq>plied  torque,  which  is  rotating  around  its  axis  of 
^mmetty,  will  remain  rotating  about  its  line  of  symmetry.  This  property  is  called  rigidity  of  the 
q;nn  axis.  An  attendant  issue  is  stability.  It  can  be  shown  t^t,  provid^  the  axis  of  symmetry  is  the 
axis  of  smallest  inertia,  the  motion  is  stable  to  small  perturbations. 


9.  Assomiiig,  of  coarse,  imifonn  densi 
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^  Forced  Rotation  of  an  Arbitrary  Body 

Another  way  to  achieve  pure  spinning  is  to  constrain  all  other  motion  physically.  A  body  whose 
angular  acceleration  is  constrained  in  this  way  must  still  follow  Newton’s  laws.  Consider  such  a 
body  mounted  in  such  a  way  that  it  is  free  to  spin  on  one  axis  only.  It  is  clear  that  such  a  body  will 
have  a  constant  direction  of  angular  velocity.  However,  if  the  spin  axis  does  not  coincide  with  an 
axis  of  symmetry,  then  the  inertia  tensor  will  not  be  constant  Hence  the  net  torque  on  the  body  will 
not  be  constant  These  changing  torques  mu^  be  supplied  by  the  constraint  tarings  which  will 
therefore  undergo  frequent  oscillatory  loading^^. 

^  Free  Precession  of  a  Symmetric  Body 

The  next  to  previous  section  had  assumed  that  axial  spin  was  the  case  and  shown  that  it  was 
consistent  with  the  equation  of  motion.  It  turns  out  that  there  are  other  motions  possible  for  a 
annning  rotor.  To  get  these,  consider  again  Euler’s  equation: 


Identifying  the  fixed  frame  with  the  inertial  frame  (subscript  i),  and  the  moving  frame  with  a  frame 
fixed  to  the  body  (subscript  b),  the  theorem  of  Coriolis  can  be  applied.  This  gives  after  some 
manipulation: 


= 


,  +  ©  X  Leo 

dt  Jb 


This  result  is  an  expression  of  Euler’s  equation  in  a  very  subtle  form.  It  is  important  to  be  very 
careful  about  what  ^  above  symbols  mean.  The  first  term  on  the  right  hand  side  is  the  body  fixed 
observer’s  derivative  of  die  inertial  observer’s  angular  momentum^^  The  first  ©  in  the  cross 
product  is  the  relative  rotation  of  the  fixed  and  rotating  frames.  The  second  ©  is  the  rotation  of  the 
body  in  an  inertial  frame.  Only  when  the  fixed  frame  is  inertial  and  the  second  is  body  fixed  do 
these  two  angular  velocities  become  equal. 


10.  This  is  die  wheds  of  cars  ate  balanced,  and  why  it  is  done  by  adding  small  weights  in  strategic  places. 

1 1 .  Yes,  diat  sounds  like  a  very  strange  dung  to  be  evainadng.  But  this  is  exaedy  what  the  theorem  of  Quiolis 
is  for. 


Modem  hienial  and  Satdlite  Navigaden  Systems 


Now,  consider  expressing  these  eouations  in  the  coordinate  system  in  which  the  inertia  matrix  if 
diagonal,  called  the  principal  axes^^: 


It  is  easy  to  vaify  that  the  two  equations  togetiier  describe  a  circle,  and  the  last  a  constant  The 
12.  By  a  ttwocem  of  hiiear  aigdna,  all  bodies  of  any  shape  possess  princ^  axes. 
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conclusion  is  that  the  unforced  symmetric  body  has  an  angular  velocity  which  describes  a  cone 
in  space.  Now  the  obvious  question  is  which  observer  sees  it  this  way.  The  angular  velocity  in  the 
formulation  is  the  angular  velocity  of  the  fixed  body  frame  with  respect  to  inertial  space.  This 
quantity  was  expressed  in  the  principle  body  coordinate  system. 

It  is  tempting  to  say  that  the  body  fixed  observer  sees  the  conical  motion.  However,  remember  that 
the  angular  acceleration  vector  is  inertially  referenced  and  only  represented  in  body  coordinates. 
In  an  inertial  frame,  for  torque  free  motion,  the  angular  momentum  is  constant,  and  the  kinetic 
energy  is  constant  The  latter  is  given  by: 


Thus,  the  angular  momentum  vector  and  the  angular  velocity  vector  make  a  fixed  angle,  and  the 
angular  momentum  itself  is  constant  Hence,  the  velocity  vector  rotates  around  the  fixed 
momentum  vector  while  maintaining  a  constant  angle  with  it  To  the  inertial  observer,  the  object 
spins  about  an  axis  which  is,  in  general,  not  the  symmetry  axis,  and  the  axis  itself  spins  in  a  cone. 
This  conical  motion  of  the  spin  axis  is  called  free  precession^^ 

fQKsd  Prectsaon 

Consider  now  what  happens  when  a  rotating  symmetric  body  is  subjected  to  a  torque.  It  is  clear 
from  Euler’s  equation  that  when  the  2q>plied  torque  is  parallel  to  the  angular  momentum,  the  body 
experiences  an  angular  acceleration  which  changes  the  magnitude  of  the  angular  velocity  but  not 
its  direction.  However,  that  is  where  die  analogy  to  the  translational  case  ends.  Strange  as  it  may 
seem,  the  equation  implies  that  the  response  of  the  rotor  to  a  torque  applied  normal  to  the  spin  axis 
is  to  rotate  in  a  direction  normal  to  both  the  spin  axis  and  the  applied  torque.  This  phenomenon  is 
also  called  precession 

That  is,  if  one  attempts  to  rotate  the  spin  axis,  it  does  not  rotate  in  the  direction  that  the  torque  is 
trying  to  force  it  to  go.  Consider  Figure  4.  The  vector  algebra  makes  it  seem  obvious,  but  the  r^ty 
seems  counterintuitive,  because  the  body  wiU  tilt  into  the  torque.  It  can  be  understood  intuitively 
by  considering  the  incremental  velocities  applied  to  particles  on  the  equator  of  the  body.  The  torque 
imparts  an  incremental  upward  velocity  to  the  particle  shown  in  the  right  side  figure,  and  an 
incremental  downward  velocity  on  its  twin  behind  the  body. 

When  added  to  its  current  velocity,  the  net  result  is  a  rotation  of  the  velocity  vector  about  an  axis 


13.  This  phenonaenoa  is  lelaiBd  to  ibe  lotatioo  tbe  eardi’s  siTO  axis,  which  is  also  called  precessicHi.  Qearly 
the  earth  itself  is  a  ^phming  rotor. 

14.  Ancient  astronomm  noticed  that  the  oonstdlation  in  which  the  son  rises  (» the  vernal  equinox  advances 
or  **ixeoedes”.  Hence  the  term. 

15.  There  is  much  confusion  about  die  terms  precession  and  nutation  in  ii»  literature.  Some  audKvs  call  the 
finee  precessum  described  above  “nutation”,  reserving  precession  for  the  gyrosct^  reqionse  to  torques.  Oth¬ 
ers  define  precesstm  as  drove  and  reserve  rmtadou  to  small  variadons  of  the  shape  (rf  the  cone  swqx  out  by 
die  qiin  axis.  This  is  die  sense  of  nutation  as  used  with  reflect  to  the  earth’s  motion.  Others  reserve  nutation 
for  die  noddiig  of  die  qsBi  axis  of  a  tcqi.  Still  others  define  three  rates  of  Eukr  smgles  called  {secession,  nuta- 
don,  and  qiin.  Ihe  auditor  has  been  miable  to  find  a  strict  equivalence  in  these  uses  of  the  terms. 
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Figure  4  Precession 


If  the  angle  between  the  new  and  old  angular  acceleration  vectors  is  defined  as  d6,  and  its  rate  is 
defined  as  Q,  then  the  gyroscope  rigidity  equation  is: _ 

^  IdCI  ^  WIdt 

iLi  |Ll 

=  H 

|l| 

Remember  that  Q  is  the  rate  of  change  of  the  angular  velocity  vector,  not  the  rotation  of  the  body. 
Considering  the  orthogonality  of  the  precession,  this  equation  is  often  expressed  as: _ 

A  A  A 

N  =  QxL 

where  only  the  direction  of  the  vectors  is  considered.  So  the  rate  of  change  of  the  angular  velocity 
is  proportional  to  the  applied  torque,  and  inversely  proportional  to  the  angular  velocity.  This  is  the 
justification  for  the  experimental  observation  that  the  faster  a  rotor  spins,  the  harder  it  is  to  disturb 
it  This  technique  is  called  spin  stabilization  and  it  is  commonly  used  for  communication 
satellites. 
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4  Motion  on  the  Earth 


The  earth  is  not  an  inertial  frame  because  objects  do  not  fall  to  it  under  the  influence  of  its 
gravitation  according  to  Newton’s  laws.  This  section  will  investigate  why  this  is  so  in  more  detail. 

4L  Eailii  as  an  In?.ilMErams 

The  motion  of  the  earth  through  space  is  a  complicated  superposition  of  several  component 
motions.  The  most  important  of  these  are: 

•  daily  rotations  about  its  own  axis 

•  aimual  revolution  about  the  sun^^ 

•  precession  and  nutation  of  the  axis^^ 

•  motion  of  the  sun  in  the  galaxy 

•  motion  of  the  galaxy  in  the  cosmos 

This  section  will  investigate  the  reasons  why  the  first  of  these  is  the  only  one  that  is  important  to 
the  problem  of  finding  an  inertial  frame. 

4.1.1  Daily  Rotation 

This  is,  of  course,  a  rotation  of  360  degrees  per  day.  However,  there  is  a  sublety  of  importance.  The 
solar  day  must  be  distinguished  from  the  sidereal  day.  The  first  is  a  complete  rotation  with  respect 
to  the  sun  and  is  24  hours^^  long.  The  second  is  a  complete  rotation  with  respect  to  the  stars  and  is 
about  4  minutes  shorter.This  is  a  consequence  of  the  fact  that,  if  the  earth  did  not  rotate  on  its  axis, 
the  sun  would  still  rise  and  fall  once  per  year  as  the  earth  orbits  it,  so  the  sun  goes  around  the  earth 
slightly  faster  than  the  stars  do  as  ^own  in  Figure  5.  If  the  stars  are  an  inertial  frame,  it  is  the 
sideted  rotation  of  the  earth  that  matters.  Also,  since  the  earth-fixed  frame  rotates  with  respect  to 
the  stars  at  a  considerable  rate,  it  cannot  be  inertial.  The  extent  of  this  effect  will  be  calculated 
shortly. 


16.  As  yoa  read  this  you  are  travelling  tiuougb  q)ace  at  70,000  in.pJi.  with  respect  to  the  sun.  Hold  on. 

17.  The  ccHistellatitnis  trf  die  zodiac  are  those  constellations  which  appear  on  the  horizon  at  sunrise.  Ancient 
Egyptians  noticed  a  correlatioa  between  the  constellations  and  the  nile  flood  and  assumed  thm  was  a  OHisal 
relatioiish^  -  hence  astrcflogy  was  bom.  Due  to  the  aimual  circuit  about  the  sun,  these  consmllatitxis  ai^iear 
to  cirde  the  hofiztm  mce  a  year.  Due  to  precession,  the  zodiac  also  rotates  once  every  26,000  years,  so, 
according  to  standard  tables,  everybody’s  sign  is  now  off  by  one  constellation  and  things  won’t  be  right  again 
until  the  year  24000.  Astrologers,  naturally,  don’t  emidiasize  this. 

18.  The  hour  is  defined  so  that  this  is  so. 
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Another  sublety  is  the  change  in  the  earth’s  orbital  speed  about  the  sun  due  to  its  elliptical  shape. 
Figure  5  demonstrates  that  the  solar  day  varies  in  length  for  this  reason.  The  familiar  24  hour  day 
is  the  mean  solar  diay.  Of  course,  the  amount  of  daylight  per  day  varies  with  the  season  as  a 
consequence  of  the  tilt  of  the  earth’s  axis  at  roughly  22  degrees^^  to  the  ecliptic^,  but  the  current 
issue  is  the  day  when  it  is  defined  as  a  complete  rotation  of  the  earth. 

4.1,2  Annual  Revolution  About  die  Sun 

While  the  earth’s  velocity  varies  with  the  season  of  the  year,  ignoring  daily  rotation,  this 
acceleration  is  a  translational  fi:ee-fall  because  it  is  completely  caused  by  the  gravitational 
attraction  of  the  sun  on  the  mass  center  of  the  earth.  The  earth  is  also  in  gravitational  firee-fall  with 
respect  to  all  other  bodies  in  the  universe,  so  a  geocentric  earth  frame  which  does  not  rotate  about 
the  earth’s  axis  must  be  inertial  by  the  principle  of  relativity. 

4.13  Precession 

However,  there  is  one  sublety.  The  earth  is  not  perfectly  roimd,  so  the  mass  attraction  of  the  moon, 
sun,  and  planets  on  the  bulge  of  the  equator  cause  precession^^  of  the  axis.  These  attractions  try  to 
pull  the  earth’s  axis  perpendicular  to  the  ecliptic,  and  they  actually  cause  precession  because  of  the 
earth’s  rotational  motion.  Precession  amounts  to  small  changes  in  the  angular  velocity  vector  of  the 
earth  about  its  own  axis,  and  as  was  shown  earlier,  gives  rise  to  an  effect  called  Euler  apparent 
force.  Lucidly,  the  precessional  rate  is  360  degrees  every  26000  years^  and  the  effect  is  so  small 
as  to  be  undetectable.  The  discussion  could  stop  here  since  an  inertial  frame  has  been  found,  but  it 
will  continue  for  reasons  that  will  be  clear  later. 


19.  Ibis  is  why  the  trtqncs  of  cancer  and  capticrsn  ate  22  degrees  tioin  the  equator. 

20.  The  plane  of  the  earth’s  orbit  about  the  sun. 

21.  The  so-called  precession  of  the  equinoxes  was  sqjparently  discovered  by  die  Babylonian  astronomerQde- 
nas  in  about  343  BC. 

22.  The  ice  ages  ate  tied  to  this  cycle.  They  are  as  predictable  as  the  seasons  of  the  year.  Buy  boots  btf ore  the 
winter  of  24000  AD. 
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4.1.4  Motion  in  the  Galaxy 

The  sun  orbits  the  center  of  the  galaxy  at  a  period  of  one  revolution  every  200  million  years  with 
a  tangential  velocity  of  roughly  30,000  m.p.h.^^.  This  is  gravitational  free-fall. 

4.1.5  Motion  in  the  Cosmos 

According  to  astronomical  observations,  all  galaxies  are  flying  apart  at  incredible  speeds.  Although 
the  speeds  are  so  remarkable  that  relativistic  ted  shifts^^  of  starlight  are  the  rule,  the  accelerations 
are  completely  determined  by  gravitation.  This  is  gravitational  free-fall. 

4. 1.6  Motion  of  the  Cosmos 

To  be  complete,  it  is  necessary  to  ask  the  question  of  whether  the  center  of  mass  of  all  the  universe 
can  be  in  motion.  If  it  can,  with  respect  to  what  would  it  be  in  motion?  There  is  nothing  else  by 
definition.  Physicists  have  pondered  this  question  deeply,  and  some  believe  that  the  paradox  is  tied 
to  the  existence  of  inertial  frames.  They  believe  that  the  existence  of  inertial  frames  is  the  only  way 
to  resolve  the  paradox.  Hence  the  mass  center  of  the  cosmos  is  the  ultimate  inertial  frame  and  the 
analysis  is  complete. 

In  summary,  this  section  has  shown  that  the  only  motion  of  the  earth  that  disqualifies  it  as  an  inertial 
frame  is  its  rotation  on  its  axis.  Frames  fixed  to  the  earth,  the  sun,  the  galactic  center,  and  the  center 
of  the  cosmos  are  ail  inertial.  Looking  as  far  as  the  galactic  center,  and  ignoring  the  earth’s  spin, 
observers  on  the  earth  are  in  an  elevator,  in  an  elevator,  in  an  elevator,  and  all  three  of  fliese 
elevators  are  inertial  frames  because  they  are  in  gravitational  free-fall  with  respect  to  each  other. 

It  is  not  the  case  that  a  nonrotating  earth  frame  is  inertial  because  the  acceleration  of  the  earth  with 
respect  to  the  sun  is  very  small.  The  qualification  as  inertial  is  independent  of  the  magnitude  of  the 
acceleration.  The  earth  could  be  spuming  around  the  sun  1000  times  a  second  at  a  radius  if  one 
meter  and  it  would  still  be  an  inertial  frame  if  it  did  not  rotate  on  its  own  axis.  On  such  a  planet,  all 
stars  but  the  sun  would  be  fixed  in  the  sky. 

^  Effetf  ftf  ths  Eaith's  RoUitigii 

In  order  to  investigate  the  effect  of  the  earth’s  rotation  on  an  accelerometer,  it  will  first  be  shown 
that  objects  do  not,  in  fact,  fall  in  straight  lines  on  earth.  Suppose  an  observer  attempts  to  reproduce 
Galileo’s  original  experiment  by  dropping  masses  from  the  top  of  a  tower  to  observe  their  motion. 
Assume  the  observer  is  positioned  on  the  equator  and  is  facing  east 

It  has  been  shown  that  the  object  is  attracted  by  the  sun  just  as  the  earth  is,  so  that  does  not  enter 
the  picture.  It  is  not  necessary  to  consider  what  happens  from  a  sun  fixed  viewpoint  or  from  a 
gala^  fixed  viewpoint  A  nonrotating  geocentric  frame  of  reference  is  inertial. 

In  this  inertial  frame,  and  not  in  the  rotating  earth  frame,  the  released  object  will  move  under  the 
influence  of  all  external  forces  arising  from  interactions  with  material  bodies  according  to 
Newton’s  second  law.  A  careful  exploration  of  the  consequences  of  this  reveals  some  very  startling 
results. 


23.  At  S  billion  years  old,  the  sun  has  only  ortrited  the  galactic  center  25  times  since  it  formed. 

24.  Hiis  is  die  dof^ler  shift  of  the  fpeqoeocy  or  color  of  light  due  to  relative  motion  of  emitter  and  receiver. 
Although  the  spe^  of  light  in  vacuum  is  constant  to  all  observers,  its  color  is  not  These  ted  shifts  were  an 
important  experimental  verification  of  the  fiact  that  the  universe  is  expanding. 
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4.2.1  Tangential  Velocity 


Establish  a  tangential/radial  cartesian  coordinate  system  to  be  oriented  according  to  the  direction 
of  motion  of  the  object  at  the  instant  it  is  released.  The  radius  of  the  earth  is  about  4000  miles,  and 
it  rotates  about  once  per  24  hour  day.  Therefore,  at  the  equator,  the  tangential  velocity  as  it  spins 
on  its  axis  is: 


V,  =  R,x£J^ 

1  2jc 

V  =  4  X  10  miles  x  ®  lOOOmph 

t  24hours 


Given  that  it  does  not  fall  straight,  it  may  seem  that  once  the  observer  releases  the  object,  it  will 
remain  stationary  in  the  inertial  frame  along  the  tangential  direction,  and  accelerate  along  the  radial 
direction,  so  the  planet  will  continue  to  rotate  around  it  until  it  reaches  the  ground  somewhere 
“behind”  the  observer.  It  actually  falls  in  front  of  the  observer^^. 


The  reason  is  that  it  remains  unaccelerated  and  not  stationary  in  the  tangential  direction.  So  it 
continues  to  translate  along  the  tangential  axis  at  its  original  considerable  speed.  Now,  because  of 
its  greater  distance  from  die  center  of  the  earth,  that  speed  is  in  fact  slightly  larger  than  that  of  the 
base  of  the  observer’s  tower,  which  is  fixed  to  the  ground.  The  object,  no  longer  constrained  by  the 
observer’s  hand,  moves  at  a  velocity  determined  by  its  slightly  larger  radius,  and  hence  moves 
faster  in  the  tangential  direction.  In  ^e  radial  direction,  it  is  accelerated,  so  the  entire  trajectory  is 
ballistic  in  the  inerfial  frame.  The  situation  is  depicted  in  Figure  6. 


25.  IbeTbwerof  Pisais  180  ft  high.  As  aresolt  ignoting  air  fiictkm,  all  of  Galileo’s  masses  would  have 
fallen  about  03  inches  to  the  east  of  the  point  diiectly  below  where  he  drcqn>6d  diem.  The  Enqxze  Stiue  Build- 
ii^  is  1200  ft  high,  giving  rise  to  a6  inch  eastward  deflection. 
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A.2J.  Coriolis  Force 


The  conclusion,  which  makes  perfect  physical  sense,  defies  intuition.  The  observer  is  not  in  an 
inertial  firame  and  observes  paradoxical  motion  if  he  assumes  that  he  is.  Often,  the  assumption  is 
mostly  valid  because  the  effect  is  relatively  small. 

From  the  point  of  view  of  the  observer  in  his  rotating  frame,  there  is  some  apparent  force  which 
has  caused  the  object  to  drift  to  the  east  This  force  is  called  the  Coriolis  force  .  It  is  one  of  two 
such  apparent  forces  which  arise  when  a  reference  frame  has  an  angular  velocity  with  respect  to  an 
inertial  frame.  The  other  is  called  centrifugal. 

For  an  accelerometer  in  an  INS,  which  measures  the  inertial  acceleration  ,  the  Euler  apparent 
force  is  zero  because  the  earth  rotates  at  constant  rate^.  The  geocentric  and  spinning  frames  are 
associated  with  the  same  object  (the  earth)  so  they  have  no  translational  interframe  position, 
velocity,  or  acceleration.  The  Einstein  apparent  force  is  zero.  Rewriting  the  earlier  acceleration 
relationship: 

if  =  +  2c5  X  +  ©  X  ®  X  f  _ 

t  m  m  m 


So,  even  when  the  vehicle  is  fixed  on  the  earth,  since  rQ|  is  measured  firom  the  center  of  the  earth, 
the  vehicle  sees  a  centrifugal  acceleration,  and  when  it  moves,  it  also  sees  a  coiiolis  acceleration. 
Accounting  for  these  apparent  forces  is  a  central  issue  in  inertial  navigation. 


26.  Also  called  sideways  force.  Visit  tbe  Faacanlt  pendolom  at  the  UN  building  in  New  Yodc  to  see  a  demon- 
sttadOD  of  it  in  action. 

27.  Not  really,  see  next  sectkm. 

28.  flooring  precession. 
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PART  Il:Inertial  Sensors  and 
Stable  Platforms 


1.  Accelerometers 

The  accelerometer  is  the  basic  sensor  of  inertial  navigation.  This  section  will  examine  its 
theoretical  operation,  its  limitations,  and  its  practical  implementation. 

IL  GraYity 

Gravity  is  defined  as  the  force  per  unit  mass  required  to  keep  a  test  mass  in  the  same  position 
relative  to  the  earth.  It  has  a  nominal  value  of  32.2  ft/sec/sec  or  9.8  meters/sec/sec.  The  practice  of 
measuring  and  modelling  gravity  is  called  geodesy.  In  the  last  section,  it  was  shown  that  an  object 
fixed  to  the  surface  of  the  earth  experiences  a  centrifugal  force  as  viewed  from  a  reference  firame 
spinning  with  the  earth.  A  plumb  bob  at  the  surface  of  the  earth  does  not  point  toward  its  center  but 
is  rather  displaced  slightly  toward  the  equator. 

The  natural  trajectory  for  an  object  positioned  at  a  radius  of  R,  given  some  initial  tangential 
velocity,  is  in  an  orbit  around  the  plank.  Its  real  traiectory  is  not  a  circle  centered  at  the  center  of 
the  earth  (an  orbit)  but  a  circle  of  constant  latitude  .  Some  force  other  than  gravitation,  imparted 
by  contact  with  the  surface  of  the  earth  must  be  the  cause  of  this. 

It  is  the  practice  in  geottesy  to  combine  both  gravitational  attraction  and  centrifugal  force  together 
into  a  single  force,  divided  by  the  mass,  and  call  the  result  gravity.  Local  density  changes  in  the 
earth  give  rise  to  local  changes  in  the  gravitational  attraction^,  and  it  also  varies  obviously  widi 
elevation.  The  centrifugal  force  varies  with  latitude  as  shown  in  Figure  7.  Henceforth,  gravity  and 


29.  TUs  is  “geostatiooai/’,  or  “geosynctirQOQas”  otbits  must  be  equatorial  and  die  satellite  pictnies 

in  yov  nigbdy  newscast  are  alw^  sboki  firom  diove  die  equator.  Geosyncbronoos  oilrit  over  Lonkn  or 

Yoik  is  inqmssible  widiont  great  e:q;)enditnre  of  energy,  despite  the  pcqnlar  science  fiction. 

30.  dmsity  meteorites  can  be  detected  in  this  way. 
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I2u  Spttific  Force 

Imagine  an  acceleiometer  constructed  from  a  test  mass  on  a  spring.  It  is  intuitively  clear  that  such 
a  device  must  measure  the  influence  of  gravitation  because  it  will  deflect  with  the  weight  of  the 
mass.  It  is  not  enough  to  simply  calibrate  this  constant  out  of  the  device,  because,  when  mounted 
on  a  vehicle  climbing  through  the  atmospb^  at  constant  speed,  the  mass  would  slowly  rise,  since 
gravitatitm  varies  with  distance,  compressing  the  spring  and  indicating  an  erroneous  acceleration 
upward.  This  is  yet  another  example  of  the  princ4)le  of  relativity,  for  an  observer  in  such  a  vdiicle 
cannot  tell  he  is  in  a  gravitational  field.  This  is  the  fundamental  reason  why  an  INS  can  only  operate 
in  a  known  (or  known  to  be  insignificant)  gravitational  field,  and  why  space  vehicles  using  an  INS 
must  know  the  gravitational  field  of  any  planets  they  encounter  to  avoid  missing  the  real 
accelerations  that  they  cause. 


Provided  a  vehicle  always  operates  near  the  surface  of  the  earth,  the  gravitational  constant 
acceleration  can  in  principle,  be  calibrated  out  However,  as  a  practical  matter,  real  vehicles  rotate 
so  the  direction  of  the  gravitational  attraction  vector  is  not  constant  in  the  vehicle  frame  of 
reference.  Some  systems  however,  completdy  avoid  the  gravitation  problem  by  maintaining  die 
accelnometer  sensitive  axis  nonnal  to  the  local  gravitation  vector.  The  real  issue  is  that  the 
deflection  of  the  mass  is  directly  influenced  by  flie  total  force  on  the  mass  and  not  by  its 
acceleration.  Whmi  stationary  in  a  gravitational  field,  a  mass  can  experience  force  even  when  there 
is  no  acceleratioii.  When  flee-falling,  a  mass  can  register  no  force  even  when  there  is  an 
accderatioiL  A  quantity  called  the  qiedfic  fmxe  is  defined  as  the  total  force  divided  by  the  mass. 


H  Effect  of  the  Barthes  Gw 


It  was  shown  in  die  last  cluqtter  that  a  nonrotating  geocentrk  frame  of  refoence  is  inertial.  This 
allows  the  expressitm  of  Newton’s  laws  in  such  a  frame.  However,  it  is  now  clear  that  choosing  an 
inertial  frame  is  only  half  the  problem.  It  is  also  necessary  to  know  the  value  of  the  gravitation 
vector  field  in  which  the  device  operates.  In  general,  diis  ^d  is  determined  by  the  superposition 
of  all  influences  of  all  bodies  in  the  universe.This  can  be  sq>proximated  easily  by  considering  only 
those  bodies  uduch  exert  significant  gravitational  force  on  die  test  mass. 

Gravitatitmal  force  is  proportional  to  the  mass  on  which  it  is  exerted,  so  it  is  completely  equivalent 
to  speak  of  an  acceloation  due  to  gravitation.  It  is  clear  that  when  operating  near  the  earth,  its 
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gravitation  is  significant,  as  is  the  gravitation  of  the  sun  or  moon  or  planets  when  operating  close 
to  than.  However,  on  the  earth,  the  sun  and  moon  give  rise  to  the  tides,  so  they  are  not  completely 
insignificant  forces. 

It  may  seem  that  ttese  effects  have  already  been  discounted  by  noting  that  the  earth  is  in  free-fall 
with  respect  to  them,  but  this  simply  established  the  earth  an  inertial  frame.  In  any  inertial  frame, 
teal  forces  arising  from  material  bodies  must  be  considered  in  the  equations.  The  i^xt  step  is  to 
account  for  all  those  real  forces,  and  gravitation  of  the  sun  and  the  moon  and  everything  else  in  the 
universe  are  real  forces.  If  they  were  not,  the  sun  would  not  orbit  the  galactic  center  once  every  200 
million  years  or  so. 

Luckily,  the  acceleration  due  to  the  sun  is  about  0.0006  g  and  that  of  the  moon  is  similarly  small, 
so  they  can  be  neglected  as  engineering  approximtuions  when  operating  near  the  earth.  The  forces 
from  all  other  matter  in  the  galaxy,  and  of  other  galaxies  can  likewise  be  neglected,  because, 
massive  though  they  are,  they  are  also  very  far  away. 

lA.  Basic  Equation  <if  Inertial  Navigatkm 

The  firame  of  reference  problem  is  finally  solved  in  principle.  The  basic  equation  of  inertial 
navigation  is  the  qjecific  force  equation  which  relates  the  accelerometer  readout  to  the 
gravitational  force  on  the  vehicle  and  its  inertial  acceleration.  Consider  a  spring  and  mass 
accelerometer  mounted  to  a  vehicle  chassis.  The  spring  expoknces  forces  T  required  to  counteract 
gravity  when  the  vehicle  is  stationary,  and  further  forces  as  die  vehicle  executes  accelerated  motion 
dragging  die  mass  along  with  it  Using  a  subscript  of  i  to  indicate  a  geocentric  inertial  frame,  T  to 
rqitesait  total  force  registered  by  the  accelerometer,  G  to  represent  gravitational  attraction,  and  m 


It  is  important  to  avoid  aMuming  a  constant  acceleration  due  to  gravitation.  The  accelerometer 
measures  the  specific  force  on  tte  vehicle  expressed  in  an  inertial  frame  in  vehkle  coordinates. 
Usually,  die  quantity  of  interest  is  the  position  on  the  surface  of  earth  and  hoice  in  the  rotating  earth 
frame.  Using  die  theorem  of  Coriolis: 

1=  (a,  +  2Sx^.  +  e5xi5xf_) -I-- 
m  ^ ® _ g  m _ 

This  equation  can  be  solved  for  die  acceleration  in  the  spinning  earth  frame  a^,  and  it  can  be 
integrated  twice.  The  earth’s  sidereal  rate  of  rotation  is  known.  The  position  and  velocity  can  be 
integrated  forward  from  die  initial  conditions,  and  there  are  standard  formulae  for  die  earth’s 
gravitation.  In  order  to  integrate,  all  accelerometa  readings  must  be  converted  to  a  common 
coordiiiate  system.  In  order  to  solve  this  coordinate  system  problon,  the  discussion  will  tom  to  the 
sectmd  inertial  sensor,  the  gyroscope,  in  the  next  section. 
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All  accelerometers  operate  similarly,  though  details  of  construction  diner.  A  pracucal 
accelerometer  is  a  precision  instrument  which  couples  a  mass  to  the  instrument  case  through  an 
elastic,  viscous,  or  electromagnetic  restraint  Typically,  the  mass  is  allowed  a  single  degree  of 
freedom  which  may  be  linear  or  rotary.  Rotary  deflecting  devices,  called  pendulous  devices,  are 
often  used.  Calibration  of  the  restraint  whatever  its  form,  provides  a  measure  of  the  specific  force 
that  is  tending  to  cause  a  deflection  along  the  free  degree  of  freedom. 

Most  devices  operate  in  a  closed  loop,  called  a  rebalance  loop,  which  prevents  the  mass  from 
actually  moving.  In  this  case,  the  effort  requited  to  retard  motion  is  a  measure  of  the  specific  force 
experienced.  Ttese  devices  are  generally  superior  to  those  permitting  error  inducing  motion.  Most 
devices  incorporate  fluids  to  damp  vibrations  or  to  support  critical  components  through  buoyancy. 

The  (tesign  of  inertial  grade  accelerometers  is  concerned  with  achieving: 

•  null  repeatability,  5  x  lO"^  g  is  typical 

•  low  threshold  (sensitivity),  10*^  g  is  typical 

•  linearity,  0.10%  over  10  g  is  typical 

•  large  dynamic  range 

•  time  constant,  few  milliseconds  is  typical 

•  low  sensitivity  to  vibration 

Accelerometers  may  be  designed  to  be  mounted  on  a  stable  platform.  These  incorporate  their  own 
small  gyroscopes  and  perform  integration  as  part  of  their  function.  The  string  accelerometer 
maintains  mass  between  two  strings  under  tension.  Beat  frequencies  between  the  vibrating  strings 
provide  the  output  There  are  many  exotic  types  including,  nuclear,  electrostatic,  cryogenic, 
gaseous,  particle  stream,  solid  state,  and  vibrating  reed  instruments. 


Modem  Ineftial  and  Sattilite  Navigaikm  Systems 


page2S. 


2  Mechanical  Gyroscopes 

While  the  accelerometer  measures  translational  motion,  the  gyroscope^^  or  gyro  is  a  device  used 
to  measure  or  maintain  orientation.  Cartesian  coordinate  systems  are  defined  by  the  orientation  of 
their  reference  axes,  so  the  gyroscope  can  be  used  to  solve  the  coordinate  system  problem.The 
properties  of  die  gyroscope  are  based  on  the  spinning  rotor  which  was  fijrst  investigated  by  Euler 
and  Faucault  in  the  eighteenth  century.  The  two  primary  properties  of  interest  to  inertial  navigation 
are  the  rigidity  of  orientation  in  inertial  space  and  precession. 

2,L  Simple  Gyroscope 

A  stable  angular  reference  might  be  constructed  based  on  the  rigidity  property  of  the  gyro  as 
follows.  Use  a  symmetric  body  (disk  or  ring),  spin  it  about  its  axis,  and  isolate  it  firom  external 
torques.  The  method  used  to  isolate  the  spun  mass  from  external  torques  is  to  use  a  system  of 


As  the  platfmm  is  rotated,  the  gyro  will  tend  to  remain  stably  pointed  in  inertial  space.  Sensors 
mount^  on  the  gimbals,  as  indicated  schematically  in  the  figure,  can  measure  the  vehicle  rotation 
in  this  way. 


31.  Fram  the  Greek  gyros  for  "code”  or  “lotaiion”,  and  skopetn  *^0  view”.  Faucault  used  the  gyroscope  to 
view  the  roittk»  of  earth  by  pointing  the  device  at  a  star. 
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Real  gyros  are  built  based  on  the  second  propemr  of  the  spinning  rotor,  precession,  and  employ 
^rvomechanisms  to  implement  a  stable  table^  on  which  the  accelerometers  are  mounted. 
Consider  Figure  9: 


Rgure9  Single  Axis  Gyro  Precession 

Any  attempt  to  yaw  the  platform  on  which  the  gyro  is  mounted  about  die  input  axis  will  give  rise 
to  a  nodding  of  the  device  about  die  ouQiut  axis.  A  rotational  transducer  such  as  a  resolver  placed 
on  the  output  axis  will  register  this  rotation  which  is  proportional  to  the  applied  torque.  Using  this 
technique,  gyros  can  be  employed  as  null  indicating  devices  in  tibe  feedbacJc  loops  of 
servomechanisms. 

Single  axis  gyros  such  as  this  one  are  caUed  restraiiied  gyros  since  they  are  prevented  from  motion 
about  all  axes  but  one,  and,  in  practice,  the  ouqiut  angles  are  never  allowed  to  become  large.  There 
are  also  two  degree  of  freedom  gyros  in  common  use.  These  are  called  firee  gyros,  incorporate  an 
extra  set  of  gimbals  and  ate  used  to  control  angular  velocity  around  two  orthogonal  axes.  Gimbals 
may  be  used  both  in  the  gyro  itself  and  in  the  stable  table,  so  sometimes  it  is  necessary  to 
distinguish  these. 


32.  A  stable  table  is  staUe  in  orieotattoD  in  a  specific  condinate  system. 
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Practical  Gyroscopes 

Gyros  can  be  used  as  rate  gyros  or  as  inte^ting  gyros.  The  fonner  are  characterized  by  spring 
restraint  of  the  output  axis  and  low  damping  and  derive  their  name  from  the  fact  that  the  output 
angle  is  proportional  to  the  input  rate.  These  are  useful  as  turn  indicators,  but  are  not  used  in 
precision  inertial  systems^^.  The  latter  are  characterized  by  no  spring  restraint  and  much  higher 
damping  of  the  output  axis,  and  derive  their  name  from  the  fact  that  the  ouQ)ut  angle  is  proportional 
to  the  integral  of  the  input  rate. 

Typical  requirements  for  a  gyro  are  the  measurement  of  one  thousandth  of  a  degree  per  hour  of 
angular  motion.  The  ability  to  manufacture  precision  gyroscopic  instruments  which  are  sufficiently 
free  of  unpredictable  drift  is  the  ultimate  limiting  factor  in  the  performance  of  inertial  navigation 
technology.  Gyro  drift^  is  the  overriding  performance  concern  and  it  is  caused  by  parasitic 
precessional  torques  which  are  outside  the  system  mechanization  loops. 

Gyroscopes  develop  imperfect  performance  for  all  of  the  following  reasons: 

•  mass  imbalance  of  the  rotor 

•  thermal  gradients 

•  anisoelasticity 

•  buoyancy  imbalance  (in  floated  devices) 


Current  technology  relies  on  the  following  techniques  to  improve  performance: 

•  flotation  in  fluid  to  buoy  components  and  damp  oscillations 

•  high  incidence  of  cylindtical  symmetry  in  components 
■  precise  manufacturing  tolerances 

•  matched  coefficients  of  thermal  expansion 

•  temperature  control 

•  {Heloaded  beatings 


Some  practical  designs  of  gyroscopes  include  the  floated  gyro,  originally  developed  at  MTT,  air 
bearing  gyros,  electrostatic  gyros,  and  optical  gyros.  All  but  the  latter  derive  their  name  from  the 
technique  used  to  isolate  the  rotary  bearings  from  loads.  This  is  because  spin  axis  bearings  are  an 
important  limiting  factor  in  performance. 


33.  See,  bowevei;  sttapdown  systems  later. 

34.  The  earth  rotates  IS  d^iees  per  boor.  Drift  is  often  measored  in  mim  earth’s  rate  units  or  “meius”  equal 
to  0.0015  d^iees  per  hour. 
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3  Stable  Platforms 

In  the  stable  table  or  stable  platform,  a  feedback  loop  is  used  to  regulate  the  table  in  a  such  a  way 
as  to  prevent  the  precession  of  the  gyroscopes.  This  important  principle  is  the  basis  of  maintaining 
a  stable  angular  reference  in  inertial  systems. 

^  Sin^  Axis  Stable  Table 


The  figure  uses  the  standard  notation  for  the  input  (lA),  output  (OA),  and  spin  axes  (SA)  of  a 
gyroscope.  The  directions  of  input  torque  (N),  angular  momentum  (L),  and  precession  (Q),  are 
associate  widi  these  axes  as  shown  in  the  figure. 

An  intuitive  explanation  of  the  operation  of  the  platfonn  is  as  follows.  Suppose  the  rotation  of  the 
ouq>ut  axis  is  measured  and  used  to  drive  an  amplifier  coimected  to  a  motor  which  is  able  to  rotate 
the  platform.  If  the  vehicle  to  which  the  rectangular  base  is  attached  performs  a  yawing  motion, 
this  is  a  rotation  about  die  input  axis  (lA)  of  the  gyro.  Liitially,  the  motor  will  trananit  this  torque 
directly  through  to  the  circular  table.  This  will  yaw  the  table  a  small  amount  which  will  cause  the 
circular  table  to  rotate.  This,  in  turn,  will  cause  the  gyro  to  nod  about  its  output  axis  (OA).  If  this 
ouqmt  signal  is  fed  back  to  the  motor,  a  suitable  servo  control  law  will  drive  the  circular  table  in 
such  a  way  as  to  null  the  gyro  output 

The  conclusion  is  that  the  platform  will  be  forced  to  execute  the  correct  motions  to  keep  the  gyro 
from  nodding  about  its  output  axis.  From  the  point  of  view  of  the  vehicle,  it  performs  a  yawing 
motion,  and  the  table  stays  behind.  The  circular  platform  stays  where  it  is  pointed  in  an  inertial 
frame  of  refermice. 
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^  iartt  AXIS  attPifi-iaBK 

It  should  be  clear  that,  because  the  accelerometers  are  mounted  on  the  stable  platform  with  their 
sensitive  axes  oriented  along  three  orthogonal  directions,  the  orientation  of  the  stable  platform 
defines  the  coordinate  system  in  which  the  inertial  specific  force  is  measured.  It  is  the  function  of 
the  stable  table  to  ensure  that  the  accelerometers  always  point  in  the  direction  required  by  the 
calculations. 

Practical  stable  tables  always  incorporate  at  least  three  gimbals  and  often  more.  Torque  motors, 
called  torquers,  and  rotary  position  sensors,  called  signal  generators  or  pickoffs  are  mounted 
directly  to  the  gimbals.  For  earth-fixed  systems,  the  torquers  are  also  used  to  rotate  the  table  in 
order  to  maintain  a  locally  level  configuration. 

Such  tables  directly  implement  the  Euler  angle  model  of  rotation.  This  has  two  important 
consequences  to  the  design  of  a  three  axis  stable  table.  First,  because  the  rotation  axes  of  the 
gimbids  are  not  always  orthogonal,  it  is  necessary  to  uncouple  the  control  loops  driving  the  gimbal 
torquers.  Second,  the  gimbals  are  subject  to  the  inherent  singularity  of  such  an  arrangement  In  the 
context  of  stable  tables,  the  condition  of  singularity^^  is  called  gimbai  lock,  and  it  occurs  when  all 
three  axes  of  rotation  fall  in  the  same  plane.  In  this  condition,  the  apparatus  looses  one  of  its 
degrees  of  freedom. 

Gimbal  lock  can  be  avoided  in  two  ways.  For  ships  and  similar  vehicles,  the  required  angular  travel 
of  the  table  with  respect  to  the  vehicle  is  small  enough  that  the  arrangement  of  gimbals  can  be 
designed  to  avoid  gimbal  lock.  For  guided  missiles  and  other  vehicles  which  require  unconstrained 
rotation,  a  fourth  gimbal  can  be  added. 

Ctfmptete  IWS 

^th  a  three  axis  stable  table,  the  coordinate  system  problem  is  solved.  The  gyros  can  be  mounted 
on  the  table  in  such  a  way  as  to  provide  orientation  information  around  three  orthogonal  axes.  The 
accelerometers  are  mounted  on  the  same  table  in  such  a  way  that  their  sensitive  axes  are  also 
oriented  along  three  orthogonal  axes.  This  alignment  permits  measurement  of  the  specific  force 
vector  in  three  dimensions  in  an  ineitially  stable  coordinate  system. 

Another  frequently  used  iq)ptoach  is  to  deliberately  drive  the  platform  to  rotate  as  the  earth  does, 
so  that  an  earth  toted  coordinate  system  is  used,  and  motion  relative  to  the  earth  is  measured. 
However,  in  this  case,  it  is  only  the  coordinate  system  which  is  earth  referenced.  The 
^elerometers  will  still  register  the  cotiolis  and  centrifugal  forces  because  specific  force  is 
inherently  an  inertially  derived  vector,  and  no  technology  can  change  this  fact  Simply  changing 
coordinates  cannot  eliminate  the  effect  of  rotiuing  at  1000  m.p.h.  around  the  center  of  the  earth. 
Here  the  importance  of  the  distinction  between  reference  frame  and  coordinate  system  is 
illustrated. 


35 .  Singularity  is  a  commoo  concetc  in  in«»rh«nisms  It  can  be  defined  in  many  ways:  loss  of  a  degree  of  free¬ 
dom,  matfaematkal  singolaiity  of  the  Jacobian  matrix,  co-linearity  of  axes  or  linear  dqteodoice  of  axes,  mtd 
coDvogence  of  redundancy  to  a  single  point  All  definitions  are  equivalent  Singularity  is  always  a  very  se- 
tkms  matter  smce  it  implies  loss  of  a  degree  of  fieedcmi  fin’  posirioning  and  infinite  rates  for  rate  controL  The 
Euler  wrist  commonly  used  in  robot  mampulators  suffos  from  the  same  problems. 
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4  Gyrocompasses 

Although  it  is  typically  not  a  component  of  the  inertial  navigation  system,  the  gyrocompass  is  an 
ingenious  device  based  on  many  of  the  same  principles.  The  gyrocompass  has  replaced  the 
magnetic  compass  on  virtually  all  large  ships  and  it  is  not  based  on  magnetism  at  all.  The  device 
can  be  used  to  implement  an  automatic  helmsman  on  a  ship  by  using  it  as  a  null  indicator  in  a  servo 
loop.  The  magnetic  compass  has  two  disadvantages  compared  to  the  gyrocompass.  First,  the  earth’s 
magnetic  field  varies  locally  and  the  poles  themselves  move  slightly  over  tu?.e.  Second,  the 
presence  of  magnetic  material  on  a  vessel  disrupts  the  operation  of  the  magnetic  compass.  The 
gyrocompass  is  immune  to  these  two  problems. 

It  will  be  shown  that  the  equation  of  motion  of  this  device  is  that  of  a  pendulum  where  the  spin  axis 
oscillates  about  the  local  north-south  direction.  Damping  this  oscillation  causes  the  device’s  spin 
axis  to  eventually  settle  into  north-south  alignment.  In  this  way,  it  functions  as  a  compass.  Several 
means  of  damping  have  been  used,  but  one  of  the  simplest  and  best  methods  is  a  eccentrically 
weighted  semicircular  ring  developed  by  Elmer  A.  Spenry^^. 

^  Equation  of  Motion 

The  theory  of  operation  of  the  device  can  be  derived  from  Euler’s  equation.  The  device  is 
constructed  firom  a  rotor  spinning  at  veiy  high  speed.  The  rotor  is  constrained  so  that  its  axis  must 
stay  in  a  plane  which  is  horizontal  at  its  current  position  on  the  earth.  That  is,  the  normal  vector  to 
the  plane  is  the  local  vertical.  This  rotor  is  mounted  inside  a  set  of  gimbals  which  pemit  rotation 
about  the  local  vertical. 

To  derive  the  equation  of  motion,  consider  Figure  11.  The  right  figure  is  a  view  from  the  equatorial 
plane  of  the  earth  spinning  on  its  axis.  The  left  view  is  a  local  overhead  view  of  the  device  with  its 
spin  axis  displaced  from  the  local  north-south  direction. _ 


Rgure  11  Gyrocompass 

Let  the  device  be  positioned  at  some  latitude  X  on  the  rotating  earth.  There  are  three  independent 
rotations  happening  simultaneously.  The  earth  rotates  about  its  axis,  the  device  outer  gimbal  rotates 


36.  Of;  you  guessed  it,  Sperry  Gyrosayw  CcHopany. 
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about  its  local  vertical,  and  the  device  itself  rotates  at  a  high  rate  about  its  symmetry  axis.  Let  the 
gimbal  rotation  be  given  by  0.  The  discussion  will  be  concerned  with  four  different  frames  of 
reference  and  their  relative  motion.  A  frame  positioned  at  the  center  of  the  earth  which  is  not 
rotating  is  inertial  (subscript  i).  The  spinning  earth  frame  will  have  subscript  e.  The  gyrocompass 
is  positioned  on  the  rotating  earth  and  rotates  about  the  local  vertical.  A  frame  of  reference  attached 
to  the  compass  but  not  spinning  around  its  symmetry  axis  will  be  called  the  rotating  frame 
(subscript  r).  Finally,  a  frame  attached  to  the  spinning  rotor  will  be  called  the  spinning  frame 
(subscript  s). 

The  problem  is  solved  by  writing  Euler’s  equation  in  the  rotating  frame.  Although  it  is  not 
spinning,  this  frame  is  chosen  because  the  inertia  tensor  is  still  constant.  Also,  it  is  the  frame  whose 
motion  will  ultimately  prove  interesting.  A  cartesian  coordinate  system  for  this  frame  is  given  in 
the  figure.  The  angular  velocity  of  each  noninertial  frame  with  respect  to  the  inertial  frame  will  be 
expressed  in  terms  of  the  angles  defined  in  Figure  11  and  expressed  in  the  chosen  coordinate 
system.  In  doing  so,  the  fact  that  angular  velocities  add  like  vectors^^  will  be  used. _ 

rotation  of  earth  frame  in  inertial  frame 
rotation  of  rotating  frame  in  inertial  frame 
rotation  of  spinning  frame  in  inertial  frame 
rotation  of  spinning  frame  in  rotating  frame 

Now  all  of  these  are  expressed  in  the  rotating  frame  of  reference. 

©-  =  (iisinXi  +  GicosXsinQ}  +  (OcosXcosQk 

w  w  w  w 

©s  =  ©j.  +  i2k  =  (©gjmX.  +  0)i +  ©gC(75A,5m6j  +  (©gC<7.yA,c<7j0  +  n)k 


Substitution  of  these  complicated  expressions  will  be  avoided  in  the  following  development 
Consider  now  Euler’s  equation  applied  to  the  rotor.  The  theorem  of  Coriolis  will  be  used  where  the 
inertial  finune  is  the  fixed  frame  and  the  rotating  frame  is  the  moving  one.  The  theorem  takes  the 


form:  ____ 
N  = 


Notice  that  the  two  angular  velocities  of  the  cross  product  are  now  different  This  is  because  the 
first  angular  velocity  is  the  relative  rotation  of  the  two  frames  and  comes  from  the  theorem  of 
Coriolis  while  the  second  is  the  rotation  of  die  body  in  inertial  space  because  it  is  related  to  angular 
momentum. 


37.  While  this  is  true,  angles  tbemsdves  do  not  add  like  vectors.  Tb^  are  not  a  vector  space  in  the  linear  alge¬ 
bra  sense.  This  is  the  basic  reason  why  inelegant  rotation  matrices  are  used  to  rqnesent  totaricmal  displace- 
moits,  while  everydiing  else  in  simple  mechanics  is  generally  a  vector. Angles  can  be  made  to  add  scnnewhat 
like  vectms  osing  a  generaliziakn  of  the  ctanplex  numbm'  called  the  quatemicm. 
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Let  denote  the  moment  of  inertia  about  the  symmetry  axis  z,  and  denote  the  moment  about  the 
other  two  axes.  Performing  the  cross  product  and  moving  subscripts  to  superscripts  to  allow 
indication  of  components  gives: 


1 :  N  =  I  cb*  +  (I  cp^  CO®  —  I  co^co* ) 

X  a  x  '  s  y  z  a  z  y^ 

2:N  =  I  (b®  +  (I  co^^co®  -  I  co^co®) 

3;N  =  I  cb®+  (I  (0*^(0® -1, 0)^(0®) 

z  sz'^axy  ayx'^ 


This  equation  looks  incorrect  because  it  is  not  the  standard  form  of  the  scalar  Euler’s  equation,  but 
this  arises  from  the  use  of  the  rotating  frame  as  the  moving  frame  and  not  the  spinning  frame. 
Consider  the  last  equation.  The  z  axis  is  gimballed,  so  it  can  exert  no  torque.  Also  since 
and  CO^y  =  <£fy,  the  bracketed  term  is  zero.  This  implies  that  (O^z^  constant.  This  constant  is 
composed  of  two  terms  proportional  to  the  earth’s  rate  of  rotation  and  the  rotor’s  rate  of  rotation. 
The  rotor  spins  at  many  orcters  of  magnitude  higher  velocity,  so  the  earth’s  rate  can  be  neglected. 
It  can  be  concluded  that: 


(0®  =  constant  =  co^co5A,co56  +  Q«Q  (since  Q»co^) 

Z  w  6 


Consider  the  first  equation.  The  x  axis  is  also  gimballed,  so  there  can  be  no  applied  torque.  Also 
the  second  term  in  brackets  is  of  order  and  can  therefore  be  neglected.  Rewriting:  _ 


I  (b®  =  -Lco’^Q 

ax  s  y 


L6  =  -l.((0.cosXsinQ)  Q, 


„  I  (CO 

e  +  -? — - sine  =  0 

^a 
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^  Interpretation  as  a  Pendulous  Oscillator 

Now  this  is  the  equation  of  a  pendulum.  The  variable  d  oscillates  like  a  pendulum,  or  equivalently, 
the  spin  axis  oscUlates  about  the  local  north-south  direction.  Addition  of  a  damping  system  will 
cause  the  device  to  settle  into  pointing  north. 

Notice  that,  like  a  magnetic  compass,  the  device  will  have  trouble  near  the  poies^^  because  of  the 
latitude  dependence,  and  that  the  rate  of  oscillation  is  proportional  to  the  spin  rate.  Also,  if  the  earth 
did  not  rotate,  the  device  would  not  work  at  all.  In  this  respect,  it  operates  as  a  compass  because 
the  earth  is  rotating.  It  is  also  interesting  to  note  that  the  device  can  be  used  to  determine  latitude, 
by  measuring  the  spin  rate  of  the  device  and  its  period  of  oscillation. 

The  operation  of  the  device  can  be  understood  on  an  intuitive  level  as  follows.  It  is  mounted  on  a 
spinning  earth^^,  so  it  is  experiencing  a  very  small  torque  transmitted  through  its  fixture  sufficient 
to  rotate  it  once  per  day  in  inertial  space.  This  torque  is  parallel  to  the  earth’s  axis  itself,  but  only 
that  component  of  it  which  is  oriented  along  the  y  axis  in  the  figure  can  be  transmitted  through  the 
gimbals  to  the  rotor.  The  device  is  a  gyroscope,  so  it  will  precess  according  to  the  gyroscope 
rigidity  equation,  but  this  precession  will  be  normal  to  both  the  spin  axis  and  the  axis  of  the  earth’s 
rotation.  The  kinematics  is  such  that  the  tendency  to  precess  is  opposite  to  the  current  direction  of 
rotation,  so  the  device  oscillates. 


Period  of  OsdUation 


Using  typical  figures  for  the  parameters  involved,  the  period  of  oscillation  can  be  calculated.  The 
inertia  ratio  I^/Ia  is  close  to  2.0  for  a  flat  circular  object  Let  the  latimde  be  the  equator,  and  the  rotor 
spin  rate  be  a  typical  value  of  3600  rpm  or  60  rev/sec: 


T  =  2jc 


I 


l.(GicosX)  Cl 


1/2 


'T  -  24  x  60  x  60 

“  4x27Cx27Cx60^ 


T  =  47  secs 


1/2 


A  very  high  rate  of  spin  is  required  in  order  to  get  a  reasonable  period  of  oscillation.This  is 
important  because,  depending  on  the  initial  conditions,  the  damping  cannot  bring  the  device  into 
alignment  any  faster  than  one  period,  and  a  device  which  takes  an  hour  to  settle  is  useless.  Notice 
that  even  for  a  rotor  spin  of  3600  rpm,  the  pendulous  oscillation  will  be  imperceptible  to  the  eye. 


38.  However,  the  gyroconqass  measures  true  north,  not  magnetic  nordi,  so  it  is  immune  to  the  complicated 
motkms  of  the  magnetic  pole. 

39.  This  is  the  fourth  time  in  this  rqxnt  Aat  a  phenomenon  that  measures  the  earth’s  rotation  has  occurred: 
the  Fancauh  pendulum,  die  eastward  deflection  of  a  dropped  object,  the  gyroscope  pointed  at  a  star,  and  the 
gyroconqiass. 
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5  Optical  Gyroscopes 


With  the  recent  arrival  and  rapid  improvement  of  optical  gyroscopes,  the  complete  disappearance 
of  mechanical  gyroscopes  in  most  applications  now  seems  inevitable  within  the  next  ten  years.  For 
this  reason,  these  devices  which  were  once  considered  impractical,  have  now  become  required 
knowledge  in  the  field  of  inertial  navigation.  This  section  will  briefly  overview  the  principles  of 
operation  of  these  devices. 

Current  commercial  optical  gyros  are  designed  to  be  strapped  down  directly  to  the  vehicle  - 
requiring  no  stabilized  platform  at  all.  In  this  respect,  they  complement  the  general  trend  toward 
strapdown  technologies.  Optical  gyros  are  fundamentally  solid  state  devices,  and  because  of  this, 
are  mote  reliable,  more  rugged,  and  are  able  to  withstand  high  g  fields  beyond  the  capabilities  of 
mechanical  devices.  They  can  be  designed  for  harsh  environments,  exhibit  large  dynamic  range^ 
and  exhibit  linearity  of  parts  per  billion  over  that  range.  These  devices  are  currently  in  use  on  the 
latest  civil  transport  aircraft,  including  the  Boeing  757,  and  767,  and  the  Airbus  A3 10,  and  have 
now  logged  millions  of  flight  hours. 

SJ.  Michelson  Interferometer 

In  the  period  shortly  before  relativity  theory  was  developed  by  Emstein,  physicists  believed  that 
light  propagated  through  an  invisible,  odorless,  tasteless,  medium  that  pervaded  all  of  space  called 
the  Ituniiiiferoiis  eAer.  In  1881,  A.  A.  Michelson  conducted  an  experiment  which  demonstrated 
this  idea  to  be  wrong  and  set  the  stage  for  the  development  of  the  special  theory  of  relativity^^. 
Michelson  was  the  first  American  to  be  awarded  the  Nobel  Prize  in  science. 

He  demonstrated  this  through  the  use  of  a  device  of  incredible  senativity  which  he  invented,  called 
the  Michelson  interferometer^^.  He  showed  that  the  “ether  wind”,  supposed  to  result  from  the 
earth’s  motion  through  space,  was  undetectable  by  devices  that  were  easily  sensitive  enough  to 
detect  it  if  it  existed. 


52.  SasnitfEl 


Consider  now  a  slightly  different  apparatus,  called  a  rii^  interferometer,  upon  which  optical 
gyroscopes  ate  based,  as  indicated  in  Figure  12.  If  the  device  is  made  to  simultaneously  propagate 
two  counterrotating  beams  of  coherent  light,  it  can  be  used  to  measure  rotation  according  to  an 
effect  discovered  in  1913  by  G.  Sagnac  and  known  today  as  the  Sagnac  effect  All  current  optical 
gyroscopes  are  based  on  the  Sagnac  effect 


40.  Dynaoiic  range  is  the  range  of  rotation  rai^  that  can  be  measured  cmrectly. 

41 .  AUioa^  Einst^  is  jusdy  credited  widi  the  thec^,  many  at  the  ideas  resulting  from  it  were  stated  earlier 
by  Lotents  and  also  by  Poincaure’ .  Einstein’s  contribotioa  was  to  show  that  all  of  the  fudge  facttxs  being  pro- 
posed  to  explain  Michelson’s  results  could  be  derived  fiom  a  very  small  set  of  simple  postulates. 

42.  Tb  diis  day,  devices  based  on  interferometry  are  among  the  most  precise  instruments  built.  This  arises 
firom  ibdr  abOity  to  measure  distance  on  the  st^  of  the  q;)ticai  wavdength. 
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5^1  Ring  Interferometer 


The  ring  interferometer,  here  composed  of  four  mirrors,  splits  the  beam  into  two  counterrotating 
components  and  they  ate  recombined  after  a  complete  transit  of  the  loop  to  observe  interference 
fringes.  Sagnac  showed  that  rotation  of  the  entire  loop  gives  rise  to  an  apparent  lengthening  of  the 
path  length  for  one  beam,  and  an  apparent  shortening  for  the  other. 

This  path  length  difference  gives  rise  to  motion  of  the  interference  fringes  which  is  a  direct  measure 
of  the  rotation  of  the  device  in  inertial  space.  Where  Michelson  showed  that  it  is  impossible  to 
measure  inertial  translation  using  light  because  the  speed  of  light  is  the  same  for  all  observers, 
whatever  their  state  of  motion,  Sagnac  showed  that  optical  measurement  of  inertial  rotation  was 


Figure  12  Ring  Interferometer 


There  are  two  intuitive  interpretations  of  the  effect,  based  on  different  propagation  path  lengths,  or 
different-times  of  transit  for  both  beams,  and  both  lead  to  the  same  conclusion,  namely  a  relative 
phase  shift  of  both  beams. 

A  rigorous  d^vation  of  the  effect  can  only  be  done  using  general  relativity,  because  in  the 
presence  of  a  gravitational  field,  even  special  relativity  caimot  be  used.  However,  there  is  a  classical 
derivation  due  to  Aronowitz,  which  gives  die  tight  answer  and  is  intuitively  clear.  To  simplify 
things,  consider  a  perfectly  circular  ting  interferometer  of  radius  r.  Let  the  beam  enter  at  point  A, 
be  split  into  a  clockwise  and  a  counterclockwise  beam,  complete  one  circuit,  and  be  recombined 
again  at  A.  Let  the  device  rotate  at  a  rate  Q  and  let  c  be  the  speed  of  light  Let  4  and  t  denote  the 
round  trip  transit  times  for  the  counterclockwise  and  clockwise  beams  respectively. 
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SJIJI  Apparent  Patii  Length  Difference 


Expressions  for  these  can  be  written  in  terms  of  themselves  by  noting  that  the  apparent  path  length 
difference  is  given  by  the  arc  length  travelled  by  the  interferometer  during  the  transit  time,  or  rflt 


27cr  +  rQt , 


Ikt  -  rQt 


U  = 


t  = 


Solving  for  the  times  and  taking  the  difference  gives: 


Wq 

At  =  t  -t.  =  ^ — 5-^ 


For  practical  values  of  r,  (rfli)^  «  c^.  so  the  expression  simplifies  to 


Thus,  the  round  trip  path  difference  is: 

AL  =  cAt  = 

c 


Real  devices  also  typically  use  a  square  or  triangular  loc^,  but  this  does  not  significantly  affect  the 
formula.  For  arbitrary  sh^)ed  devices,  the  circular  area  is  rq>laced  by  the  polygonal  area  and  the 
rotation  rate  is  projected  onto  the  normal  to  this  area.Unfortunately,  the  pad)  lengths  given  by  this 
formula  are  much  less  than  an  optical  wavelength,  ^^ch  should  be  clear  when  contemplating  how 
fast  light  traverses  a  distance  of  a  few  meters  . 

For  a  loop  area  of  one  square  meter,  and  a  rotation  rate  of  one  degree  per  second,  the  device 
measures  only  a  few  angstroms  of  deviation.  Real  tqtplications  require  much  smaller  sized  devices 
and  sensitivities  many  thousand  times  this  one  degree  per  second  figure.  For  diis  reason,  the  Sagnac 
effect  was  a  mere  scientific  curiosity  for  decades  after  it  was  discovered.  However,  as  a  result  of 
two  recent  technical  developments,  lasers,  and  fibre  optics,  two  different  ways  of  significantly 
increasing  the  sensitivity  of  the  effect  have  become  possible. 


I 


I 


43.  As  a  convenient  mnemonic,  light  tiavels  a  foot  in  a  nanosecond. 
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5.3  Lasers 

The  term  laser  stands  for  light  amplification  by  stimulated  emission  of  radiation.  The  laser  is 
basically  an  optical  oscillator.  As  such,  it  incorporates  all  of  the  features  common  to  an  electrical 
oscillator,  namely  filtering,  and  positive  feedback.  The  filtering  is  provided  by  the  fact  that  the 
resonance  cavity  sets  up  standing  waves  whose  wavelength  must  divide  die  cavity  length  by  an 
integer.  Positive  feedback  arises  from  the  use  of  a  lasing  medium,  a  gas,  liquid,  crystal  ions,  or  any 
of  a  number  of  other  possibilities.  The  lasing  medium  has  the  property  that,  under  the  conditions 
of  operation,  excited  atoms  can  be  stimulated  to  release  new  photons  of  radiation  of  precisely  the 
same  phase,  frequency,  and  directional  characteristics  as  the  energy  stimulating  the  atom.  Where 
there  was  one  photon,  a  new  i(tentical  copy  appears  as  well  Energy  is  continually  supplied,  or 
pumped,  to  the  system  to  excite  the  lasing  me^um  atoms. 

SA.  Riny  Laser  Gyroscopes 

The  ting  laser  gyro  (RLG)  combines  the  features  of  a  ring  interferometer  and  a  laser  into  the  same 
device.  Toe  lasing  medium  is  placed  directly  into  the  path  of  the  counterrotating  beams  and 
coherent  light  is  made  to  circulate  in  the  ting. 

Now,  since  the  laser  light  wavelength  is  related  to  the  cavity  length,  and  since  the  q)patoat  cavity 
length  is  different  for  both  beams  by  the  Sagnac  effect,  it  follows  that  the  frequency  of  the  two 
beams  will  slightly  differ.  For  a  path  length  of  L,  an  integer  m,  and  wavelength  X,  tite  resonance 
condition  can  be  expressed  as: _ 


mX  =  L 


Using  the  previous  result  for  die  Sagnac  effect,  the  frequency  difference  for  the  two  beams  can  be 
approximated  by: 


.  me  me  mcAL  ooAL 

Aa)  =  ©^-o) 

*  -  I-+  L  i}  L 


Since  light  frequencies  ate  extremely  high,  this  ftequmicy  differrace  is  laig^  dian  die  Sagnac  path 
differmtial  by  the  conaderable  factor  of  OO/L.  This  amplification,  combined  with  die  relative  ease 
of  measuring  the  beat  frequency,  typically  in  the  kHz  range,  makes  the  RLG  a  practical  device.  It 
is  important  to  recognize  that  the  RLG  measures  the  frequency  difference,  not  the  path  length 
difference  given  by  the  Sagnac  effect 

tl  Ring  Laser  Gyroswnts 


The  development  of  the  RLG  had  to  wait  for  the  laser  to  be  developed.  In  1962,  a  mete  two  years 
after  die  h^um-neon  laser  developed,  Macek  of  Sperry  Gyroscope  company,  produced  a  one 
meter  square  prototype  that  measured  the  earth’s  rotation  of  1/4  rad/hour.  Sagnac’s  original 
interferometer  had  a  sensitivity  of  about  2  rad/sec,  so  Macek’s  was  over  10,000  times  more 
smtative. 
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One  of  the  (tactical  problems  of  building  these  devices  is  that  the  two  beams  tend  to  compete  for 
excited  atoms,  and,  if  left  alone,  one  beam  or  the  other  will  eventually  win  the  competition,  and 
half  of  die  laser  will  be  shut  off.  One  method  available  to  counter  this  effect  is  to  use  two  different 
lasing  isotopes  which  are  individually  tuned  to  one  beam  or  the  other.  The  neon  laser  is  a  popular 
choice  for  RLCs  because  of  the  i^ar  ideal  separation  of  the  two  isotopes  neon-20  and  neon-22. 


A  more  important  practical  problem  is  the  iock4n  phenomenon,  which  is  a  property  of  all 
oscillators.  Imperfections  in  die  mirrors  used  in  the  laser  cause  a  portion  of  the  clockwise  beam  to 
be  reflected  badt  along  the  path  of  the  counterclockwise  beam  and  vice  versa.  This,  in  turn,  causes 
the  two  beams  to  t^proach  each  other  in  frequency  when  the  frequency  difference  is  small,  as  it 
will  be  for  small  rotational  rates.  The  net  effect  is  a  large  deadband  in  the  response  of  the  device. 
Rotational  rates  below  a  threshold  will  not  be  sensed  at  all  because  the  beams  will  lock  and  the  beat 
frequency  will  be  zero.  Left  unchecked,  this  mode  locking  gives  a  deadband  that  is  10,000  times 
that  required  for  accqitable  angular  rate  sensitivity. 

The  classical  ^iproach  to  any  deadband  problon  is  to  introduce  bias.  That  is,  to  artificially  push 
the  device  away  from  operation  near  zero  input  through  some  means.  For  the  RLG,  d^  bias  takes 
the  classical  form  of  mechanical  didier.  The  device  is  mechanically  shaken  at  high  frequency  and 
imperc^tible  amplitude,  so  that  the  true  rotational  rate  is  measured  as  an  increase  over  the  rate 
artificially  induced  by  the  didier.  The  practical  RLG  loses  its  theoretical  designation  as  a  solid  state 
device  b^use  of  this,  and  research  is  underway  to  discover  a  way  to  create  optical  bias. 


Unlike  earlier  modular  versions,  modem  RLGs  sacrifice  the  maintenance  omvenience  of 
modularity  for  the  structural  stability  of  a  single  monolithic  machined  block  for  all  duee  axes. 
Cavities  ate  precision  bored  to  hold  the  lasing  medium,  the  mirrors,  and  the  electrodes  into  a  single 
compact  package.  A  typical  configuration  is  a  cube  where  three  sets  of  four  faces  each  are  used  to 
implement  three  mutually  orthogonal  closed  loops.  One  such  loop  is  indicated  in  the  figure  below: 


Machined 

Block 


Figure  13  Typical  Monolithic  Ring  Laser  Gyro  Configuration 
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SJL  Ontiiai  Fibre  Gyroscopes 

Unlike  the  RLG,  the  fibre  optic  gyroscope  (FCXj).  employs  the  Sagnac  effect  directly.  It  achieves 
the  required  sensitivity  increase  by  increasing  the  path  length.  Recasting  the  earlier  Sagnac  length 
difference  result  in  terms  of  a  phase  difference  for  a  single  round  trip  between  the  counterrotating 


beams: 

, ,  27cAL 

A^  =  -^  = 

cX 

Now  by  using  n  coils  of  a  fibre  optic  light  pipe,  this  phase  difference  is  multiplied  by  n.  Substituting 
L  for  2ltr. 


4niiLrQ 

^♦  =  -7r-  = 

cX 

which  shows  that  the  phase  shift  grows  linearly  with  the  number  of  turns  of  the  fiber.  Current 
bendable  optical  fibers  allow  several  miles  of  extremely  thin  cable  to  be  coiled  into  a  compact 
package.  FOGs  may  be  simpler,  more  reliable,  and  less  costly  than  the  RLG.  They  do  not  suffer 
from  lock-in,  require  no  mechanical  dither,  and  require  no  precision  machined  blocks.  They  require 
no  critical  minors,  and  therefore  do  not  suffer  from  mirror  aging.  They  promise  to  replace  tl»  RLG 
even  as  it  continues  to  replace  the  mechanical  gyro  today. 

Practical  Optical  Fibre  Gvrosscopes 

The  feasibility  of  the  FCXj  had  to  await  the  development  of  the  optical  fiber,  and  its  development 
closely  follows  optical  communications  technology.  The  first  experimental  FCXjs  were  produced 
in  1976,  and  research  has  been  feverish  since.  The  first  all  single  mode  FOG  was  reported  in  1981. 
As  new  problems  are  discovered  and  quantified,  solutions  are  found  quickly  and  performance 
improves.  In  one  case,  the  Kerr  effect,  identified  in  1982,  two  solutions  to  die  problem  wm 
published  only  a  few  months  afterward. 

Chirrent  designs  ate  only  useful  for  modest  applications,  and  have  been  proposed  for  automobiles. 
These  employ  a  ting  interferometer  construct^  of  optical  fiber  coiled  into  a  compact  loop,  and  a 
semiconductor  laser,  similar  to  those  used  in  optical  disks.  A  beam  splitter  splits  the  beam  to  send 
it  in  both  directions  simultaneously  around  the  loop.  At  the  end,  both  beams  are  recombined  and 
sent  to  a  photodetector. 

An  ^ect  called  polarization  instability  is  now  believed  to  be  the  substantial  limiting  factor 
preventing  teal  performance  from  approaching  the  theoretical.  The  polarization  of  the  beam  in 
current  fib^  is  unstable,  so  interferometer  performance  is  degraded.  In  some  designs,  polarization 
maintaining  fibers  are  used  in  a  special  loop  before  the  main  loop,  to  stabilize  the  polarization  of 
the  beams.  With  this  enhancement,  the  research  community  now  believes  that  the  major  problons 
encountered  by  the  original  systems  are  undmstood,  and  the  FOG  is  now  entering  the  developmmit 
and  engineeting  phase.  Their  use  in  high  performance  initial  sensors  is  now  considered  to  be 
inevitable. 
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PART  III:lnertial  Navigation 

Systems'*^ 


Hie  name  inertial  navigation  or  inertial  guidance  derives  from  the  use  of  inertial  principles  by 
certain  types  of  navigation  systems.  These  systems  incorporate  a  large  number  of  components 
including  accelerometers,  gyroscopes,  compute,  and  sometimes  gimballed  stable  platforms  -  all 
of  which  are  engineered  to  perform  together  to  supply  navigation  information. 

Eariy  versions  of  the  INS  were  used  by  the  Pemiumiinde  group  in  Germany  in  World  War  n  to 
guicte  the  A3  series  rocket  and  later  the  V2  rocket  The  navigation  systems  were  enabled  during  the 
powered  ascent  but  were  shut  off  with  the  rocket  engines,  since  the  rest  of  the  vehicle  trajectory 
became  ballistic  and  therefore  predictable.  The  design  of  complete  systems  in  their  modem  form 
first  occurred  in  the  1940’s  and  is  credited  to  Qiarles  Draper  and  his  colleagues  at  MIT. 


1.  Basic  Principle  of  Inertial  Navigation 


UL 


of  the  INS 


These  devices  exhibit  a  large  number  of  ideal  characteristics.  For  the  military,  they  are  a  favorite 
navigation  aid  because  they  radiate  no  energy  that  may  assist  an  adversary  in  locating  the  vehicle. 
Also,  diey  requite  no  information  fitom  the  environment,  so  they  are  indifferent  to  weather 
conditimis,  and  the  attendant  poor  visibility  of  landmarks,  and  they  cannot  be  jammed  or  intmfmed 
with.  They  can  even  be  used  on  vehicles  travelling  at  high  enough  speeds  to  ionize  the  surrounding 
air,  and  effectively  jam  their  own  navigation  radar. 

They  can  be  used  easily  in  unsurveyed  areas  where  landmark  positions  are  unknown,  in  featureless 
terrain  where  there  are  no  landmarks,  and  in  situations  where  artificial  landmarks  cannot  be 
installed  for  any  of  a  variety  of  reasons.  They  are  ubiquitous  components  in  missiles,  tanks, 
satellites  and  long  haul  commercial  aircraft 


Initial  navigation  systems  OKS)  are  important  to  mobile  robotics  because  they  provide  the  most 
accurate  form  of  dead  reckoning  currently  available.  Also,  die  use  of  inertial  principles  implies  that 
these  devices  can  be  used  litmdly  anywhere  in  the  universe,  including  the  entire  surface  of  the 
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globe,  and  the  atmo^l^ie,  provided  the  local  gravitational  field  is  known.  They  ate  feasible 
alternatives  in  many  places  where  other  types  of  systems  cannot  be  used,  including  inteiplai^tary 
space,  underground,  and  undersea.  In  this  respect,  they  are  the  most  generally  applicable 
navigation  technology  in  existerrce. 

Some  people  now  envision  a  day  soon  when  inertial  navigation  systems  will  be  robust  enough  and 
cbe^  enough  to  become  standard  equipment  on  automobiles.  To  the  ancient  mariner,  these  little 
black  boxes  which  radiate  nothing  and  do  not  need  to  see,  touch,  or  otherwise  sense  the 
environment  would  have  likely  seemed  nothing  less  than  sorcery. 

The  ii»rtial  navigation  system  (INS)  can  be  distinguished  from  other  types  of  navigation 
instruments  by  the  following  characteristics: 

•  requiring  no  external  information,  are  jamproof 

•  radiating  nothing,  exhibit  perfect  stealth 

•  caimot  sense  accelerations  of  unpowered  space  flight 

•  can  operate  only  in  a  known  gravitational  field 

•  provide  continuous  or  sampled  position,  velocity,  and  acceleration 

•  achieve  highest  accuracies  after  warm  up  and  drift  trim 

•  most  errors  exhibit  Schuler  oscillation 

•  most  errors  are  time,  not  speed  or  distance  related 

•  requite  input  of  initial  position,  velocity,  and  attitude 

UL  ThftNMYftlMS 

The  basic  principle  employed  in  inertial  navigation  is  deduced  reckoning.  A  set  of  accelerometers 
ate  used  to  measure  acceleration  along  some  set  of  orthogonal  axes,  and  their  ouq>uts  are  integrated 
twice  to  determine  position.  The  integration  is  usually  p^otmed  numerically  in  a  computer. 
Althou^  the  above  statonent  of  the  principle  is  an  accurate  one,  it  does  not  address  a  large  number 
of  other  practical  engineering  considerations.  The  complexity  of  a  modem  inertial  navigation 
system,  incorporating  gimbals,  accelerometos,  gyros,  and  a  stable  table  results  directly  from  the 
rased  to  mediaiiixe  die  above  principle  in  teal  hardware. 

Radier  than  stating  the  genetic  design  of  an  INS,  the  discussion  will  start  by  explaining  why  all  of 
the  componmits  ate  necessary.  Most  of  the  issues  can  be  illustrated  by  examining  the  problems 
which  arise  wl^  accelerometers  ate  simply  str^rped  physically  to  a  vehicle  as  shown  in  Hgure 
14.  Although  it  is  costtnnary  to  mount  the  accelerometers  ^ong  orthogonal  axes,  it  is  not  necessary. 
Suppose  three  accelerometers  ate  fixed  such  that  their  sensitive  axes  are  oriented  with  one  along 
the  vehicle  longitudinal  axis,  one  transvase  to  the  vehicle  longitudinal  axis,  and  one  vertical. 
Furtha,  siqipose  that  the  accelerometer  ouqruts  ate  read  directly  and  integrated  twice. 

While  it  may  seem  that  this  is  a  basically  sound  approach,  die  reader  may  be  surprised  to  discover 
that  it  is  fimdamentally  flawed  in  severad  ways.  T^  naive  implementation  of  tte  basic  principle 
suffers  from  all  of  the  problems  listed  below. 
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Figure  14  Naive  INS  Concept 


1.2.1  Coordimite  Systems 

Deduced  reckoning  theoretically  consists  of  adding  up  a  number  of  differential  displacement 
vectors.  After  the  accelerometer  reading  is  integrated  over  a  time  step  twice,  the  result  is  a 
differential  displacement  of  the  vehicle  with  respect  to  the  earth.  In  this  case,  since  the  axes  of  the 
accelerometers  move  with  the  vehicle,  the  displacement  vectors  are  inherently  expressed  in  a 
coordinate  system  which  rotates  with  the  vehicle. 

In  order  to  add  vectors  of  any  kind,  they  must  all  first  be  expressed  in  a  common  coordinate 
sysiem^^  Further,  that  coordinate  system  must  be  fixed  with  respect  to  the  reference  firame  of 
interest  Th^fore,  all  displacement  vectors  must  first  be  expressed  in  a  coordinate  system  fixed 
with  respect  to  the  earth,  not  the  vehicle.  In  order  to  do  the  transformation,  a  sensor,  or  collection 
of  sensors,  is  needed  which  measures  the  total  rigid  body  orientation  of  the  vehicle  -  pitch,  yaw, 
and  rolL  The  solution  to  this  coordinate  system  problem  is  to  incorporate  a  set  gyroscopes  to 
measure  the  orientation  of  the  accelerometer  triad,  and  to  physically  or  computationally  resolve  the 
differmitial  displacements  into  an  earth  fixed  coordinate  system  before  adding  them. 

1.2.2  Referaioe  Frames 

Acceleration  is  inherently  an  inertially  lefermced  vector,  so  accelerometers  basically  measure 
acceleration  of  the  vehicle  with  respect  to  inertial  ^ace.  Also,  gyroscopes  inherently  measure 
inertial  rotation,  that  is,  rotation  with  respect  to  the  stars,  not  with  respect  to  the  earth.  These 
properties  result  firom  the  physical  principles  upon  which  the  devices  are  based.  The  solution  to  this 
reference  frame  problem  is  to  use  the  transformation  laws  of  mechanics  to  convert  the  inertiaUy 
ref^nced  quantities  to  earth  tefierenced  quantities.  In  the  case  of  the  accelerometers,  die  coriolis 
and  centrifugal  forces  caused  by  the  rotation  of  the  earth  must  be  accounted  for.  In  the  case  of  the 
gyros,  their  orimitation  refnence  must  be  somdiow  altered  so  that  it  rotates  widi  the  earth. 


45.  The  reader  is  eDComaged  to  review  the  discussioa  of  (be  distinction  between  coordmate  systems  and  ref- 
ereaoe  frames  in  an  earlier  section. 
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Gravitatioii 


Recall  that  accelerometers  are  adversely  affected  by  gravitational  fields,  because  they  actually 
measure  a  quantity  called  specific  force,  and  not  even  inertial  acceleration.  The  solution  to  this 
problem  is  to  account  for  the  efi^ect  of  the  gravitational  field  of  the  earth  in  the  mechanization 
equations. 

1^4Sttiiuiiary 

In  summary,  the  naive  implementation  of  an  inertial  navigation  system  suffers  from  the  following 
set  of  fundamental  problems: 

•  the  accelerometers  measure  die  wrong  quantity 

•  they  measure  it  in  the  wrong  reference  frame 

•  they  represent  it  in  e  wrong  coordinate  system 
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2  Generic  Mechanization 

It  is  the  engineering  solutions  to  the  problems  of  the  last  section  which  cause  the  modem  INS  to 
take  the  form  that  it  does.  The  following  discussion  will  explain  how  each  of  these  three  problems 
are  corrected  in  the  following  order: 

•  convert  specific  force  to  inertial  acceleration 

•  convert  inertial  motion  to  earth  referenced  motion 

•  express  earth  referenced  motion  in  the  navigation  coordinate  system. 

2JL  Specific  Force  to  Inertial  Acceleration 

Recall  that  the  accelerometer  functions  by  coupling  a  proof  mass  to  the  instrument  case  through  a 
calibrated  restraint  Hie  deflection  of  the  device  is  a  measure  of  the  force  or  torque  experienced  as 
a  result  of  its  accelerated  motion.  Higher  qualiQ'  devices  employ  rebalance  loops  which  prevent  the 
mass  from  moving,  and  the  effort  required  to  do  this  is  a  measure  of  specific  force. 

Such  a  device  can  be  modelled  conceptually  as  a  mass  attached  to  a  spring  which  is  constrained  to 
move  only  in  a  single  direction  within  a  frictionless  tube  as  shown  in  Figure  15.  Suppose  both  the 
tube  and  the  end  of  the  spring  are  attached  to  die  vehicle  and  the  mass  is  free  to  move  while  the 
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Let  the  vehicle  exert  a  force  T  on  the  spring  because  the  vehicle  is  moving,  and  let  the 
gravitational  force  be  W  and  the  inertial  acceleration  be  of  the  mass  m  .  The  device 
registers  the  specific  force,  ^  while  the  mass  responds  to  the  sum  of  these  two  forces,  the  spring 
responds  only  to  the  force  T  .  After  the  spring  has  finished  extending,  the  motion  of  the  mass  is 
the  same  as  that  of  the  vehicle  to  which  it  is  attached.  The  fundamental  equation  of  inertial 


navigation  is  Newton’s  second  law  applied  to  the  proof  mass; 

=  f +  W  =  mai 

Equation  1 

• 

The  specific  force  is  the  quantity  T/m.  The  inertial  acceleration  can  be  expressed  in  terms  of  the 
vitational  forces  as  follows: 


a,-  =  — 1-  —  =  t  +  w 

Equation  2 

1  m  m 

At  this  point,  it  is  possible  to  convert  the  specific  force  into  inertial  acceleration,  so  that  it  can  be 
integral.  Notice  that  explicit  knowledge  of  the  gravitational  field  strength  is  required  at  every 
position  of  the  vehicle.  Inertial  navigation  is  only  viable  when  this  field  is  known  or  known  to  be 
insignificant 

Since  most  inertial  systems  operate  near  the  earth,  the  gravitational  effects  of  the  other  heavenly 
bodies  are  relatively  small,  so  it  is  the  earth’s  gravitational  field  which  is  required.  A  practical 
gravitation  model  must  account  for  the  earth’s  oblate  spheroidal  shape  and  other  anomalies'*^,  but 
for  purposes  of  illustration,  consider  the  model  of  an  i^al  spherical,  homogeneous  earth  which  is 
given  by  Newton’s  law  of  universal  gravitation  applied  to  a  point  mass  at  the  earth’s  center. 
Let  G  be  the  universal  gravitation  constant,  M^be  the  mass  of  the  earth  and  be  the  radius 


GM, 

->  e  ^ 

W  - - i-r. 

Equation  3 

46.  Tbe  practice  of  iiuf>ping  the  earth’s  gravitation  is  called  geodesy. 
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2^ 

In  the  conversion  of  the  basic  equation  to  an  earth  referenced  form,  it  is  necessary  to  recognize  that 
all  inertial  frames  of  reference  are  equivalent  from  the  point  of  view  of  measuring  acceleration. 
Hence,  the  previous  result  can  be  considered  to  be  referenced  to  any  inertial  frame.  It  can  be  shown 
that  a  geocentric  frame  of  reference  that  is  not  rotating  with  respect  to  the  stars  is  inertial. 

In  the  following  result,  two  frames  of  reference  are  defined  to  be  situated  at  the  center  of  the  earth. 
The  geocentric  inertial  frame  is  associated  with  the  subscript  i  and  the  geocentric  frame  which 
rotates  at  the  earth’s  sidereal  rate  is  associated  with  the  subscript  e.  Another  frame  of  reference  is 
associated  with  the  moving  vehicle,  and  has  subscript  v. 

These  three  frames  are  in  r^tive  rotation.  Let  the  constant  rotation  of  the  earth  with  respect  to 
inertial  space  be  given  by  Q  ,  and  the  rotation  of  the  vehicle  frame  with  respect  to  the  earth  be 
given  by  p  .If  ^  denotes  the  inertial  rate  of  the  vehicle  frame,  it  can  be  expressed  as: 

®  =  d  +  p  Equation  4 

Let  f  ,  Vjj. ,  and  be  the  position  vector,  velocity,  and  acceleration  of  the  vehicle  measured  in 
the  frame  x.  Then,  since  the  position  vector  is  the  same  in  both  the  inertial  and  earth  frames,  the 
inertial 


The  result  of  differentiating  the  above  one  more  time  in  the  inertial  frame  is: 
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23  Vehicle  Rgfgrcnccd  Acceleratinn 

The  first  term  on  the  right  hand  side  can  be  referred  to  the  vehicle  frame  by  another  application  of 
the  theorem  of  Coriolis: 


dv 

dv 

Equation  7 

Substituting  Equation  7  and  Equation  5  into  Equation  6  yields: 


^^e  -  -  -  - 

aj  =  (^  )  +  (CO  +  Q)  xVg  +  Qx  (Qxfg) 


Equations 


Substituting  this  into  Equation  2: 


dv  _»  _» 

=  t- (©  +  Q)  X +  X  (Q xf 

til  Y 


Equation  9 


e-' 


This  is  one  of  the  most  convenient  forms  of  the  equation  of  motion  of  the  vehicle.  The  last  two 
terms  on  the  right  hand  side  are  both  functions  of  only  the  position  and  are  often  grouped  together 
and  called  gravity.  Gravity  is  part  gravitation  and  part  centrifugal  force: 


g  =  w-lix  (Qxr) 


Equation  10 


This  gives  the  vehicle  referenced  acceleration  as: 


In  precise  terms,  this  quantity  is  the  derivative  of  the  velocity  of  the  vehicle  with  respect  to  the 
earth,  where  the  derivative  is  computed  by  a  vehicle  fixed  observer.  It  is  not  the  derivative  of  the 
velocity  of  the  vehicle  with  respect  to  the  vehicle,  which  would  be  zero.  Further,  it  is  also  not  the 
derivative  of  the  vehicle  velocity  with  respect  to  the  earth  that  would  be  measured  by  an  earth  fixed 
observer,  because  it  includes  terms  for  the  rotation  of  the  vehicle  with  respect  to  the  earth. 


The  reason  for  the  use  of  such  a  quantity  is  that,  in  the  most  general  case,  the  accelerometers  are 
strapped  down  to  the  vehicle  chassis.  By  the  Theorem  of  Coriolis,  the  rotation  of  the  vehicle  with 
respect  to  the  earth  adds  an  extra  amount  to  the  the  rotation  of  the  moving  frame. 
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iaS.  vecior  int 


Solving 


or  the  quantities  of  interest: _ 

t 

Vg  =  j[t-  (c5  +  il)  xVg  +  g]dt  +  Vg^ 

^  Equation  12 


0 


Now  these  expressions  are  vector  integrals,  and  the  vectors  must  be  converted  into  a  common 
coordinate  system  before  they  are  integrated.  It  is  the  position  of  the  vehicle  with  respect  to  the 
earth  which  is  required,  so  the  integrands  must  be  referred  to  an  earth  fixed  (e  frame)  coordinate 
system  before  they  can  be  used. 


It  is  incorrect  to  integrate  them  in  a  coordinate  system  fixed  to  the  vehicle  frame.  One  way  to  see 
this  is  that  the  velocity  vector  is  often  of  constant  orientation  in  the  vehicle  frame  -  even  though  the 
vehicle  rotates  in  the  earth  frame.  Integration  of  the  forward  velocity  without  accounting  for 
rotation  leads  to  the  conclusion  that  motion  is  always  rectilinear,  which  obviously  cannot  be  right 

These  equations  can  be  solved  for  the  unknown  position,  and  velocity,  given  the  following 
information: 


•  a  model  of  the  earth’s  acceleration  due  to  gravitation  as  a  function  of  position 

•  the  earth’s  sidereal  rate  of  rotation 

•  the  vehicle  rotation  rate  with  respect  to  the  earth 

•  the  specific  forces  from  the  accelerometers 

•  the  initial  position 

•  the  initial  velocity 
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Z&.  Gmmc  MtchaniMtiftn 

The  function  of  the  INS  is  to  solve,  or  mechanize,  the  above  equations  in  order  to  detennine  the 
positional  variables  representing  the  motion  of  the  vehicle  over  the  surface  of  the  earth.  To  do  this, 
all  terms  in  the  vector  equations  must  first  be  expressed  in  a  common  coordinate  system.  This 
matter  will  be  deferred  until  the  next  section.  For  now,  consider  the  problem  in  its  coordinate 
system  independent  form. 

The  mechanization  equations  amount  to  a  need  to  do  four  things: 

•  add  gravitation  to  specific  force 

•  remove  centrifugal  force  due  to  the  vehicle  offset  from  the  center  of  the  earth 

•  remove  coriolis  force  due  to  the  vehicle’s  motion  on  the  surface  of  the  earth 

•  poform  two  integrations,  incorporating  initial  conditions 


The  process  of  solving  the  equations  in  a  computer  consists  of  initializing  the  positional  variables 
to  the  initial  conditions,  reading  the  accelerometers,  computing  gravitation,  coriolis  and  centrifugal 
terms  using  die  position  and  velocity  from  the  last  cycle,  transforming  specific  force  to  eardi 
referenced  acceleration,  and  finally  doing  the  integration  to  generate  position  and  velocity  for  the 
current  cycle.  The  solution  of  the  equations  can  be  visualized  with  tte  functional  block  diagram 
shown  below: 


Figure  16  Generic  Nfechanization 

Notice  that  the  system  provides  a  continuous  ouqiut  of  the  vehicle  position,  velocity,  and 
accderatitm  with  respect  to  the  cmiter  of  the  earth. 
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2al  Temporal  Error  Propagation 

It  is  instructive  to  whether  any  of  the  tenns  in  the  equations  are  insignificant  and  to  ask  what 
the  result  of  neglecting  it  would  be.  For  illustration  purpqs^,  let  the  vehicle  rotation  rate  match  the 
earth  rate  so  that  the  vehicle  is  always  level  and  c5  =  Oi .  The  earth’s  sidereal  rate  is  360°  in  24 
hours  or  7.27  X  10'^  rad/sec  and  the  radius  to  the  equator  is  7000  Km.  For  a  vehicle  at  the  equator, 
moving  eastward  at  a  velocity  of  10  meters  per  second,  and  accelerating  at  0.1  g,  the  following 
table  gives  the  magnitude  of  each  term; 


Term 

Nominal  Value 

Specific  Force  t 

0.1  g 

Gravitational  g 

1.0  g 

Centrifugal  Q  X  Q  X  f ^ 

3.8xlO-^g 

Coriolis  2Cl  X 

1.5xl0-^g 

Figure  17  Terms  in  tl»  Navigation  Equation 


While  it  may  seem  that  the  Coriolis  term  may  be  neglected,  rmnember  that  the  process  of 
integration  multiplies  acceleration  by  the  square  of  time,  so  that  after  1  hour,  neglecting  the 
Coriolis  term  accounts  for  almost  10  kilometers  of  accumulated  error.  Note  that  this  error  occurs 
even  when  the  system  is  stationary  for  that  time  period  -  so  it  is  truly  time  dependent 
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3  Complete  Solution  in  Geocentric  Cartesian  Coordinates 

One  of  the  reasons  for  the  paiticular  expression  of  the  navigation  equations  given  in  a  earlier 
section  is  that  the  vector  ^  is  very  close  to  tl»  direct  output  of  the  gyroscopes.  In  the  strapdown 
case,  for  an  orthogonal  arrangement  of  gyros,  each  has  its  sensitive  axis  aligned  with  one  of  the 
designated  axes  of  the  vehicle,  and  each  responds  directly  to  the  component  of  the  inertial  angular 
velocity  which  is  aligned  with  it 


ITTrTilTIl 


Let  a  cartesian  coordinate  system  called  the  navigation  coordinate  system,  indicated  by  the  letter 
n,  be  assigned  to  the  center  of  the  earth  in  accordance  with  the  conventions  set  out  in  [27].  Let  a 
set  of  axes  called  the  body  coordinate  system  be  assigned  to  the  vehicle,  and  indicated  by  the  letter 
b.  In  accordance  with  the  same  conventions,  z  points  up,  y  forward,  and  x  out  the  right  side.  In  the 
body  coordinate  system,  the  inertial  angular  velociQ^  referred  to  body  coordinates  can  be 
expressed  as: 


T 

Equation  13 

The  rotation  matrix  which  converts  coordinates  from  the  body  frame  to  the  navigation  frame  will 
be  called  the  direction  cosine  matrix  and  denoted  by  R^.  It  can  be  used  to  convert  coordinates  to 
the  earth  frame  as  follows: 


co“  =  R?® 


Equation  14 


This  quantity  represents  the  angular  velocity  of  the  vehicle  widi  respect  to  inertial  space,  referred 
to  earth  coordinates.  The  earth’s  sidereal  rotational  rate  referred  to  earth  coordinates  is  given  by: 


£f  =  [OOQ] 


Equatiim  15 


These  two  terms  give  die  dual  angular  velocity  term  of  the  navigation  equations. 


The  gyros  indicate  the  angular  velociQr  of  the  vehicle  with  resp^t  to  initial  space.  In  order  to 
determine  the  attitude  with  respect  to  the  earth,  the  eardi  rate  Q  must  be  removed.  Recall  the 
definitions  of  the  three  angular  velocities  and  then: 


^  =  ©“-if  =  R;®*’-a 


Equatkm  16 
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The  accelerometer  specific  forces  are  intrinsically  referred  to  the  body  frame  thus: 


These  are  referred  to  eaitii  coordinates  in  the  usual  way: 


t“  =  Rgt* 


Equatum  18 


M.  ikaxitx 

Gravity  is  intrinsically  referred  to  earth  coordinates: 


^  Vetodtv  Solution 

Finally,  the  initial  velocity  must  be  leferted  to  earth  coordinates,  and  all  subsequent  computed 


t 

v“  =  /[Rbt*’-  (RS5^  +  ^)  xv“  +  g“]dt  +  vj 

EqiialicHDilO 

0 

Attitnde 

From  [27],  the  rates  of  the  Euler  angles  describing  the  orientation  of  the  vehicle  ate  nonlineatly 
related  to  the  angular  velocity  vector  p°  by^^: 


T  = 


e 

♦ 


c<)>  0  s<|> 
tOstj)  1  -t6c<|> 

ce^  c0 


This  mtpiession  can  be  directly  integrated  in  order  to  determine  the  Euler  angle  expression  of  the 
attitude  of  the  vdiicle  with  respect  to  the  earth.  Thus^: 


t 

r  =  f[r(ffl“-Q“)]dt+TQ 

Eqiiatk>n22 

0 

47.1b  [31],  die  matrix  F  is  called  G. 

48.  Recall  earikr  commenis  on  the  distmgaished  nature  of  3D  rotadoas.  niis  qnaotity  is  not  a  vector  in  iilie 
seme  of  a  vector  q»oe  of  linear  algAra. 
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Once  the  Euler  angles  are  known,  the  direction 

(cyc^-sys6s^)  -sycO  (cys^ 
Rfc  *  (syc^  +  cys6s^)  cyc0  (sys^ 


(cyc^-sys6s^)  -sycO  (cys^+sys6c(^) 
(syc^  +  cys6s^)  eyed  (sys^eysOe^) 
-e6s^  s0  cOe^ 


matrix  is  given  in  [27]  as: 


Equation  23 


Finally,  the  position  solution  is  dven  b 


t 

f“  =  fv“dt  +  fQ 

Equation  24 

J 

0 
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4  Schuler  Dynamics 


Tbe  gyrocompass  oscillates  like  a  pendulum  with  a  characteristic  period  that  is  caused  by  the 
rotation  of  the  earth.  This  principle  of  pendulous  oscillation  is  central  to  the  understanding  of 
inertial  systems  and  is  (teveloped  further  here. 

A  gyrocompass  mounted  on  a  moving  ship  is  subject  to  motions  much  more  pronounced  than  the 
earth’s  rotation.  In  an  important  pi^)er  published  in  1923,  Max  Schuler,  who  was  investigating  this 
problem  of  the  gyrocompass,  pointed  out  that  if  the  gyrocompass  period  was  about  84  minutes, 
then  it  will  be  unaffected  by  motions  of  any  vehicle. 

Ai,  Sdittter  Ptirifld 

Many  ^stmn  variable  errors  in  an  inertial  navigation  system  exhibit  a  period  called  the  Schuler 
)eriod^  in  the  two  level  directions.  The  Schuler  period  is  given  by  the  following  relation  where 
^  is  the  radius  of  the  earth  and  g  is  gravity  : 

Re 

T  =  2jc(— )  «84.4miii 

g 

The  correct  application  of  the  Shuler  principle  is  an  essratial  component  of  precision  inertial 
navigation  systems.  The  Schuler  oscillation  is  the  natural  horizontal  oscillation  in  the  error  in 
computing  gravitation.  , 

SchnlerTVmiiiy 

While  precision  components  is  one  of  die  reasons  why  inertial  navigation  is  possible,  the  primary 
reason  for  the  feasibility  of  terrestrial  inertial  navigation  is  that  all  systems  which  navigate  close  to 
the  earth  experience  horizontal  oscillatoty  errors  reflecting  the  famous  Schuler  period  of  84 
minutes.  Horizontal  error,  being  oscillatory,  is  inherently  bounded. 


Sdiuier  tnirfng  can  be  described  as  die  construction  of  a  device  whose  natural  frequency  is  the 
same,  as  that  of  a  pendulum  whose  lengdi  is  the  radius  of  the  earth.  A  pendulum  attached  to  a 
vehicle  moving  on  the  surface  of  the  earth  wOl  be  deflected  because  of  the  vehicle  acceleration  as 
shown  below.  In  fact,  this  is  the  principle  of  the  pendulous  accelerometer. 


Figure  18  Accelerated  Pendulum 


49.  Whereas  a  geosyndmooiis  ofbft  has  a  period  of  exactly  24  boms,  low  eaiUi  ortnts  iqieat  sevoal  times  a 
diy.  It  tunis  cut  that  the  Stanler  period  is  also  equal  to  tbe  period  of  an  obit  at  treetopbeigbL  , 
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The  Schuler  pendulum,  having  a  radius  equal  to  Uk  earth’s  radius,  can  be  accelerated  arbitrarily 
around  the  surface  of  a  sphere  without  being  dismrbed  by  the  motion  of  the  vehicle.  In  practice, 
one  cannot  construct  such  a  pendulum,  but  any  dynamic  system  with  the  correct  natural  frequency 
will  exhibit  this  same  stability  proper^.  In  the  context  of  the  INS,  it  is  the  natural  frequency  which 
is  important 

dLl  SdmlCTLwp 

A  stable  table  for  a  locally  level  system  can  be  constructed  on  the  simple  principle  that  when  the 
table  is  not  level,  the  horizontal  accelerometers  will  measure  a  component  of  gravity.  Consider  a 
servo,  called  a  Sdiuler  loop,  which  computes  the  angular  velocity  of  the  vehicle  with  respect  to 
the  earth  and  rotates  the  table  ^propriately.  Nominally,  the  required  rate  of  rotation  of  the  table  is 
equal  to  the  earth’s  rate  with  respect  to  inertial  space  plus  the  vehicle  rate  over  the  earth’s  surface. 
The  accelerometer  ouq>ut  is  integrated  to  give  a  velocity  signal,  and  this  signal  is  then  divided  by 
tt»  radius  to  the  earth’s  center  in  order  to  get  the  angular  velocity.  This  signal  is  used  to  torque  the 
platform  to  keep  it  level.  The  velocity  signal  is  also  integrated  a  second  time  to  give  position.  These 
ideas  ate  summarized  in  the  functional  block  diagram  below. 


Figure  19  Schuler  Loop 


Ihe  feedback  path  is  ctmstructed  by  recognizing  that  the  accelerometer  specific  force  includes  any 
component  of  gravity  due  to  an  off-level  condition.  This  feedback  is  inherent  in  the  accelerometer 
-  dtere  is  no  electrical  coimection  required.  If  the  table  is  off  level,  the  accelerometer  wiU  measure 
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There  is  an  inherent  minus  sign  in  the  process  as  illustrated  below  using  a  simple  spring  model  of 
an  accelerometer: 


;  e 

1)  platform  2)  spring  /  3)  which  4)  which  cause 

tilts  off  deflects  is  apparent  the  platform  tc 

level  this  this  way  th« 

way  ^ 


o 


motion 
this  way 


rotate  this  waj 
to  stay  level 


o 


Figure  20  Accelerometer  Negative  Feedback 

The  dynamics  of  this  loop  can  be  investigated  easily  through  the  Laplace  transform  techniques  of 
linear  systems  analysis.  However,  for  the  benefit  of  readers  unfamiliar  with  such  techniques,  the 
time  domain  differential  equations  of  the  systems  will  be  used.  The  accel  'lometer  ouqiut  can  be 
written  in  terms  of  itself  by  following  the  feedback  loop  around  backwards  as  follows: 


a-g0  =  a-gja)dt  =  a-~Jvdt 

Equation  26 

a-g0  =  a-ljj[a-g0]dt 

Next  the  common  term  of  a  is  canceled,  and  the  result  is  differentiated  twice  to  obtain: 


g8  =  |(a-g0) 


Equation  27 


This  is  the  differential  equation  describing  the  motion  of  the  stable  table.  It  is  the  characteristic 
equation  of  forced  simple  harmonic  motion,  which  implies  that  the  table  will  naturally  oscillate  at 
a  period  of: _ 


T  =  27CJ—  =  84minutes 


Equation  28 
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^  Faucanlt  Effect 

The  Schuler  period  can  be  viewed  as  the  oscillation  of  a  pendulum  whose  length  is  the  radius  of 
the  earth.  It  is  the  period  of  the  pendulum  itself.  This  conceptual  earth-radius  pendulum  is  also 
subject  to  an  oscillation  of  the  plane  of  its  rotation  just  as  is  the  Faucault  pendulum.  This  period  is 
24  hours  at  the  equator  and  increases  to  infinity  at  the  poles.  For  a  period  T,  earth  rotation  rate  of 
Q  and  latitude  X,  the  equation  is: 

T  - 

~  Qsin'K 
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5  Inertial  Error  Dynamics 

A  more  refined  analysis  of  the  phenomenon  of  gravity  feedback  leads  to  the  conclusion  that  most 
errors  in  an  INS  oscillate  with  the  Schuler  period.  Such  an  analysis  would  consider  the  particular 
mechanization  of  the  system  -  since  many  errors  are  most  easily  described  in  the  mechartized 
coordinate  system,  and  would  consider  the  effects  of  all  significant  errors  simultaneously. 
Unfortunately,  such  an  analysis  consumes  more  space  than  is  available  here.  However,  it  is 
possible  to  illuminate  the  basic  ideas  by  considering  a  single  error  source  -  accelerometer  bias,  and 
developing  the  equations  describing  its  propagation  throughout  the  system  in  a  mechaiuzation 
independent  manner. 


£1  Perturbative  Analysis 

The  performance  of  the  Schuler  tuned  system  in  the  presence  of  sensor  norudealities  can  be 
investigated  by  a  technique  known  as  perturbation  analysis.  To  do  this,  consider  the  basic 
navigation  equation  -  inertial  acceleration  expressed  in  terms  of  the  specific  force  indicated  by  the 
accelerometers  a  ,  and  gravitation  g  : 

^  ^  Equation  29 

^inertial  ~  ,  2  ”  ^  8 

at 


In  perturbation  analysis,  a  perturbative  error  is  tq)plied  to  the  sensed  specific  force  and  the  effect 
of  this  on  the  system  output  is  investigated.  Let  the  indicated  specific  force  include  an  error  denoted 
Aa ,  and  cause  errors  in  the  computed  position  and  gravitation  denoted  by  Af  and  Ag .  This  is 


Where  the  subscripts  i  and  t  represent  indicated  and  true  quantities.  Substituting  this  back  into  the 


This  is  the  differential  equation  which  describes  the  propagation  of  errors  from  the  accelerometer 
to  the  position  and  gravity  computations. 
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1£1 


The  gravitational  force  is  a  function  only  of  position.  A  Taylor  series  expansion  assuming  a 
spherical  homogeneous  earth  aUows  rewriting  the  gravitation  error  in  terms  of  the  error  in  the 
position  as  follows: _ 


GM.  _ 

GM  . 

Equation  32 

6  - 

-2  — 

..3/2 

(rr) 

Ag  = 

r^gi 

LafJ 

II 

< 

GM  ..  .GM, 
- ^Ar  +  3— 

r  r 

r '  Ar)  r  by  the  product  rule 
of  differentiation 

Substituting  this  into  Equation  31  yields: 


-^At  =  Aa - x-Ar  +  3— r-  (r  •  Ar)  r 

dt^ 


Equation  33 


Solution  Near  the  Earth^s  Surface 

Further  analysis  of  this  equation  requires  that  a  coordinate  system  be  adopted.  Let  the  origin  be 
placed  at  the  center  of  the  earth,  and  the  three  cartesian  axes  be  oriented  s^itraiily.  Writing  this 
out  in  component  form: _ 


...  GM .  -GM  .  .  A  .  A  \  A 
Ax  +  — :r-Ax  -  3— r-  (xAx  +  yAy  +  zAz)  x  =  Aa 

AyH — ^Ay-3 — ^  (xAx  +  yAy  +  zAz)  y  =  Aa 
r  r  ^ 

Az  +  — ^Az -  3—^  (xAx  +  yAy  +  zAz)  z  =  Aa^ 
r 


Equation  341 


Any  particular  solution  to  these  equations  requires  knowledge  of  the  trajectory  followed  by  the 
vehicle.  Let  the  start  point  for  the  system  be  along  the  z  axis  on  the  surface  of  the  earth,  and  let  the 
vehicle  trajectory  remain  close  to  this  point  so  that  x  =  y  =  0  and  z  =  r  =  R  . 
Under  this  assumption,  the  cross  coupling  terms  in  the  equations  cancel  and  they  reduce  to: 

Ax+  (GM/R^)Ax  =  Aa^ 

Ay  +  (GM/R^)  Ay  =  Aa^  Equation  35 

Az  -  (2GM/R^)Az  =  Aa^ 
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Let  gQ  and  Rq  be  the  gravitational  acceleration  of  and  radius  to  the  local  region,  which  have 
been  assumed  to  be  approximately  constant.  If  the  accelerometer  errors  are  assumed  to  be  constant 


SA.  Error  Dynamics  of  the  Free  Inertial  System 

Hence,  the  horizontal  position  channels  oscillate  at  the  Schuler  frequency  with  an  amplitude 
proportional  to  the  accelerometer  bias.  The  analysis  is  mechanization  independent,  so  this 
oscillation  arises  in  all  systems  which  operate  near  the  earth.  The  main  conclusion  is  that  while  it 
is  unfortunate  that  gravity  needs  to  be  modelled  to  use  an  INS,  the  very  existence  of  a  gravitational 
field  inherently  stabilizes  the  process  of  dead  reckoning  from  an  accelerometer.  Although  the 
analysis  has  assumed  only  accelerometer  errors,  many  other  sensor  nonidealities  cause  similar 
behavior. 

The  vertical  channel  is  divergent  with  a  scale  factor  proportional  to  the  accelerometer  bias  and  a 
time  constant  related  to  the  Schuler  period.  Many  local  vertical  systems  possess  no  vertical 
computation  channel  at  all  for  this  reason.  If  they  do,  the  chaimel  must  be  augmented  by  other 
measurements  to  stabilize  it.  The  next  section  discusses  the  available  methods  of  stabilization. 

Notice  that  for  periods  of  time  which  are  small  compared  with  the  84  minute  Schuler  period,  the 
Taylor  approximations  for  the  cosine  and  hyperbolic  cosine  in  Equation  36  both  indicate  a 
parabolic  growth  of  the  accelerometer  bias  with  time  as  would  be  expected  without  gravity 
damping.  This  implies  that  while  the  error  may  be  bounded  in  the  horizontal  case,  the  bound  may 
be  a  large  number  depending  on  the  size  of  the  sensor  error.  Hence,  gravitational  feedback  can  be 
viewed  as  having  the  effect  of  converting  three  parabolic  error  profiles  into  two  oscillatory  ones 
and  one  exponential  one. 

These  results  are  independent  of  mechanization  since  no  particular  mechanization  was  assumed. 
The  fundamental  message  of  the  equation  is  that  errors  in  accelerometer  indication  give  rise  to 
errors  in  computed  position.  These  errors  in  position  give  rise  to  errors  in  computed  gravitation. 
The  errors  in  gravitation  tend  to  be  opposite  in  direction  to  the  errors  in  acceleration  which  caused 
them.  The  errors  may  arise  from  physical  deflections  or  considerations  of  a  more  computational 
nature.  That  is,  errors  oscillate  in  Schuler  tuned  systems  for  exactly  the  same  reason  that  a 
pendulum  oscillates  in  a  gravitational  field  -  the  mathematics  are  identical. 
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6  Implementation 

This  section  witi.  briefly  discuss  some  of  the  important  aspects  of  complete  inertial  navigation 
systems.  Today’s  commercial  systems  can  achieve  accuracies  on  the  order  of  0.2  nautical  miles  of 
error  per  hour  of  operation. 

All  inertial  navigation  systems  use  gyroscopes  to  measure  vehicle  rotation  in  inertial  space,  and 
accelerometers  to  measure  inertial  specific  force.  Two  or  three  accelerometers  are  used  to  measure 
specific  force  along  as  many  axes  and  three  single  degree  of  freedom  or  a  pair  of  two  degree  of 
freedom  gyroscopes  complete  the  sensor  package.  This  package  is  often  called  the  inertial 
measurement  unit  (IMU)  or  inertial  reference  unit  (IRU). 

Computers  perform  a  large  number  of  computations  in  order  to  convert  the  IMU  outputs  into  the 
required  earth-relative  motion  of  the  vehicle.  In  modem  systems,  digital  computers  are  used  and 
the  integrations  are  performed  using  digital  quadrature  ^gorithms.  In  many  cases,  integrating 
accelerometers  and/or  integrating  gyroscopes  are  used  which  perform  one  of  the  integrations 
through  the  addition  of  a  damping  term  in  their  transfer  functions. 

Although  simple  computational  models  of  the  gravitation  of  the  oblate  spheroidal  earth  are 
available,  many  systems  do  not  employ  them  at  present  Modem  systems  employ  sophisticated 
high  speed  integration  algorithms  and  real-time  Ktdman  filtering  to  improve  performance. 

Real  INS’s  ate  typically  designed  to  report  position  in  a  particular  frame  of  reference  and 
coordinate  system.  For  terrestrial  applications,  an  earth  fixed  frame  is  called  for,  while  for  extra¬ 
terrestrial  applications,  a  space  reference  is  more  useful  Devices  such  as  missiles,  even  though 
they  operate  substantially  in  space,  must  be  aware  that  their  target  will  have  moved  in  inertial  space 
since  the  time  they  were  launched.  All  devices  which  navigate  close  to  the  earth  experience  errors 
reflecting  the  famous  Schuler  period,  and  vertical  error  increases  without  bound  over  time  unless 
stabilized  by  periodic  external  updates 

ILL  Gimballed  and  Strapdown  Mechanizations 

The  implementation  of  a  particular  navigation  coordinate  system  involves  either  actively 
controlling  or  passively  keeping  track  of  the  orientation  of  the  accelerometers  as  the  vehicle  moves. 
A  trade-off  exists  between  mechartical  complexity  and  computational  complexity  which,  over 
time,  has  tended  to  drive  system  designs  toward  computational  complexity  as  computer  systems 
have  improved  in  performance.  Three  classes  of  systems  are  distinguished  based  upon  their  relative 
position  on  this  spectrum. 

6.1.1  Gimballed  Systems 

The  first  class,  known  as  gimballed  or  geometric  systems,  employ  a  stabilized  platform  which  is 
actively  servoed  to  the  required  orientation.  This  type  of  system  was  the  first  practical  class  to  be 
developed  historically  since  navigational  information  was  available  directly  from  the  gimbal 
angles.  Up  to  five  gimbals  ate  used  to  instrument  both  an  inertially  nonrotating  reference  firame  and 
a  local  navigation  frame.  These  systems  require  only  minimal  computational  capacity.  In  practice, 
they  have  been  replaced  by  another  class  of  systems,  called  semi-analytic. 
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6.1^  Semi* Analytic  Systems 


The  second  class  of  systems  ate  called  semi-analytic.  These  ate  gimballed  systems  which  control 
the  orientation  of  a  platform  in  order  to  instrument  the  navigation  frame  only.  Computations 
transform  the  sensor  outputs  from  the  inertial  frame  and  those  systems  which  do  not  actively 
maintain  a  north  reference  must  also  transform  coordinates  from  the  platform  axes  to  local  north 
and  east  directions.  Computations  are  also  used  to  determine  the  latitude  and  longitude  of  the 
vehicle. 

Gimballed  systems  employ  gyro  torquers  which  drive  the  platform  to  rotate  as  required.  Untorqued 
gyros  naturally  implement  an  inertial  reference.  The  total  torquer  signal  may  include  compensation 
for  the  gyro  drift  rate,  the  earth’s  sidereal  rate  and  angular  velocity  of  the  vehicle  with  respect  to 
the  earth,  called  the  vehicle  rate. 

6.1,3  Strapdown  Systems 

The  third  class,  known  as  strapdown  or  analytic  systems,  are  said  to  be  computationally 
stabilized.  Strapdown  systems  get  their  name  from  the  fact  that  they  strap  the  sensors  directly  to 
the  vehicle  chassis,  eliminating  the  platform  gimbals,  at  the  cost  of  causing  the  components  to 
suffer  the  dynamic  rotation  of  the  vehicle.  Strapdown  mechanization  places  more  stringent 
requirements  on  the  sensory  components  used,  and  on  the  computational  throughput,  so  titese 
systems  have  appeared  relatively  recently  in  the  history  of  development  of  inertial  navigation.  Such 
systems  promise  to  replace  semi-analytic  systems  in  most  applications. 

These  systems  are  analysed  mathematically  in  an  identical  manner  to  gimballed  systems  with  the 
exception  that  the  accelerometer  readouts  are  first  converted  to  a  common  coordinate  system  using 
a  matrix  transformation,  often  called  the  direction  cosine  matrix.  There  are  many  advantages  to 
the  strapdown  approach,  other  than  the  obvious  elimination  of  the  mass,  volume,  and  mechanical 
complexity  of  a  gimballed  platform. 

These  systems  use  rate  gyros  which  can  be  as  small  as  one  cubic  inch  in  volume  admit 
precessiotud  rates  up  to  300  times  that  of  gimballed  systems.  High  processional  rates  are  necessary 
to  faithfully  track  the  motion  of  the  vehicle  to  which  they  ate  attached.  Torquers  ate  required  to  be 
much  more  powerful.  Accelerometers  can  be  as  small  as  5  grams  mass  and  0.05  cu  inch  in  volume. 
The  operating  principle  is  to  use  very  high  speed  rebalance  loops  which  continually  torque  the 
gyros  back  to  their  null  vehicle  aligned  positions.  The  effort  required  to  do  this  is  proportional  to 
the  angular  rate  of  the  vehicle.  The  rebalance  loops  report  current  level  of  effort  to  the 
computational  stabilization  algorithm  which  maintains  the  current  vehicle  orientation  in  inertial 
space. 

The  test  of  the  strapdown  system  is  identical  to  a  gimballed  system.  Strapdown  systems  are  aligned 
and  trimmed  in  a  maimer  analogous  to  gimballed  systems.  Strapdown  systems  are  becoming 
particularly  suited  to  the  relatively  low  accuracy  requirements  of  land  navigating  mobile  robots. 


I 
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Perfa^>s  the  most  important  distinguishing  characteristic  of  inertial  navigation  systems  is  the 
coordinate  system  in  which  the  navigation  equations  are  mechanized.  Hie  two  categories  here  are 
space  stabilized,  and  earth  stabilized. 

<1.2.1  Space  Stabili^  Systems 

Space  stabilized  systems  employ  an  inertially  fixed  coordinate  system.  As  a  result,  the  navigation 
equations  take  a  particularly  simple  form.  Such  systems  are  very  inconvenient  for  use  in  terrestrial 
iqiplications.  Space  applications  often  use  hybrid  stellar-inertial  navigation  systems  employing 
star-trackers  since  star  sighting  is  both  easy  and  a  very  accurate  method  of  determining 
orientation. 

6.2J2  Eardi  Stabilized  Systems 

Earth  stidiilized  systems  maintain  orientation  referenced  in  various  ways  to  the  earth.  Eardi  fixed 
systems  employ  a  coordinate  system  fixed  with  respect  to  the  rotating  earth.  These  are  also  known 
as  base  p<rat  or  launch  point  systems  since  the  origin  of  coordinates  does  not  move  with  die 
vehicle  but  rather  remains  fixed  at  the  point  of  initialization. 

Local  vertical  systems  are  the  most  common  earth  stabilized  systems  and  employ  coordinates 
referenced  to  the  local  vertical  direction  with  the  other  two  axes  horizontal.  In  case,  because 
the  coordinates  rotate  and  move  with  the  vdiicle  they  are  called  locally  levd. 

Both  kinds  of  earth  stabilized  systems  require  that  an  earth’s  rate  signal  be  generated  by  the  system 
to  drive  the  platform  or  direction  cosine  transform  to  rotate  with  the  earth  and  can  ignore  the  gravity 
vector  provided  the  platform  is  sufficiently  stable.  For  local  vertical  systems,  the  compensation 
signal  must  also  account  for  vehicle  motion,  for,  depending  on  latitude,  each  nautical  mile  of  travel 
corresponds  to  up  to  one  minute  of  arc  of  change  in  the  local  vertical. 

Among  local  v^cal  systems,  there  are  three  common  variants.  Norfii-slaved  systems  maintain 
one  coordinate  axis  pointed  north  at  all  times.  Axes  are  typically  called  north,  east,  and  aamuth^^. 
Free  azimadi  and  wander  azimiith  systems  allow  the  level  coordinate  axes  to  rotate  with  respect 
to  north  about  the  vertical  axis  as  a  function  of  latitude  and  computations  are  employed  to  keep 
track  of  the  angle  of  rotation. 

A  summary  of  this  taxonomy  of  inertial  systems  is  provided  in  Figure  21: 


SO.  Azimodi  refers  to  the  angle  made  by  an  object  with  respea  to  an  aibitrary  r^eroice.  The  angle  is  mea¬ 
sured  in  the  local  earth  tangent  plane,  novation  is  the  angle  an  object  makes  widi  the  horizon. 
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Note:  any  system  above  can  be  analytic,  semi-analytic, 
or  geometric. 


Figure  21  Taxonomy  of  Inertial  Systems 
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The  three  most  common  terrestrial  INS  mechanizations  differ  only  in  the  way  that  the  platform^^ 
rotates  with  respect  to  the  earth.  For  strapdown  systems,  the  same  options  are  available  and  the 
direction  cosine  computations  are  driven  by  ^propriate  rates  to  implement  a  pseudo  stable 
platform.  In  the  notation  of  earlier  sections,  they  differ  in  the  prescription  of  the  vector  ^  .  These 
systems  will  be  summarized  below. 

6.3.1  Nortfi  Slaved,  Locally  Level  System 

This  system  is  defined  in  the  following  figure.  It  has  the  advantage  shared  by  aU  locally  level 
systems  that  compensation  of  the  horizontal  channels  for  the  gravity  vector  is  not  necessary,  since 
it  always  points  along  the  z  axis.  like  all  earth  st2d)ilized  systems,  gyro  torquers  are  required. 

It  permits  easy  calculation  of  latitude  and  longimde  and  is  commonly  used  in  long  term  cruise 
systems.  It  indicates  the  vertical  so  it  can  be  used  easily  to  assist  military  weapons  systems  and  it 
can  be  easily  connected  with  many  types  of  navigation  aids  which  can  be  used  to  reduce  errors  in 
the  inertial  system  without  transforming  coordinates. 

One  serious  drawback  of  this  system  is  that  a  singularity  exists  at  the  poles,  because  the  x  axis  must 
rotate  very  rapidly  around  the  azimuth  near  either  of  the  poles.  In  gimballed  systems,  this  marufests 
itself  as  high  pl^orm  azimuth  rates.  In  strapdown  systems,  the  high  rotation  rates  limit  the 
accuracy  of  the  computations.  This  is  one  of  two  aspects  to  the  pole  problon  for  the  INS.  The 


SI.  Hie  i^atfonn  fiame  in  gimballed  oonfiguratkms  is  equivalent  to  ibe  vdiicle  frame  in  strsqxlown  systems. 
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63^  Free  Azimuth,  Locally  Levd  System 


This  system  is  defined  in  the  following  figure.  It  solves  the  high  polar  platform  rate  problem  by 
provi(&g  no  inertial  platform  rate  about  the  vertical  axis.  The  torquing  error  of  the  z  gyro  is 
ordinarily  worse  than  the  two  other  axes,  so  it  is  advantageous  to  remove  it  in  this  way.  The 
azimuth  is  allowed  to  wander  with  respect  to  north  and  computations  keep  track  of  the  true 
direction  of  north.  £a  a  sense,  such  systems  are  hybrid  gimballed  -  strapdown  systems.  The  angle 
between  the  x  axis  and  true  north  is  called  the  wander  angle.  Wander  azimuth  is  a  variant  where 


the  z  component  of  the  inertial  platform  rate  is  equal  to  the  z  component  of  the  inertial  earth  rate. 


dA  ,  dX 
(0.^  =  (Cl  +  j-  )  cosKcosa  -  -j-  sina 
*  dt  at 

^  dA  ,  dX 

<o„  =  (Cl  +  :n  )  cosXsina—j2  cosa 
y  '  dt  '  dt  j 

<ss,  =  0  / 

For  Wairder  Azimuth:  / 

25^  =  —ClsinX  I 


center  at  present  vehicle  location  ^ 
X  axis  horizontal,  at  angle  a  with  north 
y  axis  horizontal,  at  angle  a  with  east 
z  axis  vertical,  down 


Figure  23  Free  Azimuth,  Locally  Level  System 
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Di^^t  Plane,  Locally  Level  System 

This  system  is  defined  in  the  following  tiguie.  In  this  case,  while  the  platfonn  moves  with  the 
vehicle,  its  orientation  with  respect  to  the  earth  remains  fixed  at  its  initial  value.  It  takes  advantage 
of  the  fact  that  constant  platform  rates  can  be  generated  much  more  accurately  than  vari^le  ones. 
These  rates  can  be  precisely  determined  at  start-up.  Such  systems  are  usually  employed  when 
computation  of  latitude  and  longitude  is  not  required  and  when  the  mission  excursion  is  limited. 
For  excursions  beyond  a  few  miles  from  the  base  point,  compensation  for  gravitation  must  be 
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While  the  stabilized  coordinate  system  is  the  system  in  which  the  vector  equations  of  the  first 
section  ate  rejMesented  and  manipulated  internally  by  the  computer,  the  system  ouq)ut  is  often 
converted  to  other  coordinates  for  reasons  of  convenience.  Two  commonly  used  systems  are 
latHude  and  longitude  and  the  military  grid  or  UTM  system. 

6.4.1  Latitude  and  Longitude 

The  generic  INS  mechanization  equations  of  the  previous  section  provided  the  earth  referenced 
IK)sition  vector  of  the  vehicle  as  output  For  local  vertical  systems,  it  is  more  common,  in  practice, 
to  compute  latitude  and  longitude  by  integrating  the  velocity,  so  that  the  position  vector  is  never 
explicitly  computed.  For  example,  for  local  vertical  systems,  this  is  particularly  easy  to  do. 

Let  A  be  the  vehicle  longitude  and  X  its  latitude.  The  latitude,  longitude,  and  altitude  of  the 
vehicle  relate  to  the  vehicle  velociQr  by: 


dk  _ 

dt  ”  ^north  ‘  ^ 

dA 


dt  ^down 


Equation  37 


Latitude  and  longitude  are  angular  measurements,  in  fact,  spherical  polar  coordinates.  Several 
points  about  these  coordinates  must  be  remembered^^.  Lines  of  constant  latitude,  called  parallds, 
are  normal  to  the  direction  in  which  latitude  changes.  Lines  of  constant  longitude,  called 
meridians,  are  normal  to  the  direction  in  which  longitude  changes.  The  distance  between  parallels 
is,  ignoring  the  slight  ellipdcity  of  the  earth,  constant  everywhere  on  the  globe.  The  distance 
between  meridians  varies  from  a  maximum  at  the  equator  to  zero  at  the  poles. 


52.  As  a  coDveniait  mneoionic,  one  nautical  mile  is  slightly  larger  than  the  statute  mile  and  is  defined  as  the 
arc  leagdi  subtended  at  sea  level  of  an  angle  of  <me  arc  minute  measured  from  the  earth’ s  center.  This  makes 
many  navigation  calculations  simpler. 
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6.A2  Sh^ie  of  die  Earth 

The  earth  is  not  perfectly  sp^eiical^^.  Its  polar  radius  is  about  3432  nm  and  its  equatorial  radius  is 
about  3444  nm.  So  it  is  flattened  by  atout  12  miles.  In  mathematical  terms,  its  sh^  can  be 
tqiproximated  by  an  oblate  spheroid.  This  has  many  implications  to  earth  fixed  systems: 

•  the  length  of  the  latitude  is  not  precisely  constant  over  the  globe 

•  the  radius  of  the  earth  is  not  constant  over  the  globe 

•  the  theoretical  gravity  field  is  slightly  different  than  for  a  spherical  earth 


There  have  been  many  models  of  the  earth  in  use  throughout  history,  and  for  historical  reasons, 
different  countries  continue  to  use  slightly  different  models.  The  definition  of  latitude  varies 
slightly  depending  on  which  reference  eUipsoid  is  used  by  the  country  in  which  the  system 
operates.  Some  20  slightly  different  ellipsoids  have  been  suggested  as  various  times  to  represent 
tte  sh^  of  the  earth.  They  differ  by  about  1  kilometer  in  earth  radius  at  the  equator  or  0.01%.  The 
DOD  regularly  charges  the  World  Geodetic  System  Committee  to  establish  a  standard  earth  model 
for  moping  purposes.  WGS'84  and  WGS-72  are  commonly  used  standards. 

d.43  UTM  Coordinates 

A  coordinate  system  called  Universal  Transverse  Mercator  or  UTM  coordinates  has  been  in 
common  use  by  the  army  for  some  time.  This  system  is  commonly  used  to  report  position  in  land 
navigation  inertial  systems  because  most  were  originally  developed  for  army  use.  The  utility  of 
UTM  coordinates  is  that  it  allows  the  system  output  to  be  correlated  simply  with  position  on  a  flat 
map. 

The  nearly  spherical  earth  is,  in  mathematical  terms,  a  nondevelopable  surface.  This  implies  that 
any  means  of  flattening  it  out  onto  a  plane  must  involve  some  distortion.  While  it  is  not  commonly 
£q>preciated,  the  Mercator  projection  which  is  used  tojgeirerate  a  flat  mtq>  of  the  worid  gives  rise  to 
massive  distortions  in  the  sluqre  of  the  earth’s  surface  .  The  Mercator  projection  is  constructed  by 
projecting  the  earth’s  surface  onto  a  cylinder  which  is  tangent  to  the  earth  at  the  equator. 

The  transverse  Mercator  projection  uses  a  cylinder  which  is  tangent  along  a  meridian.  UTM 
coordinates  in  the  north  and  east  directions  are  referred  to  as  noilliiiig  and  easting,  and  are  often 
reported  in  meters. 


53 .  It  was  Isaac  NewUHt  wbo  first  showed  that  die  eatb  must  bulge  at  the  equator  due  to  die  centrifugal  force 
ofitsroiatioii. 

54.  A  remarkable  example  is  tbm  Greenland  sqipears  to  be  larger  than  South  America  OD  an  everyday  mig) 
the  world.  In  fact,  Gremdand  is  mly  1/10  the  size  of  South  America.  Ihe  distortions  the  Mercator  prajec- 
don  arc  largest  near  the  ptdes,  SO  Greenland  is  particulariy  illustiative  of  this. 
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ibi  InitMiluatiftn 

When  an  INS  is  started,  a  sequence  of  operations  is  performed  which  takes  at  least  several  minutes 
of  time  before  the  system  is  operational. 

6.5.1  Start-up  and  Warm-up 

Mechanical  setup  involves  spinning  up  the  gyros  and  waiting  for  the  components  to  reach 
operating  temperature. 

6.5Jt  Trim 

Drift  trimming  is  the  process  of  calibrating  out  the  constant  components  of  gyro  drift  error  by 
trimming  the  biases  of  the  torquers.  When  the  system  is  known  to  be  stationary,  rotational  motion 
of  the  platform  is  known  based  on  the  earth’s  rate.  Any  rates  measured  above  that  expected  because 
of  the  earth’s  rotation  can  be  removed  in  this  way. 

6.5.3  Alignment 

The  determination  of  the  orientation  of  the  accelerometers  with  respect  to  the  earth  requires  the 
knowledge  of  their  orientation  with  respect  to  earth  fixed  vectors.  If  a  single  vector  were  used, 
platform  orientation  would  not  be  constrained  in  rotation  about  that  vector.  Hence,  two  noncolinear 
vectors  are  required.  These  vectors  may  be  obtained  from  star  sights,  surveyed  landmarks,  the 
gravity  and  magnetic  fields,  and  the  ea^’s  spin  axis.  In  general,  systems  which  report  relative 
position  still  requite  precise  alignment  in  order  to  avoid  incorrect  compensation  for  the  apparent 
forces  due  to  the  earth’ s  spin. 

In  tq)plications  where  the  start  point  is  always  the  same,  precision  surveyed  landmarks  provide  an 
exceUent  alternative.  Systems  that  are  required  to  align  themselves  whUe  launched  from  another 
flying  vehicle  or  on  orbit  can  make  use  of  star  sights  and  the  earth’s  horizon  as  landmarks,  or  can 
tnmsfer  aligrunent  from  the  host  vehicle.  The  variation  of  the  earth’s  magnetic  field  malces  it  an 
unreliable  heading  reference  unless  it  is  known  to  a  few  seconds  of  arc  at  the  start  point.  For 
terrestrial  mobile  robots,  self  aligrunent  using  only  the  sensors  of  the  inertial  system  itself  is  a 
convenient  method  of  aligrunent 

In  self  afignment,  the  local  vertical  is  determined  through  the  gravity  vector  and  the  earth’s  spin 
provides  a  reference  for  north.  For  an  earth  stabilized  system,  the  process  occurs  in  two  steps  called 
levdUiig  and  north  alignment  or  gyrocompassing^^. 


55.  Altbough  the  term  is  used  commonly,  it  does  ikx  mean  that  the  gyro  q;)etates  acctnding  to  the  iHindide 
the  gytocompass.  It  means  that  it  does  the  same  thing  as  the  gyioc(xiq)ass,  namely,  finds  north. 
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6i.5.4  Levdling 


During  levelling,  for  a  gimballed  system,  the  platform  is  controlled  by  a  special  levelling  servo 
which  drives  the  platform  to  rotate  into  the  gravity  vector  until  both  horizontal  accelerometers 
cease  to  measure  a  component  of  gravity.  At  this  point,  the  vertical  accelerometer  is  aligned  with 
the  local  vertical  and  the  platform  is  level. 

6.5,5  Gyrocompassing 

Gyrocompassing  is  accomplished  by  rotating  the  now  level  platform  about  the  vertical  until  one 
gyro  senses  no  component  of  the  earth’s  rotation^^.  This  rate  goes  to  zero  at  the  poles,  so  such 
systems  have  difficdty  during  initialization  near  polar  latitudes.  For  systems  operating  in  the  air, 
alignment  can  be  performed  in-flight  using  an  externally  supplied  velocity  signal.  Space  systems 
can  use  star  trackers  and,  if  in  orbit,  scatming  of  the  earth’s  horizon  to  align  themselves.  Sighting 
of  known  stars  can  also  be  used. 

Strapdown  systems  can  compute  the  vehicle  pitch  and  roll  from  the  appropriate  accelerometer 
outputs  given  knowledge  of  the  initial  position  on  the  earth  since  this  ^es  the  direction  of  the 
gravity  vector  with  respect  to  the  earth.  While  the  vehicle  is  stationary,  the  accelerometer  outputs 
measure  only  the  components  of  the  gravity  vector  and  its  resolution  onto  the  body  axes  constrains 
the  vehicle  orientation. 

For  strapdown  systems,  since  the  system  is  stationary,  the  angular  rate  vector  measured  by  the 
gyros  represents  the  resolution  of  the  earth  rate  onto  die  vehicle  axes,  so  the  vehicle  heading  can 
be  computed. 

The  gyrocompassing  step  is  complicated  by  the  low  signal  to  noise  ratios  encountered  when 
measuring  the  relatively  small  earth  rate.  Systematic  gyro  drift  must  also  be  eliminated  for  the 
technique  to  be  practical.  A  furflier  reduction  in  signal  to  noise  results  from  the  variation  of  the 
horizontal  component  of  the  earth  spin  rate  by  the  cosine  of  latitude.  This  is  the  second  aspect  of 
the  pole  problem  for  inertial  systems. 


56.  The  earth’s  magnetic  field  is  not  areliable  indicator  of  north  for  such  aprocess,  so  the  obvious  alternative 
is  not  used. 
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7  Performance 


The  relevant  performance  parameters  of  an  INS  include  accuracy  and  data  rate,  initialization  time, 
frequency  of  required  external  updates,  and  the  latitude  band  in  which  it  can  be  used.  These  factors 
and  others  are  discussed  below  from  the  point  of  view  of  providing  a  general  appreciation  of  the 
spectrum  of  systems  available. 

’LL  Ineilial  Grade  Componentry 

The  realization  of  inertial  navigation  had  to  await  the  historical  development  of  instrumentation 
technologies  which  could  supply  the  required  performance  parameters.  An  investigation  of  the 
typical  requirements  placed  on  gyros  and  accelerometers  provide  an  insight  into  the  severe 
difficulty  of  the  problems  that  have  been  solved  -  and  the  remarkable  performance  of  today’s 
systems. 

Positioning  accuracy  in  dead  reckoning  systems  is  particularly  sensitive  to  angular  drift  error^^. 
For  a  locally  level  gimballed  system,  consider  that  one  minute  of  arc  of  angular  platform  error 
corresponds  to  one  nautical  mile  of  position  error  on  the  earth’s  surface.  In  order  to  achieve  the 
modest  accuracy  requirement  of  one  mUe  of  position  error  per  hour  of  operation,  gyroscope  drift 
must  be  held  within  one  arc  minute  per  hour.  This  corresponds  to  one  thousandth  the  rotation  rate 
of  the  earth  or  one  milli  earth  rate  unit,  known  as  a  mem.  These  speciftcations  are  three  orders  of 
magnitude  more  stringent  than  those  imposed  on  gyros  used  in  conventional  aircraft 
instrumentation. 

For  accelerometers,  since  bias  is  multiplied  by  time  squared  in  deduced  reckoning,  requirements 
for  one  thousandth  of  a  g  of  sensitivity  are  common.  Yet  at  the  same  time,  the  device  must  also  be 
capable  of  measuring  the  earth’s  gravity.  Hence  wide  dynamic  range  and  extreme  sensitivity  are 
required  for  practical  systems. 

'LL  Aided  Inertial  Systems 

While  an  inertial  navigation  system  is  conceptually  a  self-contained  device,  it  is  typically  the  case, 
for  reasons  discussed  above,  that  serious  degradation  of  performance  arises  unless  redundant 
external  measurements  are  used  to  remove  the  er  s  that  build  up  over  time  in  a  marmer  typical  of 
all  deduced  reckorung  systems. 

The  use  of  external  measurements  in  this  marmer  is  referred  to  as  damping,  and  the  systems 
employing  the  technique  are  called  aided  or  augmented  inertial  systems.  In  this  view  of  an  INS, 
its  real  utility  is  in  its  ability  to  provide  continuous  positioning  information  between  external 
updates. 

Qassically,  INS  ouq)uts  were  adjusted  to  agree  perfectly  with  external  measurements  at  the  point 
in  time  when  the  external  measurement  was  made.  The  modem  technique  for  incorporating 
exterrud  measurements  is  optimal  estimation,  often  incorporating  the  Kalman  filter.  Using  this 
techrtique,  estimates  of  the  errors  in  both  the  INS  generated  navigation  quantities  and  the  errors  m 


57.  A  good  iUustratioa  of  this  is  that  a  drift  rate  of  a  mete  1/10  of  a  degree  per  hour  will  cause  a  transatlantic 
flight  to  end  up  in  Paris  rather  than  London. 
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the  external  measurements  are  combined  to  provide  the  best  overall  estimate  of  the  navigation 
quantities. 

The  damping  of  the  vertical  channel  can  be  accomplished  in  many  ways.  Measurements  of 
position,  velocity,  and  attitude  can  be  used  singly  or  in  any  combination.  The  particular  solution 
adopted  is  application  dependent 

For  flying  systems,  altitude  is  typically  available  from  other  sources.  Barometers  can  be  used  for 
only  coarse  measurement  of  altitude  since  air  pressure  varies  with  weather  conditions  as  well  as 
altitude.  Altitude  can  be  determined  to  a  few  hundred  feet  in  this  way.  Radar  altimeters  are  often 
accurate  to  a  few  percent  of  altitude.  However,  many  flying  systems  are  less  concerned  with  the 
accuracy  of  altitude  than  with  horizontal  position.  Doppler  radar  velocity  inputs  are  also  commonly 
used  for  airborne  systems. 

Seaborne  systems  can  take  advantage  of  the  extensive  radio  navigation  systems  that  have  been 
available  for  many  decades.  Space  systems  can  use  star  sights,  and  when  in  earth  orbit,  horizon 
measurement  and  navigation  satellites.  The  latter  are  also  particularly  useful  for  all  terrestrial 
systems  that  have  clear  view  of  the  satellites. 

Land  navigation  systems  can  use  satellite  systems  as  well  as  landmarks  and  map  matching  to  fix 
position.  Periodic  return  to  the  same  position  allows  position  errors  to  be  zeroed  by  telling  the 
system  that  it  has  returned  to  the  base  position.  Similarly,  the  technique  of  zero  velocity  update 
notifies  the  system  when  it  has  stopped  moving.  Odometry  damping  is  also  a  common  technique 
in  land  navigation  systems. 

lA.  Sources  of  Error 

In  practice,  the  most  significant  sources  of  errors  in  an  INS  are  limited  to  nonidealities  in  the 
accelerometers  and  gyros,  and  several  kinds  of  alignment  errors.Together,  these  often  amount  to 
90%  of  the  total  system  error.  A  more  complete  list  of  error  sources  is  as  follows: 

7.3.1  Accelerometer  Errors 

•  Bias  -  the  output  of  the  sensor  for  zero  input 

•  Scale  Factor  -  deviation  of  the  slope  of  the  device  calibration  curve  from  ideal 

•  Cross  Coupling  -  rebalance  servo  droop 

•  Vibropendulosity  -  steady  state  output  for  vibratory  input 

7.3.2  Gyro  Errors 

•  bias  -  parasitic  torques  caused  by  restraints  and  magnetic  and  thermal  effects 

•  proportional  -  torques  caused  by  mass  unbalance  and  anisoelasticity 

•  anisoelastic  -  deformations  normal  to  applied  loads 

•  torquing  error  -  errors  in  the  torques  supplied  to  rotate  the  platform 

•  random  drift  -  uncompensated  nondeterministic  errors 

7.33  Platform  At^nment  Errors 

•  initial  level  misalignment  of  platform  -  levelling  error 
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•  initial  azimuth  misalignment  of  platform  -  gyrocompassing  error 

•  dynamic  misalignment  of  platform  -  servo  droop  in  platform  controller 

•  accelerometer  misaligiunent  -  with  respect  to  platform  axes 

•  gyro  misalignment  -  with  respect  to  platform  axes 

7.3.4  Computational  Stabilization  (equivalent  to  gyro  drift) 

•  truncation  -  finite  precision  computations 

•  roundoff  -  finite  precision  computations 

•  commutation  -  small  angle  assumptions  in  forming  the  direction  cosine  matrix 

•  pickoff  -  finite  resolution  of  sensor  interfaces 

7.3.5  External  Errors 

•  errors  in  geodetic  information 

•  errors  in  damping  signals 


2A.  Configuration  Of  Inertial  Systems 

In  addition  to  the  complete  inertial  navigation  system  supplying  direct  navigation  information, 
some  vendors  supply  strapdown  IMUs  as  stand  alone  units.  These  units  output  body  coordinate 
referenced  velocity  and  orientation  information  which  must  be  converted  by  the  user  into  a  more 
useful  form.  Specifications  for  these  systems  do  not  differ  significantly  from  those  applying  to  the 
more  complete  systems. 

Many  systems  admit  velocity  and/or  odometry  damping  inputs  which  are  processed  internally  by 
the  INS  in  order  to  Kalman  filter  them  with  die  computed  motion.  Hybrid  INS/GPS  systems  are 
also  available. 

2jL  Physical 

Systems  can  weigh  as  much  as  50  lbs,  occupy  as  much  as  1  cu  ft.  of  volume,  and  consume  up  to 
100  watts  of  power.  Strapdown  systems  and  the  use  of  ring  laser  gyros  are  becoming  popular.  The 
lack  of  gimb^  structure  in  these  systems  affords  reduction  in  size,  weight,  power  consumption,  and 
cost 

Strapdown  systems  often  use  miniature  gyro  and  accelerometer  components  which  weigh  a  few 
grams  and  occupy  only  a  fraction  of  a  cubic  inch  of  volume.  The  bulk  of  the  volume  of  the 
strapdown  system  is  taken  up  by  the  processing  electronics. 

2Jl  Data  Rates 

The  rate  at  which  systems  can  supply  positioning  data  can  vary  significantly.  Up  to  several  hundred 
Hz  update  rates  are  possible,  depending  on  the  processing  electronics  and  the  data  interface. 

2iL  Initializatwn 

The  time  a  system  takes  from  power-up  to  completion  of  alignment  can  vary  from  a  few  minutes 
to  30  minutes  depending  on  ambient  temperature  and  local  latitude.  Vendors  often  quote 
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worldwide  operation  below  75  degrees  latitude.  Warm-up  may  require  elevated  power 
consumption  up  to  500  watts. 

Newer  systems  integrated  with  GPS  can  function  in  a  moving  base  alignment  mode  where 
successive  GPS  readings  are  used  to  align  the  real  or  computational  platform  as  the  vehicle  moves. 

2.1  Acgwraty 

The  accuracy  of  the  INS  can  be  specified  in  several  ways,  depending  upon  the  application  it  is 
suited  for.  For  cruise  vehicles,  the  rate  at  which  the  position  error  builds  up  over  time  is  specified. 
For  ballistic  and  space  vehicles,  the  position  error  at  the  target  or  the  velocity  error  at  propulsion 
cut-off  is  given.  Lwd  vehicle  systems  often  specify  accuracy  as  a  percentage  of  distance  travelled. 

Today’s  commercial  cruise  systems  can  achieve  accuracies  on  the  order  of  0.2  nautical  miles  of 
error  per  hour  of  operation.  Pitch  and  roll  are  often  accurate  to  0.05®  and  true  geographic  heading 
to  0.5®.  In  some  cases,  position  accuracy  along  the  trajectory  (alongtrack)  and  both  normal 
directions  (crosstrack  and  vertical)  are  distingiiished. 

Land  vehicle  navigation  systems  achieve  0.2%  to  2%  of  distance  travelled.  Pitch  and  roll  can  be 
measured  to  0.1®  and  heading  to  0.5®.  While  the  inertial  components  typically  perform  better  than 
this  by  an  order  of  magnitude,  systems  require  shock  and  vibration  isolation  mounting.  These 
isolation  systems  necessarily  permit  relative  motion  between  the  system  housing  and  the  vehicle 
chassis  which  significantly  reduces  the  accuracy  of  the  orientation  output. 

At  present,  strapdown  systems  do  not  compete  with  gimballed  ones  in  terms  of  accuracy,  but  this 
is  likely  to  change  as  component  technologies  improve. 
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PART  IV:  Satellite  Navigation 

Systems^* 


Satellite  navigation  was  developed  to  answer  the  military  need  to  precisely  determine  the  position 
of  air,  sea,  and  land  vehicles.  It  is  rapidly  becoming  a  definitive  navigation  tool  since  it  provides 
continuous,  high  accuracy  positioning  anywhere  on  the  surface  of  the  planet  and  the  near  space 
region,  24  hours  a  day,  under  all  weather  conditions.  The  receiver  systems  are  small,  lightweight, 
and  easy  to  use,  and  have  made  handheld  global  positioning  systems  a  reality. 

This  technology  is  important  for  mobile  robots  for  the  same  reasons  it  is  important  to  the  military. 
For  a  mobile  field  robot  operating  over  large  distances,  or  near  the  poles,  it  is  one  of  only  a  few 
available  good  solutions  to  the  position  estimation  problem. 

The  military  has  used  the  TRANSIT  constellation  of  satellites  since  1964  and  the  resource 
industries  have  used  the  Starfix  system  since  1986  for  such  applications  as  locating  drilling  rigs, 
seismic  surveys,  and  pipeline  construction.  Both  of  these  systems  ate  now  obsoleted  by  two  more 
recent  systems,  the  Global  Positioning  System  developed  by  the  United  States  Department  of 
Defease,  and  its  Soviet  counterpart,  called  Glonass.  The  next  generation  of  INMARSAT  satellites 
is  also  planned  to  radiate  GPS  like  signals  for  positioning. 

All  of  these  satellite  navigation  systems  are  based  on  identical  principles,  notably  range 
triangulation,  and  GPS  and  Glonass  ate  virtually  identical  from  the  perspective  in  this  report.  For 
this  reason,  only  GPS  will  be  discussed  here  but  the  reader  is  cautioned  that  virtually  ^  of  the 
information  in  this  chapter  applies  equally  well  to  Glonass. 

Since  GPS  provides  a  positioning  fix,  it  complements  inertial  navigation  or  other  dead  reckoning 
techniques.  Integration  of  both  techniques  achieves  performance  which  is  an  improvement  over 
what  would  be  possible  using  either  system  alone. 

There  are  many  navigation  signals  available  from  the  satellites  and  many  different  ways  of 
processing  lliem.  These  processing  techniques  give  rise  to  trade-offs  between  the  accuracy  of  the 
fix,  how  quickly  it  can  be  updated,  and  the  range  over  which  the  accuracy  can  be  maintained. 
Accuracies  vary  fit>m  tenths  of  a  kilometer  to  tenths  of  a  centimeter  depending  on  the  signals  used 
and  how  they  ate  processed.  The  systems  can  be  used  to  measure  position,  veloci^,  attitude,  and 
very  precise  time. 
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1.  Principles  of  Operation 

A  constellation  of  satellites  is  maintained  in  earth  orbit  and  radio  receivers  at  or  near  the  surface  of 
the  planet  are  used  to  decode  the  transmissions  of  the  satellites  to  compute  from  them  the  motion 
and  position  of  the  receiver.  GPS  was  developed  primarily  as  a  military  navigational  aid,  so  it  was 
(tesigned  as  a  one  way  broadcast  system.  Receivers  do  not  have  to  transmit  any  signals  back  to  the 
satellites  that  might  give  away  their  position. 

This  aspect  of  the  design  has  a  number  of  other  implications.  Satellites  receive  no  signals,  so  there 
is  no  satellite  capacity  limit  on  the  number  of  users.  Receivers  transmit  no  signals,  so  anyone  with 
knowledge  of  the  communications  protocol  can  use  the  system  without  being  detected.  Civilian  use 
of  GPS  has  been  provided  for  in  the  design  and  there  is  no  charge  associated  with  satellite  signal 
reception.  The  eventual  number  of  civilian  usei^  is  estimated  to  be  in  the  millions. 

Pffsitiftn  MMSurgnsnt 

Positioning  is  based  on  the  principle  of  range  trianguladon.  The  receiver  needs  to  know  the  range 
to  the  satellites  that  ate  being  used  and  the  positions  of  these  satellites.  It  can  determine  its  own 
position  from  only  this  information.  As  in  all  triangulation  situations,  the  number  of  independent 
observations  required  depends  on  how  many  dimensions  of  position  are  required  (2D  or  3D)  and 
how  much  information  is  already  known. 

The  two  dimensional  case  is  indicated  in  Figure  25.  Two  satellites  ate  required  in  general^^.  If  the 
positions  of  the  satellites  ate  known,  and  the  ranges  from  the  receiver  to  the  satellites  can  be 
measured,  the  receiver  position  is  constrained  to  be  one  of  the  two  points  of  intersection  of  the 
circular  lines  of  position^.  After  the  two  simultaneous  equations  are  solved,  a  rough  initial 
estimate,  or  a  third  satellite  can  be  used  to  resolve  which  of  the  two  solutions  applies. 


Figure  25  Two  Dimensional  Range  Triangulation 


39.  The  receiver  can  also  wait  until  one  satellite  moves  to  a  new  positioa,  bat  see  die  taler  sectioD  00  GDOP. 
60.  If  the  receiver  is  on  a  line  between  the  satellites,  only  one  point  (tf  imersection  occurs. 
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The  receiver  position  is  given  by  the  simultaneous  solution(s)  (x,y)  to  the  following  two  equations: 


where  the  satellite  positions  (xl,y  1),  (x2,y2)  and  ranges  rl  and  r2  are  considered  known  quantities. 

The  three  dimensional  case  is  similar  except  that  three  satellites  are  used  and  the  lines  of  position 
are  spheres  centered  at  the  satellites.  Two  spheres  intersect  in  a  circle  and  another  sphere  intersects 
this  circle  in  at  most  two  places.  The  ambiguity  can  be  resolved  with  another  satellite  or  an  initial 
estimate.  For  vehicles  whose  altitude  is  accurately  known,  only  two  satellites  are  required. 

The  satellite  positions  are  known  as  the  ephemeris  data.  This  information  is  broadcast  to  the 
receivers  from  orbit  The  second  required  piece  of  information,  the  ranges,  are  det^mined  by 
measuring  the  time  required  for  the  propagation  of  the  signal  from  each  satellite  to  die  receiver, 
and  multiplying  by  the  radio  wave  speed^  .  To  make  this  possible,  each  satellite  also  broadcasts 
the  time  along  wiA  its  position.Then  the  receiver  can  subtract  the  time  that  the  message  was  sent 
from  the  time  that  it  was  received. 

The  receiver  performs  the  very  difficult  task  of  measuring  a  range  of  about  50,000  kilometers  to 
the  nearest  meter  or  so.  Iherefore,  an  error  of  only  one  part  in  fifty  million  in  either  the  propagation 
time  or  the  speed  of  light  will  cause  an  error  of  one  mete:  in  the  observed  range.  Prediction  of  the 
wave  speed  and  tremendously  precise  measurement  of  time  are  both  key  design  features  of  satellite 
navigation. 

12.  Timt  Mtagyremmt 

In  order  for  a  receiver  to  compute  the  propagation  time  directly,  it  must  have  its  clock  synchronized 
with  the  clocks  used  on  the  satellites,  and  all  of  the  satellites  must  be  similarly  synchroitized  among 
tltenselves.  This  is  because  every  nanosecond  of  time  error  causes  one  foot  of  range  measurement 
error.  It  se«»ns  at  first  glance  that  all  satellites  and  all  receivers  everywhere  must  all  agree  on  the 
time  to  the  neatest  nanosecond. 

For  die  satellites,  this  synchronization  is  actually  achieved  in  practice  through  the  use  of  very 
precise  atomic  clocks.  However,  it  is  impractical  for  a  receiver  to  use  an  atomic  time  standard. 
Che^  crystal  oscillators  are  used  instead.  These  clocks  are  precise  enough  but  not  accurate 
enough,  so  they  may  be  offset  in  time  from  the  satdlites.  The  solution  to  the  problem  is  to  consider 
time  a  fourdi  unknown  m  the  ranging  process,  and  to  use  a  fourth  satellite  to  provide  an  extra 
constraint  equatioiL  If  alnuide  is  known  accurmely  (as  it  is  on  the  sea)  one  less  satellite  is  needed. 


61.  NombMlly,  ihU  is  the  qwed  of  ligbL  While  this  is  a  well  known  ooosiant  in  vacuum,  tbe  wave  q)eed 
throiqd^  the  siinoqphere  varies  slighdy  with  tme  and  atmoqiiieiic  ooDStitiients. 
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The  time  offset  of  the  receiver  is  called  the  user  clock  Mas.  The  effect  of  this  bias  is  shown  in 
Figure  26.  The  satellites  are  synchronized,  so  the  user  clock  bias  gives  rise  to  ranges  to  all  of  them 
which  are  in  error  by  the  same  amount  These  erroneous  ranges  are  called  pseudoranges. 

If  the  user  clock  is  slow,  considering  the  two  dimensional  case,  the  three  ranges  will  create  a 
roughly  triangular  region  as  shown.  Offsetting  the  sides  of  the  triangle  by  fixed  amounts  on  all 
sides,  the  receiver  can  determine  the  clock  bias  as  an  extra  unknown. 


Figure  26  Computing  the  User  Position  and  Clock  Offset 


In  duee  dimensions,  the  receiver  position  is  given  by  the  simultaneous  solution(s)  (x,y,z)  to  the 
following  four  equations. 


rl  =  J(x-xl)^+ (y-yl)^+ (z-zl)^  +  cAt 
r2  =  J(x-x2)^+  (y-y2)^+  (z-z2)^  +  cAt 

Equation  39 

r3  =  J(x  -  x3)  ^  +  (y  -  y3)  ^  +  (z  -  z3)  ^  +  cAt 

r4  =  J(x-x4)^  +  (y-y4)^+  (z-z4)^  +  cAt 

Althou^  die  equations  are  written  in  terms  of  geocentric  cartesian  coordinates,  the  answer  can  be 
ctmvermd  into  latitude,  longitude,  and  altitude  on  any  reference  ellipsoid  model  of  die  eardi^,  or 
into  grid  coordinates  on  a  map.  The  equations  are  solved  in  practice  by  a  standard  technique  for 

62.  There  are  many  models  in  use  wiiicb  represent  the  eaitb  as  an  ellqisoid.  Hence,  ibe  latitude,  loogiiiide, 
and  aldnide  (tf  a  point  depend  on  which  model  is  in  use. 
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solving  simultaneous  nonlinear  equations  called  Newton-Raphson,  and  if  more  than  four  satellites 
are  used,  a  least  squares  technique  is  used. 

Since  receivers  can  measure  their  clock  offsets  to  great  accuracy,  some  GPS  receivers  operate 
solely  as  very  precise  clocks  or  as  frequency  standards. 

L2l  WaYg  Sp^  McasBrcmfflt 

The  measurement  of  wave  speed,  or  equivalently  atmospheric  delay,  is  accomplished  in  receivers 
in  one  of  two  different  ways.  Some  receivers  use  a  mathematical  model  of  the  atmosphere  to 
compute  the  delay.  The  coefficients  in  this  model  are  broadcast  regularly  by  the  satellites 
themselves. 

Siuce  the  effect  depends  in  a  known  way  on  frequency,  some  receivers  are  able  to  actually  measure 
the  difference  between  each  of  the  two  navigation  signal  carriers  transmitted  by  each  satellite.  This 
technique  is  far  more  accurate  than  the  use  of  a  mathematical  model.  Atmospheric  delay  is  one  of 
the  most  significant  sources  of  error. 

lA.  ydftaty  Mtasuremait 

The  principle  used  to  measure  velocity  is  the  doppler  frequency  shift  The  doppler  shift  for  the 
frequency  of  each  satellite  is  a  direct  measure  of  the  relative  velocity  of  receiver  and  satellite  along 
the  line  between  them.  Each  satellite  has  a  very  high  velocity  relative  to  a  stationary  receiver 
because  of  both  the  orbital  motion  of  the  satellite  and  the  rotation  of  the  receiver  with  the  earth. 
The  velocity  solution  is  obtained  by  differentiating  the  four  dimensional  navigation  solution  of 
Equation  39.  These  are  the  equations  solved  in  the  receiver  to  deteitnine  velocity. 

LL  Orientation  Measiiroient 

The  principle  used  to  measure  orientation  of  a  rigid  body  is  to  measure  the  differential  three 
dimrasional  positions  of  different  points  on  die  body.  In  practice,  measurements  of  the  radio  carrier 
phase  are  used  to  achieve  the  positioning  accuracies  required  for  acceptable  orientation 
measurement  Very  accurate  real-time  orientation  measurement  for  rigid  vehicles  has  been 
achieved. 

The  accuracy  of  angular  measurement  increases  with  the  length  or  basdine  distance  between  the 
points  positioned.  Using  four  satellites,  the  user  clock  error  and  the  3D  position  of  each  point  can 
be  calculated.  If  baselines  are  known  accurately  enough,  less  satellites  are  required.  At  least  three 
points  must  be  measured  to  determine  all  three  angles  of  rotation. 
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2  Implementation 

The  navigation  signals  used  by  GPS  satellites  will  seem  oveily  complex  to  the  reader  unless  the 
reasons  for  this  complexity  are  presented  Hrst  The  GPS  signals  owe  their  comptexity  to  the 
requirements  that  the  system  was  designed  to  meet  For  example,  there  are  two  carrier  signals  and 
two  different  modulation  signals.  Further,  one  of  the  modulation  signals  spears  on  both  carriers 
while  the  other  appears  only  on  one.  A  third  modulation  appears  on  both  carriers  but  contains  very 
different  information  than  the  other  two.  This  section  will  explain  the  reasoning  behind  this  design. 

2mL,  Siyiial  Characteristics 

GPS  was  intended  to  meet  both  a  military  and  a  civilian  demand  for  global  all  weather  navigation. 
Satellite  radio  sources  are  particularly  suited  for  diis  purpose  since  they  can  radiate  their  signals  to 
every  point  on  the  siuface  of  the  planet  It  is  intend^  to  be  a  multiple  simultaneous  access 
positioning  utility,  available  to  as  many  users  as  there  are  receivers.  Further,  receivers  are  required 
to  be  small  enough  to  carry  in  one  hand.  The  signals  are  intended  to  provide  highly  accurate 
measurements  of  both  position  and  velocity  while  at  the  same  time  being  easy  to  acquire  quickly. 

2.1.1  Carriers 

GPS  signals  are,  like  most  radio  signals,  a  modulated  carrier  signal.  The  carrier  signals,  being 
predictable  sinusoids,  cannot  cany  information  themselves.  Information  is  carried  by  the 
modulation  signals  which  are  mixed  with  the  carriers. 

It  is  known  that  the  radio  wave  delay  through  the  atmosphere  varies  with  time  and  position  by  as 
much  as  100  memrs  in  equivalent  range.  Mathematical  models  do  not  remove  most  of  this  error, 
so  the  design  includes  two  separate  L  band  carriers,  denoted  LI  and  L2.  Atmospheric  delay  varies 
with  the  square  of  the  carrier  frequency.  With  two  carriers,  the  differential  delay  of  both  can  be 
used  to  measure  the  absolute  delay. 

The  choice  of  L  band  carriers  represents  a  trade-off  between  the  limited  available  bandwidth  in 
UHF  and  the  excessive  atmosphoic  absorption,  called  qpace  loss,  of  C  band.  The  L  band  fiiequ^icy 
is  high  enough  to  accommodate  the  20  MSlz  bandwidth  required  and  the  centimeter  wavelengdi 
enables  precision  velocity  measurement 

2.1J2  Modulatitm  Signals 

The  tequiremoit  for  small  receivers  translates  into  a  requirement  for  small  antennae  ^  on  the 
receiver,  which  in  turn  means  diat  they  must  woiit  at  very  small  signal  levels.  In  order  to  achieve 
this,  both  modulation  signals  are  based  on  a  particular  l^d  of  digital  pseudorandom  noise  code 
called  a  Gold  code.  These  PRN  codes  allow  receivers  with  antenna  oidy  a  few  inches  across  to 
extract  vary  low  power  signals  from  the  noise  by  correlating  than  with  expectations. 


63.  Ooottast  this  widi  the  large  parabolicreflectofs  needed  to  cooceotiate  satdlite  TV  broadcasts.  GPS  signals 
are  loa^ily  die  iwwie  strength  as  saidliie  TV  broadcasts.  Even  if  a  GPS  receiver  nsed  such  a  dish,  it  would 
need  fota  them  to  shnnhaoeoasly  track  four  satellites.  The  Starfix  satellite  positioning  system  uses  a  48 
inch  antenna  and  is  hence  not  vislUe  in  many  applkatioos  siqjpotted  by  GPS. 
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PRN  codes  have  two  other  advantages.  They  are  particularly  suited  to  the  multiple  access 
requiremenL  Each  satellite  transmits  a  different  code,  so  the  receiver  can  distinguish  between  the 
signals  of  all  satellites  which  arrive  at  its  antenna  simultaneously.  Also,  PRN  codes  are  particularly 
impervious  to  deliberate  or  unintentional  jamming,  so  they  support  the  need  for  secure  military 
access. 

There  are  two  modulation  signals  used.  The  civilian  signal,  called  C/A  for  clear  access  or  coarse 
acquisition,  has  a  300  meter  pulse  wavelength.  The  military  signal  is  designated  P  for  precise,  and 
is  called  Y  when  it  is  encrypted.  It  has  a  pulse  wavelength  of  30  meters.  Pseudorange  resolution  is 
directly  related  to  the  wavelength  of  the  code,  so  the  civilian  code  has  a  lower  intrinsic  precision. 
The  civilian  C/A  code  is  mixed  with  the  LI  carrier  while  the  P  code  is  mixed  with  both  carriers  and 
broadcast  redundantly. 

The  use  of  two  codes,  one  of  which  is  encrypted,  gives  the  military  access  to  a  completely  separate, 
secure,  positioning  system.  A  second  reason  for  using  two  codes  is  that,  since  the  P  code  is  difficult 
to  acquire,  the  C/A  was  designed  for  easy  acquisition.  Once  the  C/A  code  is  acquired,  the  P  code 
is  easier  to  acquire.  Code  acquisition  is  discussed  in  the  next  section. 

2.13  Nav^ation  Message 

In  addition  to  the  PRN  codes,  a  second  block  of  information  is  digitally  broadcast  on  both  carriers. 
This  navigation  message  is  transmitted  at  the  very  slow  rate  of  SO  bits  per  second,  and  is  updated 
and  repeated  every  12.5  minutes.  The  navigation  message  contains  a  large  amount  of  information 
provided  to  assist  receivers  in  their  operation. 

The  system  time  of  the  week,  known  as  the  handover  word,  is  provided  to  assist  in  acquiring  the 
P  co<te  after  the  C/A  code  is  acquired.  The  transmitting  satellite  provides  precise  ephemeris  data 
on  its  own  orbit  for  use  in  the  navigation  solution.  Also,  less  accurate  ephemeris  data  for  all  other 
satellites,  known  as  the  almanac,  is  transmitted  to  permit  receivers  to  predict  when  new  satellites 
will  rise  above  the  horizon  and  become  avaihd}le.  The  ephemeris  data  is  described  by  the  classical 
ellipse  parameters  discovered  and  defined  by  Kepler. 

Atmospheric  delay  model  coefficients  are  provided  since  the  magnitude  of  the  delay  varies 
significantly  throughout  the  day.  Satellites  also  provide  information  on  their  own  health  and  the 
precision  of  range  measurements  that  can  be  expected  if  their  signals  are  used.  In  this  way, 
receivers  can  pick  an  optimal  set  of  satellites  in  view  for  dieir  computations.  The  complete  set  of 
GPS  modulation  signals  is  summarized  in  Figure  27: 


Figure  27  GPS  Signals 
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2mZ.  l2gar..Smimt 

The  GPS  user  segment  consists  of  the  receivers  in  use  on  the  surface  and  in  the  near  space  region. 
The  basic  functions  of  the  receiver  are  to  perform  code  correlation  in  order  to  measure  propagation 
time  and  clock  bias,  to  perform  compensation  for  known  errors,  and  to  implement  the  solution  to 
the  navigation  equations.  Receivers  internally  generate  exactly  the  same  codes  as  the  satellites. 
Time  of  arrival  measurement  is  accomplished  by  shifting  the  internal  codes  in  time  and  computing 
the  correlation  with  the  externally  received  codes  as  shown  in  Figrue  28. 


Figure  28  Receiver  Code  Correlation 

The  C/A  code  repeats  once  every  millisecond  while  the  P  code  period  is  one  week  long.  Tbe  C/A 
code  is  short,  so  the  amount  of  shifting  required  is  limited,  and  a  match  can  be  found  quickly.  A 
match  occurs  when  the  correlation  function  peaks.  Once  a  match  is  computed,  the  code  is  said  to 
be  acquired.  Once  acquired,  the  time  of  signal  propagation  is  given  by  the  time  shift  that  was 


necessary.  This  time  includes  the  clock  bias  which  is  removed  later.  The  next  step  is  to  derive  the 


wave  speed. 


Receivers  can  compensate  for  atmospheric  delay  to  compute  an  accurate  wave  speed  by  using 
either  mathematical  model  coefficients  transmitted  by  the  satellites,  or  by  using  two  carriers. 
Multiplying  this  by  the  propagation  time  gives  the  pseudorange  to  the  satellite.  When  this  process 
is  repeated  for  four  different  satellites,  the  receiver  can  solve  the  navigation  equations  and 
determine  its  position. 
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2i2.  Spatt  Stsmfflt 

The  GPS  space  segment  consists  of  a  constellation  of  18  earth  satellites^  and  three  spares  that 
circle  the  earth  in  each  of  six  orbits  as  shown  in  Figure  29.  Hie  six  orbital  planes  are  each  inclined 
at  SS°  to  the  equatorial  plane,  and  separated  by  60°  in  longitude.  Orbits  are  nearly  circular,  1 1,000 
miles  in  altitu^  and  repeat  exactly  twice  per  sidereal  day.  This  implies  that  the  satellites  are  in 
exactly  the  same  place  four  minutes  later  every  day®^. 


Figure  29  Orbital  Configuration 


The  constellation  is  designed  to  ensure  that  at  least  four  satellites  are  visible  at  any  time.  Satellites 
are  considered  visible  above  the  nominal  mask  ang^  of  5°-10°.  When  the  elevation  angle  above 
the  horizon  exceeds  the  mask  angle,  the  signals  are  strmig  enough  for  good  reception.  However, 
foliage  or  other  attenuation  may  increase  the  mask  angle  significantly. 

Use  of  GPS  requires  a  clear  line  of  sight,  and  the  signals  cannot  penetrate  water,  soil,  or  walls  very 
welL  Hence,  GPS  cannot  be  used  underwater,  in  thick  forest,  or  in  mines  and  tunnels.  Signals  can 
be  obstructed  by  tall  buildings,  or  by  parts  of  a  rotating  vehicle  if  the  receiver  antenna  has  not  been 
placed  carefully. 


64.  Aoconliiig  to  NC^AD,  Ibeie  are  cnncntly  aboiu  7000  detectable  objects  in  eaidi  otbiL 

65.  It  is  possibk  to  bave  a  saieOitB  qipear  to  te  fixed  over  the  equator  by  matcbing  its  oibital  period  witb  die 
eanh’ s  lotafion.  Ihe  oririial  penod  of  a  drculsr  orbit  is  a  fimcdcm  of  only  the  altitude.  A  geosyndtmnious  al¬ 
titude  is  mncfa  higher  than  that  of  the  GPS  satellites.  Geosynduonous  space  is  cqndly  filling  iq)  with  satellites. 
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2il  Ground  Segment 

The  GPS  ground  segment  consists  of  a  series  of  five  ground  stations  spaced  in  longitude  around 
the  globe.  One  of  the  stations  is  designated  the  Master  Control  Station  (MCS).  The  function  of  GPS 
Master  Control  is  to  track  the  positions  of  all  satellites  very  precisely  and  to  maintain  the  overall 
system  time  standard. 

Space  vehicle  positions  are  very  predictable  due  to  the  negligible  effect  of  atmosphere  at  orbital 
altitude.  Nevertheless  small  perturbing  forces  exist  which  are  caused  by  the  oblateness  of  the  earth, 
the  presence  of  the  moon,  and  even  the  pressure  of  solar  radiation  accumulated  over  time.  The 
effect  of  these  forces  is  to  perturb  the  satellite  orbits  from  their  ideal  Keplerian  elliptical  shape. 

Inverting  the  ranging  process,  considering  MCS  to  be  the  user  and  the  ground  stations  to  be  the 
satellite,  makes  it  possible  to  determine  the  position  of  each  satellite  very  accurately.  MCS  updates 
the  satellites  on  their  orbital  parameters  at  least  once  per  day  so  that  precise  information  can  be  later 
transmitted  to  the  receivers. 

Satellites  must  be  synchronized  in  time  to  a  precision  of  nanoseconds  to  avoid  the  range  errors  that 
arise  from  timing  errors.  To  do  this,  MCS  maintains  an  overall  system  time  standard  known  as  GPS 
system  time  through  the  use  of  highly  accurate  cesium  atomic  clocks^.  Also,  each  satellite 
incorporates  two  cesium  and  two  rubidium  atomic  clocks  for  redundancy. 

Atomic  clocks  are  used  since  the  more  common  crystal  oscillators  are  not  stable  enough  for  this 
application.  MCS  regularly  transmits  to  each  satellite  the  offset  of  its  clock  from  GPS  system  time. 
Satellite  clock  offsets  are  then  retransmitted  to  receivers  so  that  the  precise  GPS  system  time  of 
signal  emission  is  known.  Overall  constellation  control  is  indicated  in  Figure  30. 


66.  Iheatomk;  clock  is  not  based  on  the  same  processes  as  nuclear  energy.  Rather,  the  ability  of  the  attan  to 
resonate  as  it  changes  its  quantum  oiergy  states  is  e^qtloited.  In  the  case  of  cesium,  the  magn^  dipole  of  the 
sole  outennost  electron  is  swimhed  £rt»n  being  parallel  to  and  antiparallel  to  the  dipole  of  the  nucleus.  The 
{Hinc^les  used  are  similar  to  diose  of  a  laser  or  a  microwave  oven. 
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1  P.£rfftmiants 


The  perfonnance  of  GPS  can  be  characterized  along  three  independent  dimensions:  the  precision 
and  accuracy  of  the  fix,  the  rate  at  which  it  is  updated,  and  the  range  over  which  the  accuracy  can 
be  achieved.  For  robotic  surveying  applications,  it  is  reasonable  to  allow  the  vehicle  to  dwell  at  a 
point  of  interest  for  a  considerable  time  so  that  a  more  accurate  fix  is  achievable.  It  may  also  be 
feasible  in  certain  applications  to  postprocess  positioning  data  off-line  for  a  further  improvement 
in  accuracy,  provided  a  redundant  positioning  system  is  used  for  control  of  the  vehicle  itself. 

If  GPS  is  used  to  determine  vehicle  position  in  real  time  for  navigation  purposes,  high  accuracy 
must  be  traded  against  the  update  rates  necessary  for  real  time  control.  Further,  if  the  vehicle’s 
range  of  motion  is  restricted  by  the  application,  differential  techniques  can  be  used  to  improve 
accuracy  considerably,  at  the  cost  of  a  second  receiver,  a  second  communications  link,  and  the 
limited  range. 

The  accuracy  achievable  at  any  update  rate  is  also  a  function  of  the  signal  combinations  used  and 
the  processing  performed  on  them.  The  inherent  resolutions  of  the  C/A,  P,  and  carrier  signals  vary 
from  meters,  to  decimeters,  to  millimeters.  It  is  not  the  intention  here  to  provide  a  quantitative 
assessment  of  GPS  accuracy  in  any  of  its  modes  of  use.  Rather,  this  section  will  equip  the  reader 
to  understand  the  factors  that  affect  accuracy  to  a  degree  that  will  make  informed  design  decisions 
possible. 

2iL  Range  Resolution 

The  fundamental  range  resolution  available  in  any  mode  of  operation  is  related  to  the  signals  used. 
Specifically,  achievable  range  resolution  is  a  fraction  of  the  signal  wavelength.  The  signal 
wavelengths  are  summarized  below: 


Table  1:  Signal  Wavelengths 


Signal 

Frequency 

Wavelength 

LI 

1575.42  MHz 

20  cm 

L2 

1227.60  MHz 

24  cm 

C/A  code 

IMHz 

300  m 

Pcode 

10  MHz 

30  m 
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^  Mfiasluaes  of  Acwracy 

There  are  two  measures  of  accuracy  that  are  commonly  used  to  describe  GPS  performance.  The 
military  commonly  uses  the  Circular  Error  Probable  (CEP)  in  two  dimensions,  and  the  Spherical 
Error  Probable  (SEP)  in  three  dimensions.  These  measures  are  defined  to  be  the  median  value  of 
accuracy.  Hence  50%  of  measurements  can  be  expected  to  fall  within  the  SEP  of  the  true  value  and 
50%  will  fall  outside. 

The  scientific  community  prefers  the  standard  deviation®^  as  a  measure  of  accuracy.  In  GPS 
circles,  the  2drms  value  or  twice  the  standard  deviation  is  often  used.  Between  95%  and  98%  of 
measurements  can  be  expected  to  fall  within  the  2drms  value  of  the  true  position,  depending  on  the 
ellipticity  of  the  distribution.  Both  measures  of  accuracy  are  related  by  the  approximate  formula: 


2  X  drms  »  2.5  x  CEP  Equation  40 


n  Selective  Availability 

The  nominal  precision  of  GPS  service  is  precisely  regulated  and  deliberately  controlled  by  the 
Department  of  Defense.  Two  levels  of  positiorung  service  are  distinguished  for  GPS.  The 
standard  positioning  service  (SPS)  is  designated  for  civilian  use  and  has  a  nominal  precision  of 
120  meters  2drms  in  3  dimensions.  SPS  is  based  on  reception  of  the  C/A  PRN  code.  Velocity 
measurement  is  intended  to  be  accurate  to  0.3  meters/second,  and  time  is  accurate  to  300-400 
nanoseconds. 

The  predse  positioning  service  (PPS)  is  designated  for  US  and  NATO  military  use  and  has  a 
nominal  precision  of  10  meters  2drms  in  three  dimensions.  PPS  is  based  on  the  simultaneous 
reception  of  the  encrypted  P  code  on  both  carriers.  Only  those  users  who  have  knowledge  of  the 
encryption  will  have  access  to  the  higher  PPS  precision. 

Through  a  technique  known  as  selective  availability  (S/A)  the  physically  achievable  40  meter 
precision  of  SPS  is  deliberately  degraded  to  the  120  meter  nominal  value.  This  is  achieved  by 
artificially  introducing  errors  into  the  broadcast  ephemeris  data,  and  by  dithering  the  satellite  clock. 
As  a  result,  S/A  errors  have  both  a  slowly  varying  and  a  rapidly  varying  component  For  a  fixed 
receiver,  measured  pseudorange  will  vary  apparently  rando^y  over  time  as  shown  in  Figure  31. 


67 .  The  use  of  standard  deviation  or  CEP  as  a  measure  of  “accuracy”  is  a  misnomer.  Bodi  measure  the  spread 
of  the  data,  not  the  differeoce  between  the  sanq)le  and  p(q>nlation  means.  CEP  and  2dnns  really  measure  ine- 
dsion.  Only  in  the  case  of  a  single  measurement  in  die  absence  of  measuremmit  bias  does  the  data  qnead 
idate  to  accnracy. 
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Figure  31  Selective  Availability 


H  Sottrccs  of  Error 

In  order  to  assess  the  possible  accuracies  available  in  an  application,  an  understanding  of  the 
sources  of  error  and  when  they  apply  is  necessary.  The  total  error  in  the  positioning  fix  can  be 
separated  into  two  components:  the  errors  in  the  satellite  pseudoranges,  and  the  projection  of  these 
errors  onto  the  navigation  coordinate  system.  The  first  component  involves  physics;  the  second  is 
purely  a  mathematical  transformation. 

Pseudorange  error  sources  are  additive.  The  ^ometric  dilution  of  precision,  or  GDOP  is 
multiplicative.  GDOP  can  be  as  high  as  4  or  much  more  under  certain  circumstances^,  so  it  is 
important  to  know  whether  precision  estimates  include  GDOP  or  not,  and  if  they  assume  a 
particular  value.  It  is  also  important  to  ascertain  whether  a  particular  precision  estimate  is 
horizontal  or  3D. 

3.4.1  Sources  of  Pseudorange  Error 

The  significant  sources  of  pseudorange  errors  and  their  nominal  values  are  summarized  in  Figure 
32 


68.  During  the  w*iriai  years  when  the  GPS  satellite  coostellatioD  is  incompletB,  OX>P  values  as  high  as  20 
are  possible. 
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Error  Source 

Nominal  Value  (rms) 

Selective  Availability 

8  meters 

Atmospheric  Delays 

4  meters 

Satellite  Clock  &  Ephemeris 

3  meters 

Multipath 

1-3  meters 

Reciever  Electronics  & 

Vehicle  Dynamics 

1.5  meters 

TOTAL 

10  meters 

Notes: 

a)  pseudorange  errors  only.GDOP  (of  4  -  6)  is  not  included 

b)  multiply  by  2  for  2dnns,  by  0.8  for  CEP 

c)  errors  combined  in  root  sum  square  sense 


Figure  32  Pseudorange  Error  Sources 

It  should  be  clear  from  this  table  that  the  GPS  design  reserves  the  ability  to  remove  the  two  most 
important  errors  for  exclusive  military  use.  With  the  exception  of  S/A,  the  deliberate  signal 
degradation  discussed  earlier,  atmospheric  delay  is  the  most  important  pseudorange  error  source. 
It  is  also  called  group  delay.  It  varies  with  time  and  place.  It  can  result  in  as  much  as  30  meters 
equivalmit  range  error  when  the  satellite  is  at  zenidi. 

Atmospheric  delay  consists  of  two  components:  ionospheric  and  tropospheric.  The  ionospheric 
delay  is  caused  by  diffraction  due  to  charged  particles.  It  can  vary  by  a  factor  of  S  from  day  to  night, 
being  highest  during  the  day,  and  by  a  factor  of  3  due  to  the  elevation  angle  of  the  satellite.  This 
means  t^t  each  satellite  has  a  different,  time  varying  pseudorange  error.  It  depends  on  solar 
magnetic  activity  and  geomagnetic  latitude,  being  particularly  affected  during  magnetic  storms  and 
it  is  usually  greatest  at  the  poles,  and  at  the  equator. 

Tropospheric  delay  is  caused  by  water  and  other  atmospheric  constituents  giving  rise  to  local 
changes  in  the  index  of  refraction.  It  can  be  modelled  very  well  by  a  simple  mathematical  formula. 
It  also  varies  with  time  and  place,  giving  about  2.3  meters  delay  vertically  and  ten  times  this  at  the 
horizon. 

Multipath  error  arises  when  some  of  the  radio  signal  reaches  the  receiver  through  reflection  in 
addition  to  the  portion  which  arrived  along  the  direct  line  of  sight.  The  result  of  both  signals 
arriving  at  the  antenna  is  destructive  mterference.  This  error  source  can  be  substantial  above  water. 
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since  water  is  a  perfect  radio  reflector  under  certain  conditions,  but  it  is  usually  less  than  1  meter. 
As  indicated  in  Figure  33,  the  time  delay  for  ground  reflection  depends  on  altitude.  Multipath  error 
is  more  pronounced  when  the  receiver  is  close  to  the  surface.  Therefore,  it  is  usuaUy  advisable  to 
mount  GPS  antennae  as  high  as  possible. 


Other  less  important  sources  of  error  include  the  unintentional  component  of  satellite  ephemeris 
and  clock  errors,  imperfections  in  receiver  circuitry,  computational  approximation  and  truncation, 
and  vehicle  dynamics.  The  effect  of  vehicle  motion  can  be  substantial  if  receiver  specifications  are 
exceeded. 

3.4,2  Geometric  Dilution  of  Precision 

The  final  GPS  positioning  accuracy  is  strongly  dependent  on  the  geometry  of  the  satellite  lines  of 
position  used.  The  transformation  from  four  pseudorange  errors  into  positioning  errors  is 
accomplished  by  multiplying  by  a  geometry  f^tor  known  as  the  geometric  dilution  of  precision. 
This  GDOP  occurs  in  aU  tiiangulation  problems  and  is  defined  as: 


DOP  - 

Equation  41 

RangeError 

where  the  errors  are  represented  as  standard  deviations.  The  DOP  scales  the  range  error  into  the  fix 
error,  so  small  values  indicate  favorable  geometry. 

The  concept  is  easy  to  visualize  in  two  dimensions.  Consider  Figure  34.  Range  errors,  shown  as 
shaded  regions,  combine  to  form  a  roughly  rectangular  region  of  uncertainty  in  the  fix.  The  greater 
the  angle  between  the  satellites  used,  the  smaller  is  the  region  of  uncertainty.  It  can  be  formally 
shown  that  there  is  a  direct  relationship  between  uncertainty  and  the  volume  of  this  region. 
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Figure  34  Geometric  Dilution  of  Precision 


It  the  context  of  GPS,  five  terms  are  defined: 

•  TDOP  -  time  dilution  of  precision  (range  equivalent) 

•  PDOP  -  position  dilution  of  precision  (3D) 

•  HDOP  -  horizontal  dilution  of  precision 

•  VDOP  -  vertical  dilution  of  precision 

•  GDOP  -  geometric  dilution  of  precision 

These  are  related  by: 


GDOP  =  , 

1  (PDOP)  ^  +  (TDOP)  ^  Equation  42 

PDOP  =  *1 

(HDOP)^+ (VDOP)^  Equation  43 

For  all  us^  averaged  over  the  globe,  and  averaged  over  time,  typical  values  of  PDOP,  HDOP,  and 
TDOP  are  2.6,  l.S,  and  1.2.  Good  receivers  will  compute  the  G1X)P  for  particular  combinations  of 
satellites  and  choose  the  best  set 
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Several  modes  of  use  are  available  in  mobile  robotics  applications.  All  applications  can  be 
classified  first  according  to  whether  the  GPS  receiver  is  used  as  the  real-time  motion  control 
sensor,  or  as  a  redundant  sensor  for  site  surveying  applications.  In  the  latter  case,  off-line  data 
reduction  potentially  combined  with  long  measurement  dwell  times  can  be  exploited  to  remove 
most  of  the  random  error  in  measurement.  This  is  a  highly  specialized  application,  so  the  reader  is 
referred  to  the  bibliography  for  further  reading.  This  section  will  consider  real-time  positioning 
accuracy  only. 

Receivers  can  compute  a  positioning  fix  at  a  maximum  rate  of  about  once  per  second,  so  this  is  the 
closest  to  real-time  performance  available.  This  implies  that  certain  highly  dynamic  applications 
will  require  at  least  a  rudimentary  dead  reckoning  system  to  augment  the  GPS  receiver.  Below  this 
maximum  update  rate,  the  whole  spectrum  of  averaging  times  can  be  used  to  eliminate  random 
error  provided  vehicle  dynamics  permit  it  However,  note  that  civilian  SPS  service  users  cannot 
remove  the  relatively  low  frequency  selective  availability  error  unless  dwell  times  on  the  order  of 
hours  are  used.  Only  the  higher  frequency  random  errors  can  be  removed  quickly  enough  for  real¬ 
time  positioning. 

Having  eliminated  long  dwell  times  from  consideration,  the  following  operational  modes  remain. 

3.5.1  Absolute  GPS 

This  mode  involves  simply  adopting  an  earth  centered,  earth  fixed  coordinate  system  for  high  level 
motion  control.  Receivers  report  position  in  this  system  most  naturally.  In  this  mode,  nominal  SPS 
performance  applies  to  civilian  users.  It  has  the  advantage  that  this  accuracy  applies  over  the  entire 
surface  of  the  earth. 

3.5.2  Relative  GPS 

It  may  seem  that  for  robotic  applications,  the  system’s  inherently  higher  relative  accuracy  can  be 
exploited  by  simply  adopting  die  initialization  point  as  the  origin  of  coordinates.  However,  this 
technique  can  only  remove  bias  that  is  independent  of  time  over  the  mission  duration,  and 
independent  of  position  over  the  mission  excursion.  S/A  errors  vary  slowly  over  time,  and  they  are 
the  dominant  component  of  error,  so  relative  accuracy  will  degrade  to  levels  close  to  the  absolute 
accuracy  over  a  time  period  of  only  a  few  minutes. 

3.53  Repetitive  GPS 

It  may  also  seem  that  the  system’s  inherent  repeatability  can  be  exploited  in  applications  where  the 
vehicle  can  simply  repeat  the  acquisition  of  previously  surveyed  waypoints.  Again,  the  selective 
availability  error  limits  the  improvement  in  accuracy  achievable  with  this  approach  because  GPS 
will  report  significantly  different  positions  for  the  same  physical  point  on  the  earth  over  a  period 
of  just  a  few  minutes. 
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3.5.4  Coded  Differential  GPS 


This  technique  uses  a  second  fixed  receiver  as  the  origin  of  coordinates.  If  the  earth  centered,  earth 
fixed  location  of  this  reference  is  known  to  high  accuracy,  then  high  accuracy  absolute  positioning 
of  the  moving  receiver  is  possible.  If  not,  ^en  only  high  accuracy  positioning  relative  to  the 
reference  is  possible.  This  second  case  is  sufficient  for  most  robotics  applications.  This  mode  of 
GPS  use  is  particularly  important,  so  some  space  will  be  devoted  to  acquainting  the  reader  with  the 
issues. 

Since  the  position  of  the  reference  receiver  is  known  in  the  navigation  coordinate  system,  it  can 
compute  die  difference  between  this  and  the  GPS  reported  position  in  real  time.  Differential 
techniques  rely  on  a  second  communications  channel  which  allows  the  reference  to  broadcast  these 
differences  to  the  moving  receiver  at  a  high  rate.  If  omnidirectional  radio  is  used  to  transmit  the 
position  error,  many  mobile  receivers  can  benefit  from  a  single  reference^^.  Tethered  vehicles  can 
transmit  differential  corrections  through  the  tether.  The  general  situation  is  indicated  in  Figure  35. 


In  implementation  it  is  not  advisable  to  transmit  cartesian  positioning  error  from  the  reference  since 
the  moving  sets  may  use  a  different  navigation  solution,  or  may  even  use  a  different  group  of 
satellites.  Rather,  because  the  reference  has  no  knowledge  of  which  satellites  will  be  used,  it  must 
compute  the  pseudorange  error  for  all  satellites  in  view  and  transmit  these  to  the  moving  receivers. 


69 .  Local  oomineicial  differentiai  GPS  transmiogs  will  be  operating  in  urban  centers  soon.  The  standard  pro¬ 
tocol  for  differential  service  is  known  as  RTCM-SC-104. 
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Differential  techniques  can  remove  only  that  component  of  error  which  is  common  to  both  the 
reference  and  the  moving  set.  Pseudorange  error  varies  with  both  position  at  a  fixed  time  and  with 
time  at  a  fixed  position.  Position  variation  is  primarily  due  to  the  dependence  of  atmospheric  delay 
and  ephemeiis  errors  on  elevation  angle.  Luckily,  this  geometric  decorrelation  is  relatively 
insignificant  because  of  the  small  change  in  elevation  angle  associated  with  moderate  excursion 
from  the  reference.  This  is  a  direct  result  of  the  substantial  altitude  of  the  spacecraft  as  indicated  in 
Figure  36.  The  remaining  uncancelled  error  is  roughly  linear  in  distance.  Useful  differential 
pseudorange  corrections  can  be  obtained  over  distances  as  large  as  1,000  kilometers. 


Figure  36  Spatial  Variation  of  System  Error 

While  spatial  variation  is  insignificant,  variation  over  time  is  a  major  concern.  The  finite  time 
required  to  transmit  and  process  pseudorange  errors  imply  that  the  moving  receiver  will  always 
have  a  slightly  dated  estimate  available.  The  S/A  range  error  rate  is  such  that  pseudorange  error 
data  is  useless  after  one  minute  of  elapsed  time.  Accordingly,  the  update  rate  of  the  secondary  data 
link  is  a  critical  parameter.  For  a  stationary  receiver  near  the  reference,  2drms  error  will  vary  over 
time  as  indicated  in  Figure  37. 
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Figure  37  Temporal  Variation  of  System  Error 

3.5.5  Codeless  Djfferential  GPS 

This  technique  relies  on  bypassing  the  PRN  codes  entirely  and  using  only  the  high  resolution 
carriers.  Differential  carrier  phase  tracking  consists  of  measuring  the  phase  shift  between  the 
same  signal  received  at  two  different  locations.  The  phase  shift  is  a  direct  measure  of  the  difference 
in  path  length  from  the  satellite  to  each  antenna.  Also  known  as  interferometric  GPS,  the 
principle  is  indicated  in  Figure  38. 


Such  measurements  are  subject  to  the  phase  ambiguity  problem  and  special  techniques  are  used  to 
resolve  it  Also,  since  a  continuous  lock  must  be  kept  on  the  carriers  at  all  times,  any  interruption 
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of  the  line  of  sight  requires  that  the  system  be  reinitialized.  Receivers  may  sometimes  miss  a  carrier 
cycle  in  their  computations,  and  special  techniques  are  employed  to  eliminate  this  cycle  slip 
problem.  Accuracies  increase  by  an  order  of  magnitude  at  the  cost  of  limited  excursion,  higher  cost 
equipment,  and  sometimes  reduced  update  rate.  Relative  accuracies  of  a  few  parts  per  million  of 
excursion  are  the  rule.  For  systems  measuring  rotation  on  a  rigid  vehicle,  1  milliradian  of  accuracy 
is  achievable  over  a  baseline  of  1  meter. 

The  advantages  and  disadvantages  of  each  mode  of  operation  are  summarized  in  Figure  39 


Mode 

Accuracy 

Range 

Why 

Absolute 

120  m 

global 

Relative 

100  m 

local 

remove  position  and 
time  independent  bias 

Repetitive 

100  m 

global 

removes  time 
independent  bias 

Coded  Differential 

3-5  m 

20  Km 

remove  position  indep 
and  time  dependent  bias 

Codeless  Differential 

2  ppm 

10  Km 

bypasses  codes  and 
uses  interferometry 

Notes; 

a)  All  coded  modes  using  C/A  code,  update  rate  ~  1  Hz 

b)  Bias  is  error  varying  slower  than  1  Hz 


Figure  39  GPS  Modes  of  Operation 


The  generic  GPS  receiver  consists  of  an  antenna  and  associated  preamplifier,  an  RF  section, 
tracking  loops,  a  microprocessor,  and  optionally  includes  its  own  power  supply,  a  command  and 
display  unit,  and  an  external  data  and  control  port.  Receivers  exist  which  are  designed  to  suit  a 
number  of  diverse  applications  including  aviation,  land  vehicles,  marine,  military,  surveying, 
timing,  and  even  handheld  varieties.  Some  of  the  important  features  of  commercial  sets  are 
summarized  below. 

3.6.1  Application 

For  robotics  applications,  a  unit  which  is  designed  as  a  component  of  a  larger  navigation  system  is 
appropriate.  This  often  implies  facilities  for  external  power,  communication  with  an  external 
computer,  an  external  antenna,  and  significant  internal  storage  may  or  may  not  be  required. 
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3.6Jt  Configuration 


Some  units  are  available  for  stand-alone  use;  others  are  configured  in  a  printed  circuit  board  set. 
There  are  also  many  complete  navigation  systems  available  incorporating  integrated  GPS  and 
inertial  navigation  systems.  Volume,  weight,  power  consumption,  port  baud  rates,  and 
environmental  hardness  can  vary  significantly. 

3.6.3  InputiOutput 

Some  sets  report  position,  velocity,  and  time  and  some  report  only  position.  Special  purpose  multi 
antenna  sets  can  measure  orientation.  It  may  be  necessary  to  accept  differential  range  corrections 
according  to  a  standard  protocol.  Some  differential  sets  are  intended  for  postprocessed  surveying 
applications  only  and  postprocessing  software  may  or  may  not  be  supplied. 

Of  particular  importance  is  the  coordinate  systems  in  which  position  can  be  reported.  Options 
include  cartesian  earth-centered  coordinates;  latitude,  longitude,  and  altitude  on  a  reference 
ellipsoid;  easting,  northing  and  altitude;  and  range,  bearing,  speed,  and  crosstrack  error  with 
respect  to  a  stored  sequence  of  waypoints. 

3.6.4  Antenna 

Several  types  of  antenna  are  available  including  monopole  or  dipole  configurations,  quadrifilar 
helices,  spiral  helices,  and  microstrips.  Provision  for  an  external  antenna  is  not  always  available. 
For  precision  applications,  an  important  parameter  is  the  stability  of  the  anterma  phase  center,  the 
electrical  center  of  the  antenna  to  which  the  receiver  output  actually  refers. 

3.6.5  Channels 

Most  receivers  will  operate  solely  on  the  C/A  PRN  code,  and  the  number  of  channels  of  radio  input 
can  vary  significantly.  Some  sets  integrate  GPS  and  Glonass  tracking.  Precision  codeless  sets  use 
orrly  the  carrier  signal.  A  set  with  several  channels  can  track  several  satellites  simultaneously, 
rather  than  sequentially,  and  since  it  will  exhibit  better  signal  to  noise  ratios,  it  can  track  them 
closer  to  the  horizon. 

3.6.6  Update  Rate 

Update  rates  and  the  acceptable  vehicle  velocity  and  acceleration  increase  with  the  number  of 
chaimels.  Some  sets  require  an  initial  estimate  of  position  to  start,  and  others  can  cold  start.  The 
time  to  the  first  fix  can  vary  from  seconds  to  30  minutes  depending  on  the  initial  information 
supplied  and  the  receiver  mode  of  operation. 
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4  Glossary 

itiwftial  systems  -  inertial  systems  which  make  use  of  redundant  external  measurements  of  navigation  quantities. 
alongtradc  -  the  direction  parallel  to  the  direction  of  motion. 

analvrir .  a  designation  for  systems  which  compute  rather  than  instrument  the  navigation  frame  quantities. 

ha.se  point  -  designation  for  system  which  mainrain  the  platform  orientation  at  the  initial  orientation  widi  respect  to  the 
earth. 

crosstrack  -  the  direction  perpendicular  to  the  direction  of  motion,  measured  horizontally. 

<hminii^y  .  the  process  of  filtering  redundant  external  measurements  with  the  INS  output  in  order  to  reduce  the 
magnitude  of  system  error. 

directicMi  cosine  mairiT .  the  mafrix  which  (Ascribes  the  transformation  of  the  body  axes  into  the  navigatirm  frame, 
earth  fixed -see  base  point. 

eartii-stabilized  -  designation  for  systems  which  maintain  anavigation  frame  which  is  defined  with  respect  to  die  eardi 
in  some  way. 

easting  -  eastern  coordinate  in  UTM  coordinates. 

free  azimuth  -  a  designation  for  systems  providing  no  platfram  rate  about  the  azimuth. 

ge<an«ric  -  designation  for  systems  which  perform  geometric  stabilization  utilizing  stable  platform. 

gimballed  •  designation  to  systems  incorporating  a  stable  platform. 

gimballodc  -  singularity  of  the  platform  gimbal  mechanism. 

gravity  •  vector  sum  of  gravitation  and  centrifugal  force. 

gnd  -  a  regular  coordinate  grid  placed  on  a  map.  Unlilm  latitude  and  longitude,  grid  lines  are  parallel  and  equally 
placed. 

gyrocompas-sing  •  the  {uocess  of  using  the  earth’s  spin  to  find  die  direction  of  north. 

inertial  measurement  unit  -  the  sensor  sirite  consisting  of  the  gyros  and  accelerometers. 

inertial  mfwence  unit  -  see  inertial  mea.«aiiement  imit 

laritiide .  angle  of  a  point  on  the  earth  measured  from  the  equatorial  plane. 

levelliny .  the  process  of  levelling  a  stable  plattom  during  alignment. 

inraiiv  level .  desigoadon  for  systems  which  maintain  the  platform  level. 

local  vertical  -  see  locally  level. 

longitnde  -  angle  of  a  point  on  the  earth  measured  frran  the  prime  meridian. 

-  the  innx%ss  by  whidi  the  navigation  equadons  are  implemented  in  hardware  and  software. 
m*ir«dian  •  a  Une  d  constant  longinide. 
n«w«biny .  northern  coordinate  in  UTM  cotudinates. 

north-slaved  -  designation  for  systems  which  physically  trade  the  direction  of  north  by  rotating  the  platform, 
parallel  -  a  Hue  of  ennstanf  laritiide. 

OKkoS  -  the  transducer  wbidi  measures  a  gimbal  angle. 

nnle  nmhiein .  the  d^radadtm  of  petfonnanoe  e:q)erienced  by  inertial  systems  as  they  apisoach  the  eaith’ s  poles, 
rehatanre.  loops  -  a  contiol  loop  wfaich  actively  prohibits  deflectimi  of  a  compliaa  t  transducer. 

refeMioe  elllmnid  .  an  lA^aHawd  marhemartcal  mndel  of  the  earth 
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Schuler  loop  •  the  loop  which  enq>ioys  gravity  feedback  in  order  to  level  a  table. 

Schuler  period  >  die  period  of  oscillation  of  the  Schuler  pendulum,  or  84  minutes. 

Schuler  tuning  <  the  fmxxss  of  designing  a  system  to  exhibit  the  Schuler  namial  frequency, 
self  alignment .  the  process  of  aligning  an  INS  without  the  aid  of  external  measurements, 
semi-analvrir  -  a  designation  for  systems  which  instrument  only  the  navigadon  frame, 
signal  geMtator  -  see  pickoff . 

space-stahilized  -  designation  for  systems  which  maintain  constant  inertial  platfonn  orientation. 
specific  force  -  the  net  inertial  acceleration  minus  gravitational  attractions. 

■Stable  nlatfftrm .  a  platform  whose  orientation  is  stable  in  a  particular  reference  frame. 
strapdown  -  see  analytic. 

Star-tracker  -  a  device  which  tracks  ad  reports  the  beating  of  a  star/ 
tocQuer  -  device  for  applying  torque  to  a  gyro  or  gimbal. 

vehicle  rate  -  the  angular  velocity  of  the  vehicle  with  respect  to  the  center  of  the  earth, 
wander  angle  -  the  rotation  of  the  coordinate  axes  with  respect  to  north. 

wander  azimiirti .  a  designation  for  a  system  which  tcnques  the  azimudi  to  follow  the  earth  rate  only. 

almaoac  -  imprecise  satellite  qdiemeris  data  broadcast  to  allow  receivers  to  predict  when  .sateHites  will  rise  above  the 
horizon. 

baseline  •  the  distance  between  two  antomae  in  differential  GPS. 

codeless  •  a  designation  for  receivers  which  compute  position  firtxn  only  die  signal  carriers. 

cycle  slip  -  a  temporary  loss  of  lodt  in  the  carrier  tracking  loop  of  a  receiver. 

differeniial  carrier  pha.se  fraclring  .  a  technique  which  improves  accuracy  by  measuring  the  carrier  jdiase  shift. 

enhemeris  data  -  the  data  specifying  die  orbital  position  nf  a  satedite 

geometric  demrrRistfinn .  the  ttegree  to  nhich  system  error  is  a  noounifonn  function  of  position. 

geometric  dilution  of  luecision  -  the  dilution  of  precision  which  occurs  when  pseudoranges  are  transftmned  into 
cartesian  positKKi. 

GPS  svsrem  time .  the  time  Standard  to  which  all  satellites  are  actively  synchronized  by  MCS. 
group  delay  -  the  pnqiagatkm  delay  when  radio  passes  through  the  atmoqihete. 

hmdnverwqed .  a  ctmqxmeQt  of  the  data  message  which  facilitates  the  transfer  of  tracking  from  C/A  code  to  P  code. 
hiierfiMtimetric  ca*S  -  see  differential  carrier  phase  traddng. 

™wif  anyii» .  the  elevatioQ  angle  to  the  horizon  above  which  a  satellite  is  practically  visible. 

naviyatifui  mesaape.  -  a  distal  data  message  modulated  onto  both  carriers  which  provides  epbemeris  data,  satellite 
clock  offsets,  and  many  other  usefril  pieces  of  information. 

rmtPT .  the  point  on  the  antenna  whidi  is  actually  being  positioned  by  the  recdver. 

nrecige  noriiioBing  .warvity .  the  most  precise  dynamic  positioning  service  available  with  GPS. 

psendotange  -  saielliiB  to  receiver  range  errms  whidi  indude  the  user  dodc  Inas. 

ESl!i.£f]da  *  digittl  codes  which  have  characteristics  shnilar  to  noise. 

MtterriWL  availflhilify  .  rtv> 

SUpDCaLdax '  ti>e  time  required  for  the  eartii  to  rotate  once  with  respect  to  the  stars.  It  is  shorter  than  die  mean  solar 
day  by  four  minutes. 
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stendairi  positioning  service  -  tie  less  accurate  posidoniag  service  available  to  civilians. 
space  loss  -  tbe  absorption  of  radio  energy  as  it  passes  through  the  atmosphere. 
user  clodc  bias  -  the  offset  of  a  receiver’s  clock  from  GPS  system  time. 
zenith  -  directly  overhead. 
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